Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,084,441 A→T 5.3% intergenic (‑252/+132) 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,4410AT 5.3% 735.3 / 22.3 226intergenic (‑252/+132)sdaC/amaCserine:H+ symport permease/D‑lysine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (78/134);  new base T (7/5);  total (87/139)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 229 total aligned reads displayed.

GGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAG  >  NC_002947/4084293‑4084589
                                                                                                                                                    |                                                                                                                                                    
  ttGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                                                   >  2:1493948/1‑150 (MQ=255)
     cACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGc                                                                                                                                                <  2:860294/150‑1 (MQ=255)
     cACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGc                                                                                                                                                <  2:5625410/150‑1 (MQ=255)
      aCTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTTTGCGGcc                                                                                                                                               <  2:423167/150‑1 (MQ=255)
      aCTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcc                                                                                                                                               <  2:1888554/150‑1 (MQ=255)
      aCTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcc                                                                                                                                               <  1:1906870/150‑1 (MQ=255)
        ttGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGcc                                                                                                                                               >  2:687550/1‑148 (MQ=255)
        ttGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcccc                                                                                                                                             <  1:556571/150‑1 (MQ=255)
           cGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTTTGCGGccctatc                                                                                                                                          <  1:3538288/150‑5 (MQ=255)
           cGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccctatc                                                                                                                                          <  2:4139045/150‑5 (MQ=255)
            gAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGCTATGGCCGCAAAGCGGCCCCAGCa                                                                                                                                         >  2:3512323/1‑150 (MQ=255)
            gAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCGCa                                                                                                                                         >  2:6007414/1‑150 (MQ=255)
            gAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCGCAAAGCCGCCCCAGCa                                                                                                                                         >  1:7268749/1‑150 (MQ=255)
             cccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCaa                                                                                                                                        <  2:693759‑M1/49‑1 (MQ=255)
              aGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGGAGCGGCCTGGGCGCCGGTGGGGCGGCCTGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCAGCAAt                                                                                                                                       <  1:2537119/150‑1 (MQ=255)
                cccTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAttt                                                                                                                                     <  2:1541395/150‑1 (MQ=255)
                    tttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGc                                                                                                                                 <  1:5536216/150‑1 (MQ=255)
                     ttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATTGGGCCGCAAAGCGGCCCCAGCAATTTCAGCg                                                                                                                                <  2:829950/150‑1 (MQ=255)
                      tCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCgg                                                                                                                               >  2:3487497/1‑150 (MQ=255)
                      tCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCgg                                                                                                                               <  2:5930678/150‑1 (MQ=255)
                          aggagCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCaaa                                                                                                                           <  1:2189835/150‑1 (MQ=255)
                          aggagCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCaaa                                                                                                                           <  1:5909112/150‑1 (MQ=255)
                          aggagCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGCCCCAGCAAAATAAGCGGCaaa                                                                                                                           >  1:6711602/1‑150 (MQ=39)
                             agCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCt                                                                                                                        <  2:1817028/150‑1 (MQ=255)
                              gCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGTGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTg                                                                                                                       <  1:6736324/150‑1 (MQ=255)
                              gCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTg                                                                                                                       <  2:4515835/150‑1 (MQ=255)
                              gCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGACAAAGCTg                                                                                                                       <  1:1229536/150‑1 (MQ=255)
                              gCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGAGAGGGCCGCACAGCGGCCCCCCCAAGATCAGCGGCGAAGcag                                                                                                                       >  2:5026390/1‑148 (MQ=255)
                               cGGTTGCACTTGACTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGc                                                                                                                      <  2:2984749/150‑1 (MQ=255)
                                 gTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCaa                                                                                                                    >  2:371829/1‑150 (MQ=255)
                                 gTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCaa                                                                                                                    <  2:5573033/150‑1 (MQ=255)
                                   tGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATc                                                                                                                  <  1:6134300/150‑1 (MQ=255)
                                   tGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATc                                                                                                                  <  2:3111221/150‑1 (MQ=255)
                                     cgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGt                                                                                                                <  1:3430972‑M1/73‑1 (MQ=255)
                                     cACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGt                                                                                                                <  1:2844644/150‑1 (MQ=255)
                                      aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTg                                                                                                                       <  1:6894482/142‑1 (MQ=255)
                                      aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTg                                                                                                                       >  2:6894482/1‑142 (MQ=255)
                                       cTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGCGTCGCGGTTGGGCCGCAAAGCGGCCCCACCAAATTAAGCCGCCAAGCAGCAAACCTAg                                                                                                              >  1:4760764/1‑150 (MQ=255)
                                       cTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAg                                                                                                              <  2:5208823/150‑1 (MQ=255)
                                       aagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAg                                                                                                              >  1:4821189‑M1/76‑150 (MQ=255)
                                        ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTTCGTCGCGATAGGGGCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAAACGTAgg                                                                                                             >  2:3594201/1‑150 (MQ=255)
                                        ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                                                   <  2:3152845/112‑1 (MQ=255)
                                        ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgg                                                                                                             <  2:2477184/150‑1 (MQ=255)
                                          gTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgggg                                                                                                           <  1:83443/150‑1 (MQ=255)
                                              ttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGTGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCt                                                                                                       <  2:6964159/150‑1 (MQ=255)
                                                aCAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTcc                                                                                                     <  1:3074202/150‑1 (MQ=255)
                                                aCAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCAAAGCGGCCCCCGCAAGTTCAGCGGCAAAGCGGCAATCGGAGGGGCGCTcc                                                                                                     >  2:1686474/1‑150 (MQ=39)
                                                   aGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCa                                                                                                  <  2:4598063/150‑1 (MQ=255)
                                                   aGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCa                                                                                                  <  2:912594/150‑1 (MQ=255)
                                                     cGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCAcc                                                                                                <  1:2996629/150‑1 (MQ=255)
                                                       gAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAATTCAGCGGCAAAGCTGCAATAGTAGGGGCGCTCCGCAccca                                                                                              >  2:6528250/1‑149 (MQ=255)
                                                        aaGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCtt                                                                                             <  1:7345751/150‑1 (MQ=255)
                                                         agaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCttt                                                                                            >  1:694042‑M1/58‑150 (MQ=255)
                                                          ggCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCtttt                                                                                           <  2:257961/150‑1 (MQ=255)
                                                          ggCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCtttt                                                                                           <  1:6516305/150‑1 (MQ=255)
                                                          ggCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGGCGCCATAGGGCCGCAAAGCGGCCCCAGCAAATTAAGCCGCAAAGCAGCAATAGTAGGGGCGCTCCCCACCCattt                                                                                           >  2:4652991/1‑150 (MQ=38)
                                                            cATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcg                                                                                         <  1:846270/150‑1 (MQ=255)
                                                               tGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGa                                                                                      <  2:5351861/150‑1 (MQ=255)
                                                                gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                                     <  1:371829/150‑1 (MQ=255)
                                                                  tGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCATTCCTCGCCATATGGCCGCAAACCCGCCCCACCCATTTCAGCCGCAAAGCTGCAATACTAGGGGCCCACCGCACCCATTTCCCGAAGc                                                                                   >  2:2378769/1‑150 (MQ=255)
                                                                    gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCaa                                                                                 <  1:1794246/150‑1 (MQ=255)
                                                                    gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCaa                                                                                 <  1:4956313/150‑1 (MQ=255)
                                                                     ccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAg                                                                                >  2:6550665‑M1/46‑150 (MQ=255)
                                                                      gtCGATCGGTGGAATGTATCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAgg                                                                               <  1:1079196/150‑1 (MQ=255)
                                                                      gtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAgg                                                                               <  2:3204128/150‑1 (MQ=255)
                                                                      caacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCACCGCACCCTTTTCGCGACGCAAgg                                                                               >  2:3430972‑M1/45‑150 (MQ=255)
                                                                          aTCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCCCAAATTCAGCGGCAAAGCGGCCATCGTAGGGGCGCTCCCCCCCCGGTTCGCGACGCACGGCCGc                                                                           >  2:3896314/1‑150 (MQ=37)
                                                                              gTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCCCAAAGCCGCCCCAGCAAATTAAGCGGCCAACCAGCAAACGTAGGGGCGCGCCCCACCCGTTTCGCGACCCAAGGCCGCTccc                                                                       >  2:2125099/1‑150 (MQ=255)
                                                                               tGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCATCGGCAAAGCTGCCATCGTAGGGGCGCGCCGCACCCATTTAGCGACGCACGGCCGCTCCca                                                                      >  2:7219936/1‑150 (MQ=37)
                                                                                 gAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTTCGTCGCCGTAGGGCCGCAAAACGGCCCCAGCAAAATCAGCGGCCAAGCGGCAAACGGTGGGGCGCGCCGCCCCCGTTTCGCGGCGCAAGGCCGCTCCcaca                                                                    >  2:6403659/1‑150 (MQ=255)
                                                                                   aTGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTCGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACaaa                                                                  <  2:5381797/150‑1 (MQ=255)
                                                                                             cTGGGCGCCTGTGGGGCGGCGTGTGGGAGGGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCAGCAATTTCAGCGGAAACGCTGCATTCTTAGGGGCGCTCCGCACCCTTTTCGCGCCGCACGGCCGCTCCCACAAAGTCCAACcgc                                                        <  2:3836609/150‑1 (MQ=255)
                                                                                                    cggggggggCGGCCTGGGGGAGCGGCCTGGCGCCGCGATCGGGCCGCAAAGCGCCCCCAGAAATGCCCGCGGCAAAGCTGCAATCGGAGGGGCGCTCCGCACCCTTTCCGCGCCGCATGCCCGCTCCCACAAAGTCCAACCGCGCCCcac                                                 <  2:6378488/145‑1 (MQ=255)
                                                                                                      tgtgGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCt                                                                                                       >  1:81995/1‑94 (MQ=255)
                                                                                                           ggcggcCTGTGGGAGCGGCCATTCGTCGCCATATGGCCGCAAACCGGCCCCACCAAATTAAGCGGCAAAGCTGCAATAGGAGGGGCGCGCCCCACCCATTTCGCGACGCCAGGCCCCGCCCACAAAGTCCAACCCCCCCCCCCACGAggg                                          >  2:5091672/1‑150 (MQ=255)
                                                                                                                    TGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGtt                                 >  1:4541554‑M1/1‑150 (MQ=255)
                                                                                                                    TGGGAGCGGCCTTGCGTCGCGAGCGGGCCGCACAGCGGCCCCCGCAATTGCAGCGGCAACGCTGCAATCGTAGGGGCGCTCCGCCCCCATTGCGCGCCGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGtt                                 >  2:4799619‑M1/1‑150 (MQ=255)
                                                                                                                             CCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTg                        <  1:3299734‑M1/150‑1 (MQ=255)
                                                                                                                              CTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCACTCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTgg                       <  2:6343932‑M1/150‑1 (MQ=255)
                                                                                                                              CTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTgg                       <  2:6343497‑M1/150‑1 (MQ=255)
                                                                                                                               TTGCGTCGCGATCGGGCCGAAACGCGCCCCCAGAATTTCCCGCGCAAAAGCTGAAATCGTTGGGGCGCTCCGCACCCTTTTCGCGACGAACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc                      >  2:3332972‑M1/1‑150 (MQ=255)
                                                                                                                                 GCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGa                    <  1:2346064‑M1/150‑1 (MQ=255)
                                                                                                                                    TCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGtt                 <  1:4226727‑M1/150‑1 (MQ=255)
                                                                                                                                      gctatcgggCCGCAAAGCGGCCCCCGCAATTTCCGCGGCAAAGGTGCAATCGTCGGGGCGCTCCGCACCCTTTTCGCGCCGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTAc               <  2:810011/144‑1 (MQ=255)
                                                                                                                                        gctggggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACtt             >  1:3760582/5‑150 (MQ=255)
                                                                                                                                         ATAGGGCCGCAACGCGGCCCCAGCAATTTCAGCGGCAACGCGGCAAGCGTCGGGGCGCTCCGCACCCTTTTCGCGCCGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTc            >  2:6711602‑M1/1‑150 (MQ=255)
                                                                                                                                         ATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTc            <  2:6068200‑M1/150‑1 (MQ=255)
                                                                                                                                         ATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTc            <  1:1001392‑M1/150‑1 (MQ=255)
                                                                                                                                          tggggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCt           >  2:1621302/3‑150 (MQ=255)
                                                                                                                                          tggggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCt           >  2:750910/3‑150 (MQ=255)
                                                                                                                                          tggggCCGCTTTGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCt           >  2:2854325/3‑150 (MQ=255)
                                                                                                                                          tAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCt           <  2:4226727/150‑1 (MQ=255)
                                                                                                                                             ggCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGaa        >  1:6832364/1‑150 (MQ=255)
                                                                                                                                              gCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAAc       >  1:5137154/1‑150 (MQ=255)
                                                                                                                                               ccGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACg      >  1:1993755/1‑150 (MQ=255)
                                                                                                                                               ccGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACg      >  2:2942300/1‑150 (MQ=255)
                                                                                                                                                cGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACaaa                                                                  >  1:1404285/1‑89 (MQ=255)
                                                                                                                                                cGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACgg     <  2:1001392/150‑1 (MQ=255)
                                                                                                                                                cGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACgg     >  1:4350841/1‑150 (MQ=255)
                                                                                                                                                cGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACgg     >  2:1079891/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTCTGAACGGCAG  >  NC_002947/4084293‑4084589

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: