Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,134,668 T→A 6.8% intergenic (+111/‑61) flgF → / → flgG flagellar basal‑body rod protein FlgF/flagellar basal‑body rod protein FlgG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,6680TA 6.8% 319.2 / 9.8 103intergenic (+111/‑61)flgF/flgGflagellar basal‑body rod protein FlgF/flagellar basal‑body rod protein FlgG
Reads supporting (aligned to +/- strand):  ref base T (50/46);  new base A (4/3);  total (54/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 104 total aligned reads displayed.

AACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCA  >  NC_000913/1134519‑1134809
                                                                                                                                                     |                                                                                                                                             
aaCGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGtt                                                                                                                                               <  1:949278/150‑1 (MQ=255)
  cGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGTGTCTATGGGGGCAATCGCAATTTATCAATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTa                                                                                                                                             >  1:469055/1‑150 (MQ=255)
    cAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                           <  1:978057/150‑1 (MQ=255)
      ggCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                         >  2:342259/1‑150 (MQ=255)
              caaccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                 >  1:529286/1‑150 (MQ=255)
              caaccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                 >  2:1576334/1‑150 (MQ=255)
                accaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                               <  1:1514342/150‑2 (MQ=255)
                 ccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACCAGTATAAGTTGCCCGAGGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcct                                                                                                                              <  2:1155783/150‑3 (MQ=255)
                    acTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaagaag                                                                                                                          <  2:1118253/150‑1 (MQ=255)
                     ctgctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaagaag                                                                                                                          >  1:481891/1‑150 (MQ=255)
                     ctgctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaagaag                                                                                                                          <  1:1249350/150‑1 (MQ=255)
                      tgctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcctacaag                                                                                                                         <  1:1527602/150‑8 (MQ=38)
                       gctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGAc                                                                                                                        >  1:1566258/1‑150 (MQ=255)
                         tgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACat                                                                                                                      >  2:94659/1‑150 (MQ=255)
                            cGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATg                                                                                                                   >  2:20947/1‑150 (MQ=255)
                            cGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATg                                                                                                                   >  1:283716/1‑150 (MQ=255)
                                   ttAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTg                                                                                                            >  1:223373/1‑150 (MQ=255)
                                   ttAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTg                                                                                                            >  2:1298955/1‑150 (MQ=255)
                                   ttAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTg                                                                                                            <  1:1312395/150‑1 (MQ=255)
                                   ttAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTg                                                                                                            >  2:238610/1‑150 (MQ=255)
                                    tAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGt                                                                                                           <  1:94659/150‑1 (MQ=255)
                                     aaTTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTc                                                                                                          >  2:148729/1‑150 (MQ=255)
                                      aTTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCa                                                                                                         <  2:869162/150‑1 (MQ=255)
                                       ttGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCAc                                                                                                        <  2:856642/150‑1 (MQ=255)
                                        tGAAAGGATACCTGAACATTATAATGTGCCCGATGCGCAAGTTTATCGTGGCTATGGGGGCAATCGAAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGCCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACt                                                                                                       <  2:653178/150‑1 (MQ=38)
                                        tGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACt                                                                                                       <  2:70351/150‑1 (MQ=255)
                                        tGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACt                                                                                                       >  1:363969/1‑150 (MQ=255)
                                         gAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTa                                                                                                      >  1:318375/1‑150 (MQ=255)
                                         gAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTa                                                                                                      >  1:502444/1‑150 (MQ=255)
                                         gAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTa                                                                                                      >  1:696189/1‑150 (MQ=255)
                                          aaaGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTaa                                                                                                     >  1:1623076/1‑150 (MQ=255)
                                             ggATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATcc                                                                                                  >  2:1037253/1‑150 (MQ=255)
                                                  cATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTAc                                                                                             <  2:324699/150‑1 (MQ=255)
                                                  cATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTAc                                                                                             >  2:273679/1‑150 (MQ=255)
                                                    tGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAg                                                                                           >  2:1586914/1‑150 (MQ=255)
                                                    tGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAg                                                                                           >  2:128741/1‑150 (MQ=255)
                                                       cAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGCCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGAc                                                                                        <  2:920379/150‑1 (MQ=255)
                                                         aGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACAt                                                                                      >  1:1268691/1‑150 (MQ=255)
                                                          gTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATTGGGGCAATCGCAATTTATCGATTTTGCGATCACTTGTAGGCCGCATAAGGCGTTTACGCCGCATCAGGCAAGAAGACATATGCACTTTGTCAGTAATCCACTACAGGAAAtt                                                                                     >  1:1056573/1‑150 (MQ=255)
                                                           tataAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACAttt                                                                                    <  1:1584496/150‑1 (MQ=255)
                                                             taAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTa                                                                                  <  2:283716/150‑1 (MQ=255)
                                                             taAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTa                                                                                  <  1:1425036/150‑1 (MQ=255)
                                                               aGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATg                                                                                <  1:342259/150‑1 (MQ=255)
                                                                gTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGa                                                                               <  1:1479001/150‑1 (MQ=255)
                                                                gTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGAAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGa                                                                               >  2:1272213/1‑150 (MQ=255)
                                                                  tGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATc                                                                             <  1:858785/150‑1 (MQ=255)
                                                                   gCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCa                                                                            <  1:614103/150‑1 (MQ=255)
                                                                   gCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCa                                                                            >  1:12835/1‑150 (MQ=255)
                                                                     ccGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGt                                                                          >  1:320343/1‑150 (MQ=255)
                                                                      cGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGtt                                                                         >  1:605265/1‑150 (MQ=255)
                                                                      cGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGtt                                                                         <  1:238610/150‑1 (MQ=255)
                                                                      cGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGtt                                                                         >  1:1303206/1‑150 (MQ=255)
                                                                          gcgcAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGGACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCAtt                                                                     <  2:12835/150‑1 (MQ=255)
                                                                           cgcAACTTTATCGGTAATATGGTGCCAAACGCTATTTATCGATTTTGCCAGGACATGTAAAACTAATAAGGCGTTTACGCCGCATCCGGCAAGAAGGCATATGCACTTTGTCACTAGTCCACTACAGGACAGTTTATGATCAGTTCAttt                                                                    >  2:177218/1‑149 (MQ=255)
                                                                             cAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATg                                                                  >  1:753595/1‑150 (MQ=255)
                                                                              aaGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATgg                                                                 >  1:834008/1‑150 (MQ=255)
                                                                              aaGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATgg                                                                 <  1:95359/150‑1 (MQ=255)
                                                                                 tttATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATACATATGCAGTTTGTCACTAATCCACTACAGTAGATTTTATGATCAGTTAATTATGGAta                                                              <  2:1625316/150‑2 (MQ=38)
                                                                                 tttATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATc                                                              <  1:1073697/150‑1 (MQ=255)
                                                                                 tttATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATc                                                              >  1:178948/1‑150 (MQ=255)
                                                                                  ttATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCg                                                             >  1:67760/1‑150 (MQ=255)
                                                                                         gTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAAc                                                      <  1:1165910/150‑1 (MQ=255)
                                                                                         gTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAAc                                                      >  2:77689/1‑150 (MQ=255)
                                                                                           cTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATc                                                                             >  1:168611/1‑125 (MQ=255)
                                                                                           cTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATc                                                                             <  2:168611/125‑1 (MQ=255)
                                                                                            tATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                   <  2:409363/150‑1 (MQ=255)
                                                                                            tATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                   <  2:1552679/150‑1 (MQ=255)
                                                                                             atcggggCAATCGAAATTTATCGATTTTGAGAGCACTTGTAGGCCGGATAAGGCGTTTACGCAGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGc                                                  >  2:1162805/4‑150 (MQ=255)
                                                                                             aTGGGGGCAATCGCAATTTATCAATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGc                                                  <  2:469055/150‑1 (MQ=255)
                                                                                               gggggCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCt                                                >  2:1406515/1‑150 (MQ=255)
                                                                                                   gCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGAc                                            <  2:14675/150‑1 (MQ=255)
                                                                                                       tCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGccc                                        <  1:1415488/150‑1 (MQ=255)
                                                                                                       tCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGccc                                        >  2:564935/1‑150 (MQ=255)
                                                                                                        cGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCa                                       >  2:934766/1‑150 (MQ=255)
                                                                                                         gCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAg                                      <  1:645381/150‑1 (MQ=255)
                                                                                                         gCAATTTATCGATTTTGCGAGAACTTGTAGGCCGGATAACGCGATTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCCCTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAg                                      >  2:222378/1‑150 (MQ=255)
                                                                                                               tATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAAcc                                <  1:1327624/150‑1 (MQ=255)
                                                                                                               tATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAAcc                                <  2:1268691/150‑1 (MQ=255)
                                                                                                                aTCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGc                                     >  2:1266367/1‑144 (MQ=255)
                                                                                                                aTCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGc                                     <  1:1266367/144‑1 (MQ=255)
                                                                                                                aTCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCa                               <  2:1566258/150‑1 (MQ=255)
                                                                                                                aTCGATTTTGCGAGCACTTGAAGGCCGGATAACGCGTTTACTCCGCATCCAGCAAGAAGAAATATGCACTTTGTCACTAATCCACTACAGGACAATTTATGATCAGTTCATTATGGATATCCAAAACGGGCCTTGACGCCCAGCAAACCa                               >  2:1453816/1‑150 (MQ=255)
                                                                                                                   gATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAata                            >  2:940993/1‑150 (MQ=255)
                                                                                                                   gATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAata                            >  2:1456129/1‑150 (MQ=255)
                                                                                                                     ttttGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATg                          >  1:316137/1‑150 (MQ=255)
                                                                                                                      tttGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATgg                         >  1:33260/1‑150 (MQ=255)
                                                                                                                       ttGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTCGACGCCCAGCAAACCAATATGGa                        <  1:77689/150‑1 (MQ=255)
                                                                                                                            aGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATAATATATATTCACTTTGTCACTAATCCAATACAGTAAATTTTATTATCAGTTCATTATGTAGAGGCAAAACGGTTCTTGACGCCCAGCAAACCAATATGGACGTCa                   <  2:696189/150‑1 (MQ=37)
                                                                                                                              cACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCAtt                 <  2:363969/150‑1 (MQ=255)
                                                                                                                                cTTGTAGGCCTGATAATGCGTTTACGCCGCATCCGGCAAGAAGACATATGCAGTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTgg               <  1:1369636/150‑2 (MQ=255)
                                                                                                                                cTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGc               >  2:1113099/1‑150 (MQ=255)
                                                                                                                                  tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCca             <  2:1236592/150‑1 (MQ=255)
                                                                                                                                   tcttgCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaa            >  2:131922/5‑150 (MQ=255)
                                                                                                                                    tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaac           <  2:605265/150‑1 (MQ=255)
                                                                                                                                    tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaac           <  1:934766/150‑1 (MQ=255)
                                                                                                                                    cttgCCGGATGCGGCGTAAAAGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaac           >  2:717847/4‑150 (MQ=255)
                                                                                                                                      tgCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaacaa         >  1:1393372/2‑150 (MQ=255)
                                                                                                                                       gCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCaacaac        >  2:1249365/1‑150 (MQ=255)
                                                                                                                                        ccGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAAcc       <  2:502444/150‑1 (MQ=255)
                                                                                                                                             tAAGGCGTTTACGCCGCATCCGGCAAGAATACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCa  <  2:67760/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
AACGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCA  >  NC_000913/1134519‑1134809

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: