Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 150,537 A→C 2.3% D139E (GAT→GAG yadC ← fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_007779150,5370A→C 2.3% 3.4 ‑2.1 0.0327D139E (GAT→GAGyadCfimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  new base (9/2):  ref base (180/136):  total (189/138)
Fisher's exact test strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

ACCCTTAGTTTGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACCTGAATAGA  >  NC_007779/150438‑150637
                                                                                                   |                                                                                                    
aCCCTTAGTTTGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTAt                                                                                                     <  1:1226425/101‑1 (MQ=255)
 cccTTAGTTTGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATc                                                                                                    >  1:3245795/1‑101 (MQ=255)
   cTTAGTTTGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGt                                                                                                  <  1:12249519/101‑1 (MQ=255)
    ttAGTTTGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTa                                                                                                 <  1:1991788/101‑1 (MQ=255)
     tAGTTTGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTAt                                                                                                <  1:7908962/101‑1 (MQ=255)
     tAGTTTGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTAt                                                                                                <  1:8752986/101‑1 (MQ=255)
      aGTTTGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATg                                                                                               >  1:2642115/1‑101 (MQ=255)
       gTTTGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGc                                                                                              <  1:1187668/101‑1 (MQ=255)
        tttGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGca                                                                                             <  1:10161132/101‑1 (MQ=255)
         ttGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGcag                                                                                            <  1:1271103/101‑1 (MQ=255)
          tGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGcagt                                                                                           >  1:10614267/1‑100 (MQ=255)
              gATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGcagaccag                                                                                       >  1:13349054/1‑101 (MQ=255)
              gATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGcagaccag                                                                                       >  1:6105557/1‑101 (MQ=255)
              gATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTATCTCGTATGccgtcttc                                                                                       >  1:1834536/1‑94 (MQ=255)
              gATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTATCTCGTATGccgtcttc                                                                                       >  1:9938003/1‑94 (MQ=255)
               aTCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGcagaccaga                                                                                      <  1:8934940/101‑1 (MQ=255)
               aTCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGcagaccaga                                                                                      >  1:2878763/1‑101 (MQ=255)
               aTCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGcagaccaga                                                                                      >  1:10428360/1‑101 (MQ=255)
               aTCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGcagaccaga                                                                                      >  1:4721366/1‑101 (MQ=255)
                tCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGcagaccagac                                                                                     <  1:9027356/101‑1 (MQ=255)
                tCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGcagaccagac                                                                                     <  1:9382917/101‑1 (MQ=255)
                  gCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGAccc                                                                                   >  1:11174118/1‑101 (MQ=255)
                  gCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGAccc                                                                                   <  1:5664270/101‑1 (MQ=255)
                  gCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCACAccc                                                                                   >  1:11311959/1‑101 (MQ=255)
                   ctgtctgtGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCt                                                                                  <  1:2556461/101‑1 (MQ=255)
                   ctgtctgtGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCt                                                                                  >  1:11690665/1‑101 (MQ=255)
                    tgtctgtGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCtt                                                                                 >  1:3062314/1‑101 (MQ=255)
                     gtctgtGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTg                                                                                >  1:5777711/1‑101 (MQ=255)
                     gtctgtGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTg                                                                                <  1:6270745/101‑1 (MQ=255)
                      tctgtGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGa                                                                               <  1:9684499/101‑1 (MQ=255)
                        tgtgACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGata                                                                             <  1:11051515/101‑1 (MQ=255)
                         gtgACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGatat                                                                            >  1:981580/1‑101 (MQ=255)
                            aCGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTaa                                                                         <  1:179915/101‑1 (MQ=255)
                             cGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAAc                                                                        >  1:4245501/1‑101 (MQ=255)
                               gaaaagaaaaATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACAt                                                                      >  1:3774168/1‑101 (MQ=255)
                                     aaaaaTTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTa                                                                >  1:12892191/1‑101 (MQ=255)
                                      aaaaTTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTaa                                                               >  1:476749/1‑101 (MQ=255)
                                      aaaaTTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTaa                                                               <  1:7994345/101‑1 (MQ=255)
                                         aTTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTGACATGGTGTAATAc                                                            >  1:10121951/1‑101 (MQ=255)
                                            cTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGa                                                         <  1:11808766/101‑1 (MQ=255)
                                            cTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGa                                                         <  1:6685434/101‑1 (MQ=255)
                                             ttggttggAAGGATCTCTGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGAc                                                        >  1:13056774/1‑101 (MQ=255)
                                              tggttggAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTAGCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGAcc                                                       >  1:8541113/1‑101 (MQ=255)
                                               ggttggAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCa                                                      >  1:12620769/1‑101 (MQ=255)
                                                 ttggAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAgg                                                    <  1:5992796/101‑1 (MQ=255)
                                                 ttggAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTATATCGTATGCCGTCTTCTGCTTGAAAAAAAACATGGTGTAATACAGACCAgg                                                    >  1:7485380/1‑101 (MQ=255)
                                                  tggAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTTTAATACAGACCAGGa                                                   >  1:9481587/1‑101 (MQ=255)
                                                   ggAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTATCTCGTATGCCGTCTTCTGCTTTGAAATTAACATGGTGTAATACAGACCAGGaa                                                  >  1:7669062/1‑101 (MQ=255)
                                                    gAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAAc                                                 <  1:6353947/101‑1 (MQ=255)
                                                    gAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAAc                                                 >  1:13583529/1‑101 (MQ=255)
                                                     aaGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCACACCCTTGATATTAACATGGTGTAATACAGACCAGGAACa                                                >  1:359415/1‑101 (MQ=255)
                                                        gATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGAt                                             >  1:7127938/1‑101 (MQ=255)
                                                         aTCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATg                                            <  1:5357464/101‑1 (MQ=255)
                                                           ctcCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTa                                          <  1:2245764/101‑1 (MQ=255)
                                                           ctcCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTa                                          <  1:8079136/101‑1 (MQ=255)
                                                           ctcCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTa                                          <  1:7260298/101‑1 (MQ=255)
                                                            tcCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTAt                                         >  1:3912302/1‑101 (MQ=255)
                                                            tcCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTAt                                         >  1:12464034/1‑101 (MQ=255)
                                                            tcCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTAt                                         <  1:2697336/101‑1 (MQ=255)
                                                            tcCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTAt                                         >  1:2274970/1‑101 (MQ=255)
                                                              cGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTa                                       >  1:4828633/1‑101 (MQ=255)
                                                               gATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTaa                                      >  1:9544044/1‑101 (MQ=255)
                                                               gATATAGATTCCTGGCGATTGAATGTCAGTTATTATCTCGTATGCCGTCTTCTGCTTTGAAATTAACATGGTGTAATACAGACCAGGAACAGATGTATTaa                                      >  1:1381606/1‑101 (MQ=255)
                                                                atatAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTaaa                                     >  1:662341/1‑101 (MQ=255)
                                                                  ataGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATa                                   >  1:10270778/1‑101 (MQ=255)
                                                                  ataGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATa                                   >  1:9078069/1‑101 (MQ=255)
                                                                     gATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGtt                                <  1:6520431/101‑1 (MQ=255)
                                                                      aTTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGttt                               <  1:1158684/101‑1 (MQ=255)
                                                                      aTTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGttt                               >  1:3011903/1‑101 (MQ=255)
                                                                       ttCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTa                              <  1:6593842/101‑1 (MQ=255)
                                                                       ttCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTa                              <  1:5504475/101‑1 (MQ=255)
                                                                         ccTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATg                            >  1:7341896/1‑101 (MQ=255)
                                                                         ccTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATg                            >  1:3472521/1‑101 (MQ=255)
                                                                           tGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTc                          <  1:4333052/101‑1 (MQ=255)
                                                                             gCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCg                        <  1:2220931/101‑1 (MQ=255)
                                                                              cgcttGAATGTCAGTTATTGTATCGTATGCAGACCACACCCTTGATATTAACATGGTGTAATACAGACCAGGACCAGATGTATTAAATAGTTTATGTCCGc                       >  1:13226783/4‑101 (MQ=255)
                                                                              cGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGc                       >  1:12790403/1‑101 (MQ=255)
                                                                               gATTGAATGTCAGTTATTGTATCGTATGCAGACCTGACCTTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGcc                      >  1:8199912/1‑101 (MQ=255)
                                                                               gATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGcc                      >  1:13188207/1‑101 (MQ=255)
                                                                                aTTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCa                     >  1:5255129/1‑101 (MQ=255)
                                                                                 ttGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCAt                    <  1:154194/101‑1 (MQ=255)
                                                                                  tGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATa                   >  1:82741/1‑101 (MQ=255)
                                                                                    aaTGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAAt                 <  1:10819657/101‑1 (MQ=255)
                                                                                    aaTGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAAt                 >  1:2384669/1‑101 (MQ=255)
                                                                                    aaTGTCAGTTATTGTATCGTATGCAGACCAGACCCTCGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAAt                 <  1:2813142/101‑1 (MQ=255)
                                                                                     aTGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATc                >  1:7414113/1‑101 (MQ=255)
                                                                                        tCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCttt             >  1:10493955/1‑101 (MQ=255)
                                                                                        tCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCttt             >  1:5604059/1‑101 (MQ=255)
                                                                                         cAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTa            >  1:291433/1‑101 (MQ=255)
                                                                                         cAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTa            <  1:9151886/101‑1 (MQ=255)
                                                                                          aGTTATTGTCTCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTAc           <  1:11369431/101‑1 (MQ=255)
                                                                                          aGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTAc           >  1:8655067/1‑101 (MQ=255)
                                                                                           gTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTAcc          <  1:7792034/101‑1 (MQ=255)
                                                                                           gTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTAcc          >  1:8842475/1‑101 (MQ=255)
                                                                                             tATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACCTg        >  1:4256449/1‑101 (MQ=255)
                                                                                              aTTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACCTGa       <  1:3120250/101‑1 (MQ=255)
                                                                                                tGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACCTGAAt     >  1:1399111/1‑101 (MQ=255)
                                                                                                  tATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACCTGAATa  >  1:1653792/1‑100 (MQ=255)
                                                                                                  tATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACCTGAATAg   <  1:1596852/101‑1 (MQ=255)
                                                                                                  tATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGAACAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACGTGAATAg   >  1:5105407/1‑101 (MQ=255)
                                                                                                   |                                                                                                    
ACCCTTAGTTTGTAGATCGCTGTCTGTGACGGAAAAGAAAAATTCTTGGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACCTGAATAGA  >  NC_007779/150438‑150637

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ 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.