Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,549,480 C→A 35.0% 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,549,4800CA35.0% 204.5 / 190.8 180intergenic (+48/+66)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base C (42/75);  new base A (34/29);  total (76/104)
Fisher's exact test for biased strand distribution p-value = 2.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

Only 100 of 196 total aligned reads displayed.

CTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGA  >  NC_002947/1549331‑1549628
                                                                                                                                                     |                                                                                                                                                    
 tgcATTTTGCCGTGGCGGGGGGGGTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                     <  2:107081/150‑1 (MQ=255)
 tgcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                     >  2:3528241/1‑150 (MQ=255)
 tgcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCGTCGCGGGTAAAACCGCTCCGAaaca                                                                                                                                                     >  1:3703388/1‑150 (MQ=255)
 tgcATTTTGCCCGGGGGGGGGGGGGGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                     <  2:107117/150‑1 (MQ=39)
        tcccgtGGCGCCGGGGTTGGCGGGGGCTTTGGCGCTTTGGCCGATGTGTTTTGTGGTGGGTGGCCGCCATGGGCAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgc                                                                                                                                              <  1:2573760/145‑1 (MQ=37)
        tGCCCTGGGGGGGGGGGTGGGGGGGGCGTTGGCGCGGTGGCCGCGGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgc                                                                                                                                              <  2:2680358/150‑1 (MQ=39)
        tGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgc                                                                                                                                              <  1:668143/150‑1 (MQ=255)
         gCCCGGGCGGGGGGGGTGGCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTAGTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgcg                                                                                                                                             <  2:4613811/150‑1 (MQ=255)
             gggcggcGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCggcgcggccc                                                                                                                                         <  1:569047/149‑5 (MQ=255)
                gggggggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCttt                                                                                                                                      <  1:2270347/146‑1 (MQ=255)
                cgggggGGGGGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCggcgcggccctgt                                                                                                                                      <  2:2016110/150‑8 (MQ=37)
                 ggcggGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGTTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTg                                                                                                                                     <  2:5058679/150‑1 (MQ=255)
                 ggcggGGGGGGGGGGGGCGTTGGCGGTGGGGCGGCTGTGGTTTTTGGTGGTGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTg                                                                                                                                     <  1:888186/150‑1 (MQ=38)
                  gggggggTGGGGGGGGCGTTGGCGTTGGGGCCGGTGTGTTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTgtgt                                                                                                                                    <  1:2560029/148‑1 (MQ=25)
                  gcggGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                    <  1:2955966/150‑1 (MQ=255)
                   gggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                   <  1:4843148/149‑1 (MQ=25)
                   ggggggTGGGGGGGGGGTTGGCGGTGTGGCCGGTGGGGTTTGGGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                   <  1:1964288/149‑1 (MQ=25)
                   ggggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                   <  2:5400456/149‑1 (MQ=255)
                   ggggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                   <  1:4616563/149‑1 (MQ=255)
                   ggggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                   <  1:1030334/149‑1 (MQ=37)
                   cggGGGGGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                   <  2:1464420/150‑1 (MQ=37)
                    ggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                  <  2:5581356/150‑1 (MQ=255)
                    gggggTGGGGGGGGCGTTGGCGCTGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGACGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                  <  2:1222032/150‑1 (MQ=255)
                    gggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAg                                                                                                                                  <  2:1331640/150‑1 (MQ=37)
                    gggggTGGCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                  <  1:4150261/150‑1 (MQ=255)
                    gggggTGGCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAg                                                                                                                                  <  2:2392284/150‑1 (MQ=37)
                    gggggTGGCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAg                                                                                                                                  <  2:917109/150‑1 (MQ=37)
                      gggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                <  2:729228/150‑1 (MQ=255)
                      gggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGa                                                                                                                                <  2:1758721/150‑1 (MQ=37)
                      gggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGa                                                                                                                                <  1:1795838/150‑1 (MQ=37)
                        gTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTTAAACCCCGGCGAAAAAGGGCCGCGCCGCGTTggatgtgc                                                                                                                              >  2:2125122/1‑143 (MQ=255)
                        gTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGc                                                                                                                              <  2:1269847/150‑1 (MQ=255)
                         tGGCGGGGGCGTTGGCGCTGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTAGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCg                                                                                                                             <  1:456714/150‑1 (MQ=255)
                          ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCttt                                                                                                                                      <  1:1906399/140‑1 (MQ=255)
                          ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                            <  2:5233222/150‑1 (MQ=255)
                          ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCgg                                                                                                                            <  1:1027405/150‑1 (MQ=37)
                           gCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTGTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGGAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCggg                                                                                                                           <  2:4096126/150‑1 (MQ=25)
                            ggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGt                                                                                                                          <  2:323756/149‑1 (MQ=255)
                                ggggtgggCGCTGTGGCCGCTGTGGTTTTTGTGGGGGGGGGGGCCTTGGGACGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTAc                                                                                                                      <  1:160993/144‑1 (MQ=37)
                                ggggTTGGCGCTGTGGCCGCTGTGGTTTGGGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTAc                                                                                                                      <  1:2366839/147‑1 (MQ=255)
                                 gCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTAcc                                                                                                                     <  1:3528241/150‑1 (MQ=255)
                                  cGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTAccc                                                                                                                    <  1:1709410/150‑1 (MQ=255)
                                   gTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCCCGACAACACCGGGCCGCGCCGCGTTGTAGGAGAGGGTTTAcccc                                                                                                                   >  2:2301069/1‑149 (MQ=255)
                                         ggtgtgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAaga                                                                                                             <  2:2233500/148‑1 (MQ=255)
                                          cgggggCCGCTGTGGTTGGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                            <  2:3643375/146‑1 (MQ=37)
                                           tgtgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTgtgt                                                                                                                                    <  2:4533470/125‑1 (MQ=25)
                                           tgtgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTgtgt                                                                                                                                    >  1:4533470/1‑125 (MQ=25)
                                           tgtgGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCACTTAGCGGGTAAACCCCCTCCAACACAGGGCCGCCCCGCGTTGTAGGAGCGGGTTTACCCCCGCACAgg                                                                                                           >  2:4118744/1‑150 (MQ=255)
                                               gCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGa                                                                                                                <  1:3063223/141‑1 (MQ=255)
                                               gCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCGTCGCGGGTAAAACCCCACCAAAACAGGGCCGCCCCCCTGTTTAGGTGCCGGTTTACCCCCGAAGAGGCCac                                                                                                       >  2:2727161/1‑149 (MQ=255)
                                                 cGCTGTGGTTTGTGGTGGGGGGGCGGCGTGGGGATGGGGCTGTCGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACACCGGGCCGCGCCGCTTTGTAGGTGCGGGTTTACCCGCGAAGAGGCCAGca                                                                                                     <  2:838205/150‑1 (MQ=37)
                                                      ggttttGTGGGGGGGGGGCGGCCTGGGGATGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTGCGCGGGTAACCCCGCTCCTACACCGGGCCGCGCCGCTTTGTCGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGGc                                                                                                <  2:1718070/147‑1 (MQ=38)
                                                          ttGTGGGGGGGGGGCGGCGGGGGGACGGGGCTGTTGCTCGATAATCGGTATTGCGGGGGCGGGCCTCTGCGCGGGTAACCCCGCTCCTACACCGGGCCGCGCCGCTTTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATa                                                                                            <  1:4118744/150‑1 (MQ=37)
                                                               gtgggtggGCGGCGGGGGGACGGGGCTGTAGCTCTATAATCGTTATTGCCGGGGCCGGCCTCTTCGGGGTTAAACCCGCTCCTACACTGGGCCGCGCCGCTTTGTCGGTGGGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAg                                                                                       <  1:5528736/150‑1 (MQ=255)
                                                                      gggggcCTGGGGCCGGGGCTGTCGCTCGATACTGGGTATTGCGGGGGCCGGCCTCTTCGCGGGTAACCCCGCTCCTCCACAGGGCCGCGCCGCTTTGTCGGCGGGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAgg                                                                                <  1:5057803/147‑1 (MQ=37)
                                                                                 gAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAAAAAAAAGCTTCAGGGCCGCGCAGcc                                                                     >  2:4660072/1‑150 (MQ=255)
                                                                                    ggggCTGTAGCTCTATATGGGGTTTGGCGGGGGCCGGCCTCTGGGCGGTTAACCCCGCTCCCACACGGGCCCGCGCCGCTTTGTCGGAGGGTTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCAcc                                                                  <  1:3996683/150‑1 (MQ=255)
                                                                                        cTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCGCCTACACAGGGCGGCGCCGCTTTGTAGGCGCGGGTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGt                                                              <  2:2575843/150‑1 (MQ=255)
                                                                                          gTAGCTCGATAATCGGTATGGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg                                                            <  2:4560466/150‑1 (MQ=255)
                                                                                                 gATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg                                                                                                                           <  1:4280014/80‑1 (MQ=255)
                                                                                                 gATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg                                                                                                                           >  2:4280014/1‑80 (MQ=255)
                                                                                                            tggcgggggCCGGCCTCTTCGGGGTTAAACCCGCTCCTACACTGGGCCGCGCGGCTTTGTAGGAGGGTTTTTCCCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTc                                          <  1:4437569/145‑1 (MQ=255)
                                                                                                              gCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTccc                                        <  2:3694903/150‑1 (MQ=255)
                                                                                                              gCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTccc                                        <  2:3695340/150‑1 (MQ=255)
                                                                                                                    gccggccTCTGCGCGGGAAACCCCCCTCCCCCCCCGGGCCGCGCCGCTTTGTCGGCGGGGGTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgccc                                  <  1:3794062/150‑1 (MQ=255)
                                                                                                                        gccTCTGCGCGGGTAACCCCGCTCCTACCCAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCt                              <  2:250096/150‑1 (MQ=255)
                                                                                                                          ctctTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc                            >  2:5529785/1‑150 (MQ=255)
                                                                                                                           cctgcgcgGGTAAACCGGCTCCTACACCGGGCCGCGCCGCTTGGTCGGAGCGGGTTTACCCGCGACGAGGCCAGCACAGGCAATACAACGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTcc                           <  2:4299833/146‑1 (MQ=255)
                                                                                                                                  gggTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACCAAGCTTCCGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGcc                    >  2:579433/1‑150 (MQ=255)
                                                                                                                                  gggTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGcc                    >  1:26991/1‑150 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCg                                               <  2:2938466/121‑1 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGa                        >  1:1248222/1‑144 (MQ=255)
                                                                                                                                        aCCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCa              >  2:3097189/1‑150 (MQ=25)
                                                                                                                                          ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg            >  1:5124655/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg            >  1:2935173/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGa                                                                <  2:156384/98‑1 (MQ=17)
                                                                                                                                          ccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg            >  2:390439/1‑150 (MQ=37)
                                                                                                                                          ccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg            >  1:1830950/1‑150 (MQ=37)
                                                                                                                                          ccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg            >  2:130764/1‑150 (MQ=37)
                                                                                                                                            gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc          >  2:653209/1‑150 (MQ=37)
                                                                                                                                            gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc          >  2:2259507/1‑150 (MQ=37)
                                                                                                                                              tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga        >  2:4829743/1‑150 (MQ=255)
                                                                                                                                              tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga        >  2:4584503/1‑150 (MQ=21)
                                                                                                                                              tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga        >  1:53535/1‑150 (MQ=37)
                                                                                                                                               ccTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat       >  2:2823652/1‑150 (MQ=37)
                                                                                                                                                cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGc                               >  1:3239572/1‑125 (MQ=255)
                                                                                                                                                cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGc                               <  2:3239572/125‑1 (MQ=255)
                                                                                                                                                cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc      >  1:4771674/1‑150 (MQ=255)
                                                                                                                                                cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc      >  2:405647/1‑150 (MQ=255)
                                                                                                                                                cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc      >  2:3375534/1‑150 (MQ=37)
                                                                                                                                                 tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  2:4385797/1‑150 (MQ=255)
                                                                                                                                                 tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  1:620511/1‑150 (MQ=255)
                                                                                                                                                 tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  2:1550220/1‑150 (MQ=255)
                                                                                                                                                 tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  2:1812327/1‑150 (MQ=255)
                                                                                                                                                 tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  2:1105937/1‑150 (MQ=37)
                                                                                                                                                 tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  1:4433265/1‑150 (MQ=37)
                                                                                                                                                 tACAAAGCGGCGCGGCCCTGTGTAGGAACGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca     >  2:4753756/1‑150 (MQ=37)
                                                                                                                                                   cacaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcagg   >  1:4965266/1‑150 (MQ=255)
                                                                                                                                                   cacaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcagg   >  2:990299/1‑150 (MQ=255)
                                                                                                                                                    acaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcagga  >  2:3404055/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGA  >  NC_002947/1549331‑1549628

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: