Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,084,500 T→C 6.7% intergenic (‑311/+73) sdaC ← / ← amaC serine:H+ symport permease/D‑lysine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,084,5000TC 6.7% 552.1 / 4.5 186intergenic (‑311/+73)sdaC/amaCserine:H+ symport permease/D‑lysine aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (108/57);  new base C (2/11);  total (112/74)
Fisher's exact test for biased strand distribution p-value = 5.43e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 191 total aligned reads displayed.

GGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGCACACAC  >  NC_002947/4084351‑4084646
                                                                                                                                                     |                                                                                                                                                  
 gCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTc                                                                                                                                                   <  1:4870818/150‑1 (MQ=255)
   aTTTCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcgc                                                                                                                                                 <  1:2721924/150‑1 (MQ=255)
   aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcgc                                                                                                                                                 <  1:3555008/150‑1 (MQ=255)
     aactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTTCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGa                                                                                                                                               >  1:2009887‑M1/52‑150 (MQ=255)
     aactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGa                                                                                                                                               >  2:4131829‑M1/52‑150 (MQ=255)
      gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTTCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCGTTTAGCGAc                                                                                                                                              >  2:1773544/1‑150 (MQ=255)
      gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                                                                                              <  2:2519145/150‑1 (MQ=255)
      gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                                                                                              <  2:2519122/150‑1 (MQ=255)
      gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                                                                                              <  1:4888643/150‑1 (MQ=255)
      gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                                                                                              <  1:2385770/150‑1 (MQ=255)
          gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCaa                                                                                                                                          <  2:4835786/150‑1 (MQ=255)
           ccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAg                                                                                                                                         >  1:101790‑M1/46‑150 (MQ=255)
            caacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAgg                                                                                                                                        >  1:3012773‑M1/45‑150 (MQ=255)
             tCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGc                                                                                                                                       <  1:876621/150‑1 (MQ=255)
               gATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCg                                                                                                                                     <  1:5491412/150‑1 (MQ=255)
                 tCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGGCGCCCTAGGGCCGCAAACCGGCCCCAGCACTTTAAGCCGCAAAGCCGCCATCGTATGGGCGCGCCGCCCCCGGTTCGCGACGCAAGGCCGCt                                                                                                                                   >  2:4329703/1‑150 (MQ=255)
                  cGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGAAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCACGGCCGCTc                                                                                                                                  <  2:1593995/150‑1 (MQ=255)
                    gTGGAATGTTGCGGCCTGGGCGCCGGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGAAATGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATGGTTGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTccc                                                                                                                                <  2:3987308/150‑1 (MQ=37)
                         tgggcgcggcgTGGGCGCCTGGGGGGGGGCGTGGGGGCGCGGCGTTGCGGCGCGATTGGGCCGAAACGCGCCCCCAGCATTTTCCGCGGAAAAGCTGAAATCTTTGGGGCGCTCCGCACCATTTGCGCGAGGCAAGGCCGCTCCCACaaa                                                                                                                           <  2:1501290/145‑1 (MQ=255)
                         gcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGGCGCAAGGCCGCTCCCACaaa                                                                                                                           >  2:3083494‑M1/32‑150 (MQ=255)
                         aggtcgcggcgGGGGCGCGTGTGGGGCGGCGTGTGGGCGCGGCCTTGCGTGGCTATCGGGCCGCACCGCGGCCCCAGAATTTTCCGCGGCAAAGCTGCAAGGTTTGGGGCGCCCCGCACCCTTTTCGCGACGCACGGCCGCTCCCACaaa                                                                                                                           <  2:2370348/145‑1 (MQ=255)
                                  ccTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTTCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAATTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCGTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACcg                                                                                                                  >  2:4129142/1‑150 (MQ=255)
                                              tgGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCAAAGCGGCCCCCGCAAATTGAGCGGCAAAGCAGCAATCGTAGGGGCGCTCCGCACCCTTTTAGCGACCGAAGGCCGCTCCCACAAAGTCCAACCCCCCCCCCCagga                                                                                                      >  2:3089840/1‑150 (MQ=25)
                                                  gcggtCTGTGGGTGCGGCCTGGCGTCGCGATAGGGCCGCATAGCGGCCCCCGCAATTTCAGAGTATTAGCTGAAATCGTCGGGGCGCTCCGACCCCTTTTCGCGATGCTAGGCCGCTCCAAAAAAGTCCAACCGCGCCCCACAGGAgggg                                                                                                  <  2:5614332/150‑1 (MQ=25)
                                                  gcggcCGGTGGGAGCGGCCTGGCGGCGCGAGCGGGCCGCAACGCGGCCCCAGCACTTCCCGCGGCAACGCTGCAATCGTCGGGGCGCTCCGCACCCTGTGCGCGCCGCAAGGCCGCTCCCACAAAGTCCACCCGCGCCCCACAGGAgggg                                                                                                  <  2:3681727/150‑1 (MQ=25)
                                                      ccTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCACCAAATTCAGCGGCAAAGCGGCAATCGTAGGGGCGCGACGCACCCCTTTCGCGACGCAAGGCCGCACCCACCAAGTCCAACCGCGCCCCACCGGAGGGGCCGc                                                                                              >  2:870584/1‑150 (MQ=25)
                                                      ccTGTGGGAGCGGCCTAGCGTCCCGATAGGGCCCCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGGGCTCCGCACCCTTTTGGCGACGCAAGGCCGCTCCCACAAAGTCCAAACGCGCCCCACAGGAGGGGCCGc                                                                                              >  2:5456203/1‑150 (MQ=37)
                                                       cTGTGGGAGCGGCCATGCGTCGCGATAGGGCCGCAAATCGGCCCCACCCAATTAAGCGGCAAATCAGCAAAACTATGGGCCCTCCGCACCCATTTAGCCAACCAAAGCCCCTCCCAAAAAATACAAACCCGCCCCACACGAGGGGCCCCa                                                                                             >  2:4293941/1‑150 (MQ=12)
                                                               GCGGCCTTGCGTCGCCATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAATGCTGCAATCGACGGGGCGCACCGCACCCTTTCCGCGACGCATGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGccc                                                                                     >  2:5299493‑M1/1‑150 (MQ=255)
                                                                 GGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGt                                                                                   <  2:4517473‑M1/150‑1 (MQ=255)
                                                                     TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCATCCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTtggg                                                                               <  2:1174201‑M1/150‑2 (MQ=255)
                                                                     TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc                                                                               <  2:1888535‑M1/150‑1 (MQ=255)
                                                                     TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTGTCGCGACGCAAGGCCGCTCCCACACAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc                                                                               <  2:4721949‑M1/150‑1 (MQ=255)
                                                                      TGCGTCGCTTTAGGGCCGCAAAGCGGCCCATGCATTGTCAGCGGCAAAGCTGAAATCGTCGGGGCGCTCCGCACCCTTTGCGCCACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCg                                                                              >  2:634105‑M1/1‑150 (MQ=255)
                                                                        CGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGAGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAg                                                                            <  2:2804740‑M1/150‑1 (MQ=255)
                                                                          TCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGtt                                                                          <  2:732415‑M1/150‑1 (MQ=255)
                                                                             CGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACt                                                                       <  1:5355628‑M1/150‑1 (MQ=255)
                                                                               ATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTc                                                                     <  2:1321993‑M1/150‑1 (MQ=255)
                                                                                tggggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCt                                                                    >  2:4971687/3‑150 (MQ=255)
                                                                                tggggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCt                                                                    >  2:4971678/3‑150 (MQ=255)
                                                                                tggggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCt                                                                    >  1:3296735/3‑150 (MQ=255)
                                                                                tggggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCt                                                                    >  1:3736357/3‑150 (MQ=255)
                                                                                tggggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCt                                                                    >  1:3861815/3‑150 (MQ=255)
                                                                                   ggCCGCAACGCGGCCCCAGCAATTTCCGCGGAAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa                                                                 <  2:3927344/150‑1 (MQ=255)
                                                                                   ggCCGCAACGCGGCCCCAGCAATTTCAGCGGCAACGCTGCAATCGTTGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa                                                                 <  2:314143/150‑1 (MQ=255)
                                                                                   ggCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCaa                                                                                                                     <  2:2560260/98‑1 (MQ=255)
                                                                                   ggCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCaa                                                                                                                     >  1:2560260/1‑98 (MQ=255)
                                                                                    gCCGCAACGCGGCCCCCGCAATTTCCGCGGCAAAGCTGAAAGGGTCGGGGCGCGCCGCACCCTTTGCGCGACGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAAc                                                                <  2:4704147/150‑1 (MQ=255)
                                                                                    gCCGCAAAGCGGCCCCAGCGATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAAc                                                                >  1:1550156/1‑150 (MQ=255)
                                                                                       gCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGc                                                             >  1:192419/1‑150 (MQ=255)
                                                                                       gCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGc                                                             >  2:3120289/1‑150 (MQ=255)
                                                                                        ctacgCGGCCCCTGCTTTTTCAGCGGCTTAGCTGCAAACGTTGGGGGGCTCGGCACGCTTTTCGCGACGAAAGGCCGCTGCCTCAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCTGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCa                                                            <  2:1817770/146‑1 (MQ=255)
                                                                                        cAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCa                                                            >  2:2329087/1‑150 (MQ=255)
                                                                                         aaaGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAg                                                           >  1:3107548/1‑150 (MQ=255)
                                                                                         aaaGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAg                                                           >  1:993770/1‑150 (MQ=255)
                                                                                         aaaGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAg                                                           >  1:2027762/1‑150 (MQ=255)
                                                                                           aGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCg                                                         >  1:4157596/1‑150 (MQ=255)
                                                                                            gCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGa                                                        >  2:3844855/1‑150 (MQ=255)
                                                                                               gCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACg                                                     >  1:233534/1‑150 (MQ=255)
                                                                                                ccccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAg                                                                                            >  2:1797446/1‑110 (MQ=255)
                                                                                                ccccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAg                                                                                            <  1:1797446/110‑1 (MQ=255)
                                                                                                  ccAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGTCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCtt                                                  >  1:1646994/1‑150 (MQ=255)
                                                                                                   cAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGa                                                                             >  1:2171513/1‑122 (MQ=255)
                                                                                                   cAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTc                                                 >  1:1874867/1‑150 (MQ=255)
                                                                                                   cAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTc                                                 >  1:2880421/1‑150 (MQ=255)
                                                                                                   cAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTc                                                 >  1:3139190/1‑150 (MQ=255)
                                                                                                   cAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTc                                                 >  1:3767027/1‑150 (MQ=255)
                                                                                                     gCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCgg                                               >  1:4447031/1‑150 (MQ=255)
                                                                                                         tttCCGCGGCAAAGCTGAAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAAc                                           <  2:2200817/150‑1 (MQ=255)
                                                                                                         tttCAGCGGCACAGCTGCACTCGTCGGGGCGCTCCGCACCCCTTTCGCGCCGCAAGCCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAAc                                           <  1:3743032/150‑1 (MQ=255)
                                                                                                         tttCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAAc                                           >  2:4305903/1‑150 (MQ=255)
                                                                                                         tttCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAAc                                           >  2:2738564/1‑150 (MQ=255)
                                                                                                           tCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAAccc                                         >  1:4451292/1‑150 (MQ=255)
                                                                                                            cAGCGGCAATGCTTCAAGCGTCGGGGCGCTCCCCACCCTTTCCGCCCCCCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCg                                        <  1:1773544/150‑1 (MQ=255)
                                                                                                              gCGTCAAAGCTGCACTCGTTGGGGCGCTCCCCCCCCCTCCCGCCCCCCACCCCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGtt                                      <  1:4999403/150‑1 (MQ=255)
                                                                                                               gggcaccgctgctCTCGTCGGGGCGCCCCCCACCCTCCCCGCCACGCCCGCCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTg                                     <  1:4427164/143‑1 (MQ=255)
                                                                                                                 gCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCa                                   >  2:2134310/1‑150 (MQ=255)
                                                                                                                   aaaGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAgg                                 >  2:4089899/1‑150 (MQ=255)
                                                                                                                   aaaGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAgg                                 >  2:910037/1‑150 (MQ=255)
                                                                                                                        tGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTgtc                            <  2:203950/150‑3 (MQ=255)
                                                                                                                        tGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGc                            >  2:5597679/1‑150 (MQ=255)
                                                                                                                        tGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGc                            >  1:121366/1‑150 (MQ=255)
                                                                                                                           aaTCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCgg                         >  2:3390442/1‑150 (MQ=255)
                                                                                                                           aaTCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCgg                         >  1:2801597/1‑150 (MQ=255)
                                                                                                                               gTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGgcgc                     <  2:4639875/150‑1 (MQ=255)
                                                                                                                               gTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGgcgc                     >  1:2501270/1‑150 (MQ=255)
                                                                                                                               gTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGgcgc                     >  1:4206382/1‑150 (MQ=255)
                                                                                                                               gTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGgcgc                     >  1:2961346/1‑150 (MQ=255)
                                                                                                                                tcggggCGCTCCGCCCCCTTTTCGCGACGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCa                    <  2:2925353/148‑1 (MQ=255)
                                                                                                                                  ggggCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAgg                  >  1:1935068/1‑150 (MQ=255)
                                                                                                                                    ggCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGcc                >  1:4732920/1‑150 (MQ=255)
                                                                                                                                          ccGCACCCGTTTCGCCACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCGCGGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATg         <  2:5475857/150‑1 (MQ=255)
                                                                                                                                            gCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGc        >  2:974793/1‑150 (MQ=255)
                                                                                                                                            gCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGc        >  1:4364733/1‑150 (MQ=255)
                                                                                                                                             cACCCTTGTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGca       <  2:4218419/150‑1 (MQ=255)
                                                                                                                                               cccTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGcaca     >  2:4396904/1‑150 (MQ=255)
                                                                                                                                                ccTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGcacac    >  1:4499670/1‑150 (MQ=255)
                                                                                                                                                ccTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGcacac    >  2:3421452/1‑150 (MQ=255)
                                                                                                                                                  ttttCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGCacacac  >  1:5259416/1‑150 (MQ=255)
                                                                                                                                                  ttttCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGCacacac  >  2:4227372/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAGCGAACGCTTCGGCAACCCGTTGCAGGTTGGCCGGGCGCAGGCCGGACATGCACACAC  >  NC_002947/4084351‑4084646

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