Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,451,708 C→A 34.6% M1243I (ATG→ATT)  rpoB ← RNA polymerase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,451,7080C→A34.6% 0.0 ‑0.2 ‑0.0358M1243I (ATG→ATT) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  new base (61/63):  ref base (112/122):  total (173/185)
Fisher's exact test strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.74e-01

Only 100 of 363 total aligned reads displayed.

CCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGG  >  NC_007779/3451608‑3451807
                                                                                                    |                                                                                                   
cccACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGc                                                                                                     <  1:9184829/101‑1 (MQ=255)
 ccACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCa                                                                                                    <  1:14317221/101‑1 (MQ=255)
  cacaCCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAAt                                                                                                   >  1:5081795/1‑101 (MQ=255)
   acacCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATc                                                                                                  >  1:11497966/1‑101 (MQ=255)
    cacTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCt                                                                                                 <  1:11333998/101‑1 (MQ=255)
    cacCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCt                                                                                                 >  1:10585393/1‑101 (MQ=255)
      cctCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCttc                                                                                               <  1:2676973/99‑2 (MQ=255)
       ttCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgt                                                                                              <  1:8488960/101‑1 (MQ=255)
       ttCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgt                                                                                              >  1:11284311/1‑101 (MQ=255)
       ttCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgt                                                                                              <  1:11877164/101‑1 (MQ=255)
        tCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtc                                                                                             <  1:1817758/101‑1 (MQ=255)
         ccATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtcg                                                                                            <  1:373249/101‑1 (MQ=255)
          cATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtcgt                                                                                           >  1:6698479/1‑101 (MQ=255)
          cATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTctcgt                                                                                           >  1:15138947/1‑101 (MQ=255)
           aTCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTctcgtc                                                                                          <  1:3249457/101‑1 (MQ=255)
            tctcCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGtcgtcg                                                                                         >  1:14567297/1‑101 (MQ=255)
            tctcCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGtcgtcg                                                                                         <  1:13023328/101‑1 (MQ=255)
             ctcCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGa                                                                                        <  1:292462/101‑1 (MQ=255)
             ctcACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGTTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGa                                                                                        <  1:542765/101‑1 (MQ=255)
             ctcACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGa                                                                                        >  1:11143052/1‑101 (MQ=255)
               ccccGAAACGCTGACCACCGAACTATGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGAcc                                                                                      <  1:8849716/101‑1 (MQ=255)
                accGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCa                                                                                     >  1:14998513/2‑101 (MQ=255)
                 ccGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAg                                                                                    <  1:12120914/101‑1 (MQ=255)
                   gAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAggt                                                                                  <  1:1670121/101‑1 (MQ=255)
                     agcGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAggtgg                                                                                >  1:717607/3‑101 (MQ=255)
                     aaCGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtgg                                                                                <  1:8016921/101‑1 (MQ=255)
                      aCGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtggt                                                                               <  1:7017468/101‑1 (MQ=255)
                       cGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGtt                                                                              <  1:2469482/101‑1 (MQ=255)
                       cGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAggtggtg                                                                              <  1:8405793/101‑2 (MQ=255)
                        gCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAggtggtgc                                                                             <  1:8990584/101‑3 (MQ=255)
                         cTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCa                                                                            >  1:10191905/1‑101 (MQ=255)
                         cTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCa                                                                            >  1:9372232/1‑101 (MQ=255)
                          tGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAg                                                                           >  1:3362655/1‑101 (MQ=255)
                          tGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAg                                                                           <  1:3661027/101‑1 (MQ=255)
                           gACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGt                                                                          <  1:7076978/101‑1 (MQ=255)
                            accaccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGtt                                                                         >  1:8586959/1‑101 (MQ=255)
                             ccaccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGttt                                                                        <  1:14892323/101‑1 (MQ=255)
                             ccaccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGctt                                                                        <  1:6443706/101‑1 (MQ=255)
                              caccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACTAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCCGGTGGTTCAGTTTc                                                                       <  1:10271417/101‑1 (MQ=255)
                              caccGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTc                                                                       <  1:6381281/101‑1 (MQ=255)
                                ccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAg                                                                     >  1:12259282/1‑101 (MQ=255)
                                 cGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGc                                                                    >  1:9436878/1‑101 (MQ=255)
                                  gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                                                   >  1:8392202/1‑101 (MQ=255)
                                  gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                                                   <  1:13156732/101‑1 (MQ=255)
                                  gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                                                   >  1:6329484/1‑101 (MQ=255)
                                   aaCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcat                                                                  >  1:6243924/1‑101 (MQ=255)
                                    aCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatg                                                                 <  1:14824831/101‑1 (MQ=255)
                                      tgtgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGcatgta                                                               >  1:13052361/1‑101 (MQ=255)
                                      tgggCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgta                                                               >  1:9616076/4‑101 (MQ=255)
                                          ccTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacatg                                                           <  1:7827935/101‑1 (MQ=255)
                                            ttACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacatgta                                                         >  1:6116082/1‑101 (MQ=255)
                                            ttACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacatgta                                                         <  1:9836264/101‑1 (MQ=255)
                                             tACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtacatgtag                                                        >  1:9783810/1‑100 (MQ=255)
                                             tACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTaa                                                        >  1:13573798/1‑101 (MQ=255)
                                              accaccCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAAc                                                       >  1:14172669/1‑101 (MQ=255)
                                                caccCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGccc                                                     >  1:7676872/1‑100 (MQ=255)
                                                caccCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCa                                                     >  1:4883551/1‑101 (MQ=255)
                                                 accCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCaa                                                    <  1:7393966/101‑1 (MQ=255)
                                                  cccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGATCGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAAc                                                   <  1:10884200/101‑1 (MQ=255)
                                                  cccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAAc                                                   <  1:10837943/101‑1 (MQ=255)
                                                   ccAGCGGCTGCTGGGTAATCATGGAGTACGGGCCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAcc                                                  >  1:8105707/1‑100 (MQ=255)
                                                   ccAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAcc                                                  >  1:5167211/1‑100 (MQ=255)
                                                   ccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACg                                                  >  1:1285489/1‑101 (MQ=255)
                                                     aGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccga                                                >  1:1265599/1‑98 (MQ=255)
                                                       cGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTa                                              <  1:14703304/101‑1 (MQ=255)
                                                        ggCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTAc                                             >  1:4037800/1‑101 (MQ=255)
                                                           tgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCgg                                          >  1:1609854/1‑101 (MQ=255)
                                                              tgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACg                                       >  1:10421470/1‑101 (MQ=255)
                                                                  tAACCATGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg                                   >  1:8727979/1‑101 (MQ=255)
                                                                  tAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg                                   >  1:6495134/1‑101 (MQ=255)
                                                                   aaCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGa                                  <  1:12786146/101‑1 (MQ=255)
                                                                   aaCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGa                                  >  1:11665858/1‑101 (MQ=255)
                                                                     ccAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAAc                                <  1:4317000/101‑1 (MQ=255)
                                                                     ccAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAAc                                <  1:6045920/101‑1 (MQ=255)
                                                                          cTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTc                           >  1:69153/1‑101 (MQ=255)
                                                                           tGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCa                          >  1:1868423/1‑101 (MQ=255)
                                                                             tAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCAcc                        >  1:3093039/1‑101 (MQ=255)
                                                                              aaGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCa                       >  1:13699774/1‑101 (MQ=255)
                                                                               aGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAg                      <  1:8288600/101‑1 (MQ=255)
                                                                                gAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGt                     <  1:2756469/101‑1 (MQ=255)
                                                                                gAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGt                     >  1:12012881/1‑101 (MQ=255)
                                                                                gAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGt                     >  1:2205089/1‑101 (MQ=255)
                                                                                    cGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGa                 >  1:7247919/1‑101 (MQ=255)
                                                                                    cGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGa                 <  1:5935384/101‑1 (MQ=255)
                                                                                       tGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCa              <  1:5642363/101‑1 (MQ=255)
                                                                                       tGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGccg              >  1:9937658/1‑100 (MQ=255)
                                                                                       tGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGccg              >  1:4212530/1‑100 (MQ=255)
                                                                                        ggAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCAt             >  1:4827013/1‑101 (MQ=255)
                                                                                        ggAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCAt             <  1:13936059/101‑1 (MQ=255)
                                                                                        ggAACGAGCGTCAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGccgt             >  1:5041902/1‑99 (MQ=255)
                                                                                         gAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATc            >  1:10770124/1‑101 (MQ=255)
                                                                                         gAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATc            >  1:1897187/1‑101 (MQ=255)
                                                                                         gAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATc            <  1:1559722/101‑1 (MQ=255)
                                                                                           aCGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCga          >  1:4130545/1‑100 (MQ=255)
                                                                                             gcgcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTAc        >  1:12386711/1‑101 (MQ=255)
                                                                                              cgcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACa       >  1:14621074/1‑101 (MQ=255)
                                                                                               gcgTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAg      <  1:13284728/101‑1 (MQ=255)
                                                                                               gcgTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCACCGAAACCGGGAACGGGAACTTTTCGCCAGAGCGGCCGTCGAACAg      >  1:1926835/1‑101 (MQ=255)
                                                                                                 ttGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGc    >  1:4599958/2‑101 (MQ=255)
                                                                                                   gCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCgg  <  1:441455/101‑1 (MQ=255)
                                                                                                    |                                                                                                   
CCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGTTCACCAGTGCGACCATCGTACAGGCGG  >  NC_007779/3451608‑3451807

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.