Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,531 C→A 28.6% intergenic (+48/+66) mdfA → / ← ampG multidrug MFS transporter/muropeptide permease AmpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,544,5310CA28.6% 190.7 / 91.5 125intergenic (+48/+66)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base C (28/61);  new base A (25/11);  total (53/72)
Fisher's exact test for biased strand distribution p-value = 1.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 132 total aligned reads displayed.

TGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGA  >  NC_002947/1544383‑1544679
                                                                                                                                                      |                                                                                                                                                     
tgcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                      <  2:3297650/150‑1 (MQ=255)
    tttgGCCGTGGGGGGGGGGGGGGGGGGGGGGTTGGCGCGGGGGCCGGTGTGGTTTGGGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgc                                                                                                                                                  <  1:2973227/150‑1 (MQ=38)
         ccgtGGGGGGGGGGGTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgcgc                                                                                                                                             <  1:3519633/147‑1 (MQ=255)
           cTGGCGGGGGGGGGGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgcgccg                                                                                                                                           <  1:4986185/150‑1 (MQ=255)
              gcggcgGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGgccgcgccgcgg                                                                                                                                        >  2:3476210/1‑148 (MQ=255)
               ggggggggggggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCttt                                                                                                                                       <  1:4694296/146‑1 (MQ=255)
                 ggggggggggggggggggTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCACCGCTTTGt                                                                                                                                     <  1:4233120/148‑1 (MQ=255)
                 ggggggggggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                     <  1:4876553/148‑1 (MQ=255)
                 ggggggggggggggggCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTgtgt                                                                                                                                     <  1:5105081/148‑1 (MQ=25)
                 gggggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                     <  1:243421/148‑1 (MQ=255)
                 gcggGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCGCTTCGCGGGGAAACCCGCTCCTACACAGGgccgcgccgcggtgt                                                                                                                                     >  2:1095522/1‑145 (MQ=255)
                  gggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                    <  2:606158/149‑1 (MQ=21)
                  gggggggggCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                    <  1:1746850/149‑1 (MQ=255)
                  ggggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                    <  1:2054935/149‑1 (MQ=255)
                  ggggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                    <  1:3465165/149‑1 (MQ=255)
                  ggggggTGGCGGGGGCGTTGGCGCTGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                    <  2:2068339/149‑1 (MQ=255)
                  cggGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                    <  2:3499199/150‑1 (MQ=37)
                     ggggggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                 <  1:4161717/150‑1 (MQ=255)
                     gggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                 <  1:754123/150‑1 (MQ=255)
                     gggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                 <  1:1434842/150‑1 (MQ=255)
                      ggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAg                                                                                                                                >  2:3045870/1‑150 (MQ=255)
                       ggggggggggCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGc                                                                                                                               <  1:3322779/145‑1 (MQ=255)
                       ggggggggggCGTGGGCGCGGGGGCCGCTGTGTTTTGGGGGGGGTGGGCGGCCTGGGGAAGGGTCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGc                                                                                                                               <  1:3339026/145‑1 (MQ=38)
                       gTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGc                                                                                                                               <  2:1068705/150‑1 (MQ=37)
                       gTGGCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGc                                                                                                                               <  1:3322753/150‑1 (MQ=255)
                         gggggggggGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                             <  1:1835738/147‑1 (MQ=255)
                         ggggggggCGTGGGGGCTGTGGCCGCTGTGGTTTGGGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                             <  1:4618011/147‑1 (MQ=38)
                         ggCGGGGGCGTTGGCGCTGTGGCCGCTGGGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                             <  2:1433258/150‑1 (MQ=255)
                          gggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg                                                                                                                            <  2:560919/148‑1 (MQ=255)
                          gCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg                                                                                                                            <  2:328723/150‑1 (MQ=255)
                             ggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGc                                                                                                                               <  1:965506/144‑1 (MQ=255)
                             ggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGc                                                                                                                               >  2:965506/1‑144 (MQ=255)
                                  gTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCACCTCGCCGGGTAAACCCCACCAACAAACGGGCCCCCCGCCTTGTTGGTGCGGGTGTTAccc                                                                                                                    >  1:2120903/1‑150 (MQ=255)
                                   ttGGCGCTGTGGCCGCGGTGGTTGGTGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCcgc                                                                                                                  <  1:2638728/150‑1 (MQ=21)
                                    gggCGCTGGGGCCGGGGTGGTTGGGGGGGGGGGGGCGCCGTGGGGACGGGGCGGTAGCTCGATACTCGGTATGGCCGGGGCCGGCCTCTTCGGGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCcgcg                                                                                                                 <  2:5157925/149‑1 (MQ=9)
                                      gcgcTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCCCCACAACAAACGGCCCCCCCCCGGTGTAGGGGCGGGGTTAccccccca                                                                                                               >  1:5412449/1‑145 (MQ=255)
                                        ggtgtgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAaga                                                                                                             <  2:2177054/148‑1 (MQ=21)
                                        gcTGTGGCCGCTGGGGTTGGTGGTGGGGGGGCGGCCGGGGGACGGGGCTGTCGCTCGATAATCGTTATGGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACTGGGCCGCGCCGCTTTGTCGTAGGGGTTTTACCCGCGAaga                                                                                                             <  2:1134062/150‑1 (MQ=37)
                                                    gtgGTTGGTGGGGGGTGGGGGGCGTGGGGATGGGGCTTTAGCTCTATAATCGTTATGGCGGGGGCCGGCCTCTTGGGGGTTAAACCCGCTCCTACACCGGGCCGCGCCGATTTGTAGGAGGGTTTTTACCCGCGAAGAGGCCAGCACAgg                                                                                                 <  2:4268387/150‑1 (MQ=255)
                                                           gtggtgGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACa                                                                                          <  2:4081260/150‑1 (MQ=255)
                                                             ggggggggggCGGCGGGGGGACGGGGCTGTTGCGCGATATGGGGTTTTGCGGGGGCCGGCCTCTACGCGGGTAAACCGGCTCCTCCACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACaaa                                                                                        <  2:3939726/147‑1 (MQ=21)
                                                                   ggggCGGCCTGGGGCAGGGGCTGTAGCTCGATCATCGGTATGGCCGGGGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCa                                                                                  <  2:1696734/149‑1 (MQ=255)
                                                                      gcggcCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCCCACCTACACAGGGCCGCGCCGCGTTGTAGGGGCCGGGTTACCCCCGGAGAGGCCACCACAGCCAATACAAAGCGTCAggg                                                                               >  2:3498855/1‑150 (MQ=37)
                                                                           cTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCCCTCCTAAACACGGCCCCGCCGCGGTGTAGGGGCGGGTTTAACCCCCAAGAGGCCCGCACCGCCAAAAAAAAACGTAAGGGgcgcg                                                                          >  2:3912123/1‑149 (MQ=255)
                                                                            tGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                 <  1:1603671/95‑1 (MQ=255)
                                                                                       cggTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGt                                                              <  2:830093/148‑1 (MQ=255)
                                                                                             cTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCCCACCGACACAGGGCCCCGCCGCGTTGTAGGATCGGGTTTACCCGCGCAGAGGCCCGCACAGGCAAAAAAAAGCATCAGGGCCGCGACGCCACCGACTGggcgt                                                        >  2:1736280/1‑150 (MQ=25)
                                                                                             cTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGGAAAACCGCGACTACACAGGGCCCCCCCGCGTTGTTGGGGCGGGGTTAACCGCGAAGAGGCCAGCACAGGCCATACAAAGCATAAGGGCCGCGCAGCCACCGAGTGgtcgg                                                        >  2:2255329/1‑149 (MQ=255)
                                                                                             cTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGGAAAACCCCTCCAAAAAACGGCCCCGCCCCGTTGTTGGAGCGGGTTTTACCGCGAAGAGGCCACCACACGCAATACAAAGCAAAAGGGCCGCGCAGCCACCCAGTGgtcgt                                                        >  2:1736201/1‑150 (MQ=255)
                                                                                              tctaTAATCGTTTTTGCGGGGGCGGGCCTCTTCGCGGGAAACCCCGCCCCTACACCGGGCGGCGCGGCTTTGTCGGTGCGGGTTACCCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGgtcgtc                                                       <  1:784245/147‑1 (MQ=255)
                                                                                              tCGATAATCGGTAGGGCCGGGGCCGGCCTCTGCGCGGGTAACCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGgtcgtc                                                       <  1:4057313/150‑1 (MQ=255)
                                                                                                    aTCGGTATTGCCGGGGCCGGCCTCTTCGCGGGGAAAACCCCACCCACACAGGGCCGCCCCCCGTTGTAGGAGCGGGTTTAACCCCGCAGAGGCCAGCACAGGCAAAACAAAGCTTCAGGGCCGCGCACCCCCCCACTGGTCGTCCGcccc                                                 >  2:3595357/1‑148 (MQ=255)
                                                                                                         ttgggcgggggCCGGCCTCTGCGCGGCAACCCCCGCCCCCCCCCCGGCCCGCCCCGCTTTGGCGGCGCGTTTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCg                                            <  1:4631120/143‑1 (MQ=255)
                                                                                                          aTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg                                                                                                                            <  2:4997244/70‑1 (MQ=14)
                                                                                                          aTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg                                                                                                                            >  1:4997244/1‑70 (MQ=14)
                                                                                                           tggCCGGGGCCGGCCGCTGCGGGGGTAACCCCGCGCCTCCACCGGGCCGCGCCGGTTTGTCGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTc                                          <  2:2748267/148‑1 (MQ=38)
                                                                                                             gcgggggCGGGCCTCTGGGGGGCAAACCCCGCCCATACCCTGGGCCGCGCGGGTTTGGAGGGGGGGGTTCACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTccc                                        <  1:597282/147‑1 (MQ=255)
                                                                                                                 gggCCGGCCTCTTCGCGGGTAAACCGGCTCCTACACCGGGCCGCGCCGCTTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgc                                    <  2:2924367/150‑1 (MQ=255)
                                                                                                                 gggCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTGGTAGGAGCGGTTTTCCCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgc                                    <  2:194111/150‑1 (MQ=255)
                                                                                                                      ggccTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGc                               >  2:422591/1‑150 (MQ=255)
                                                                                                                      ggccTCCTGGCGGTAAACCCGGCACATACACTGGCCCGCGCGGTTTTGTAGGCGGGGTTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGc                               <  1:3728544/150‑1 (MQ=255)
                                                                                                                         cactTCGCGGGAAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc                            <  1:2195603/148‑1 (MQ=255)
                                                                                                                               gcgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATg                      >  1:1084470/1‑150 (MQ=255)
                                                                                                                                  gggTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCa                                                                                                        >  1:3503430/1‑66 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCAc                  >  2:1778389/1‑150 (MQ=255)
                                                                                                                                    ctAAACCCGCCCCTCCACCGGGCCGCGCCGCTTTGTTGGCGCGGGTTTCCCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCAc                  <  1:3718299/149‑1 (MQ=255)
                                                                                                                                      aaaCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg                >  2:4176649/1‑150 (MQ=255)
                                                                                                                                       aaCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGc               >  1:625551/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg            >  1:185729/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg             <  2:3674118/149‑1 (MQ=255)
                                                                                                                                           cGCTCCACCACCGGGCCGCGCCGGTTTGATGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg           <  1:1095522/150‑1 (MQ=39)
                                                                                                                                             gCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc          >  1:3398083/1‑150 (MQ=255)
                                                                                                                                             gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc          >  1:1156157/1‑150 (MQ=37)
                                                                                                                                               tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga        >  2:4043353/1‑150 (MQ=255)
                                                                                                                                               tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga        >  1:2260766/1‑150 (MQ=37)
                                                                                                                                               tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga        >  2:2427226/1‑150 (MQ=37)
                                                                                                                                               tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga        >  1:1756532/1‑150 (MQ=37)
                                                                                                                                               tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga        >  1:3886408/1‑150 (MQ=37)
                                                                                                                                                ccTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat       >  1:3643757/1‑150 (MQ=255)
                                                                                                                                                ccTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat       >  1:4289398/1‑150 (MQ=37)
                                                                                                                                                 cTACACTGGGCCGCGCCGCTTTGTCGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc      <  1:3476210/150‑1 (MQ=255)
                                                                                                                                                 cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc      >  1:3666533/1‑150 (MQ=25)
                                                                                                                                                 cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGgtcgtc                                                       >  1:5241839/1‑101 (MQ=18)
                                                                                                                                                 cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGgtcgtc                                                       <  2:5241839/101‑1 (MQ=18)
                                                                                                                                                 cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc      >  2:1150551/1‑150 (MQ=37)
                                                                                                                                                  tACACCGGGCCGCGCCGCTTTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     <  2:1583675/150‑1 (MQ=255)
                                                                                                                                                  tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCCCCGGCAATAccga                                                                                        >  1:4661626/1‑64 (MQ=17)
                                                                                                                                                  tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCCCCGGCAATAccga                                                                                        <  2:4661626/67‑4 (MQ=17)
                                                                                                                                                  tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  1:3816534/1‑150 (MQ=255)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  1:844907/1‑150 (MQ=25)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  1:2942267/1‑150 (MQ=25)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  2:2821313/1‑150 (MQ=25)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  1:1316952/1‑150 (MQ=37)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  2:3279545/1‑150 (MQ=37)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  2:3250578/1‑150 (MQ=37)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  1:5439978/1‑150 (MQ=37)
                                                                                                                                                   acacTGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcag    <  2:2871145/150‑1 (MQ=255)
                                                                                                                                                    cacaGGGCCGCGCCGCTTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcagg   >  2:4948028/1‑150 (MQ=255)
                                                                                                                                                     acaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                                         >  1:1996294/1‑47 (MQ=255)
                                                                                                                                                     acaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                                         <  2:1996294/47‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGA  >  NC_002947/1544383‑1544679

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