Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,549,480 C→A 34.7% 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,4800CA34.7% 294.0 / 263.6 251intergenic (+48/+66)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base C (51/112);  new base A (46/41);  total (97/154)
Fisher's exact test for biased strand distribution p-value = 1.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01

Only 100 of 268 total aligned reads displayed.

CTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGAT  >  NC_002947/1549331‑1549629
                                                                                                                                                      |                                                                                                                                                      
 tgcTTTTTGCCCGGGGGGGGGGGGTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTGGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                       <  2:7525387/150‑1 (MQ=39)
 tgcTTTTGGCCGTGGCGGGGGGGGTGGCGGGGGCGTTGGCGCGGTGGCCGCTGTGGTTGGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                       <  2:7483257/150‑1 (MQ=39)
 tgcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCCCACCGacacc                                                                                                                                                       >  2:2606434/1‑149 (MQ=255)
  gcATTTGGCCCGGGGGGGGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAg                                                                                                                                                      <  1:7147601/150‑1 (MQ=255)
   cATTTTGCCCTGGCGGGGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAgg                                                                                                                                                     <  2:77146/150‑1 (MQ=255)
    aTTTTGCCCGGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAggg                                                                                                                                                    <  2:6603018/150‑1 (MQ=255)
                ggggggggggggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCttt                                                                                                                                       <  1:4725374/146‑1 (MQ=255)
                gggcggGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCttt                                                                                                                                       <  2:5764647/149‑1 (MQ=255)
                 cggggggggggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTg                                                                                                                                      <  1:7032499/147‑1 (MQ=255)
                  gggggggTGGGGGGGGCGTTGGCGCTGTGGCCGCCGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTgtgt                                                                                                                                     <  2:6706416/148‑1 (MQ=37)
                  gcggGGGTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                     <  2:2095950/150‑1 (MQ=255)
                  gcggGGGGGGCGGGGGCGTTGGCGCTGTGGCCGCTGGGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                     <  2:6491560/150‑1 (MQ=255)
                   gggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                    <  2:2185565/149‑1 (MQ=37)
                   gggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                    <  1:5879270/149‑1 (MQ=37)
                   ggggggTGGGGGGGGCGTTGGCGCTGTGGCGGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                    <  1:5613486/149‑1 (MQ=21)
                   ggggggTGGGGGGGGCGTTGGCGCTGTGGCCGGTGTGTTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                    <  2:8597654/149‑1 (MQ=255)
                   ggggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                    <  1:1154151/149‑1 (MQ=37)
                   ggggggTGGGGGGGGCGTTGGCGCTGGGGCCGCTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                    <  1:3295449/149‑1 (MQ=25)
                   ggggggTGGCGGGGGCTTTGGCGCTGTGGCCGTTGTGTTTTGTGGTGGTTGGGCGGCCTGGGTAAGGGGCTGTAGCTCGATAATCGGTATTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                    <  1:3919077/149‑1 (MQ=18)
                   ggggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGGGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                    <  1:4405836/149‑1 (MQ=255)
                   cggGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                    <  2:5479815/150‑1 (MQ=37)
                    ggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                   <  1:8522982/150‑1 (MQ=255)
                    ggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                   <  2:2693795/150‑1 (MQ=255)
                    ggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                   <  2:292008/150‑1 (MQ=255)
                    gggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                   <  1:874615/150‑1 (MQ=255)
                    gggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                   <  2:8682200/150‑1 (MQ=255)
                    gggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAg                                                                                                                                   <  2:6791301/150‑1 (MQ=37)
                    gggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                   <  2:6049780/150‑1 (MQ=255)
                    gggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCGTGGCGGGGAAACCCGCGCCGACACAGGgcccccccgccgtgggg                                                                                                                                    >  2:3364466/1‑143 (MQ=255)
                      ggggggggggggCGTTGGCGCTGTGGCCGCGGTGGTTTGGGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                 <  1:7626327/150‑1 (MQ=255)
                      gggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                 <  2:2832923/150‑1 (MQ=255)
                      gggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                 <  2:4488830/150‑1 (MQ=255)
                          ggggggggCGTTGGGGCTGTGGCCGCGGGGGTGGGTGGGGGGTGGGGGGCCTGGGGAAGGGGCTGTAGCTCGAGAATCGGTAGGGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCgg                                                                                                                             <  2:7803077/147‑1 (MQ=18)
                          ggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCgg                                                                                                                             <  1:2804132/147‑1 (MQ=25)
                          ggggggggCGTGGGGGCTGTGGCCGGTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCgg                                                                                                                             <  2:7389925/147‑1 (MQ=25)
                          ggCGGGGGCTTGGGCGGTGTGGCCGGTGTGGTTTGTGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                             <  1:2438516/150‑1 (MQ=38)
                          ggCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCgg                                                                                                                             <  2:777286/150‑1 (MQ=25)
                           gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCggg                                                                                                                            <  2:1272428/150‑1 (MQ=25)
                            ggggggCGTTGGCGCTGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGt                                                                                                                           <  1:7349680/149‑1 (MQ=255)
                            ggggggCGTTGGCGCTGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGt                                                                                                                           <  2:7990294/149‑1 (MQ=255)
                             gggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCggg                                                                                                                            <  2:427421/148‑1 (MQ=25)
                             gggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCggg                                                                                                                            >  1:427421/1‑148 (MQ=25)
                               gggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCCCGCCGCATTGTAGGAGCGGGTTTa                                                                                                                        >  2:8756056/1‑150 (MQ=255)
                                ggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCggaaatca                                                                                                                       >  2:7424193/1‑144 (MQ=25)
                                    tgggCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCcgc                                                                                                                   <  2:2804509/148‑1 (MQ=255)
                                     gggCGCTGGGGCCGCTGTGGTTTGGGGTGGGTGGGCGGCCTGGGGAAGGGGCGGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCACAGGGCCGCGCCGCTTTGTCGGAGCGGGTTTACCcgcg                                                                                                                  <  2:5827652/149‑1 (MQ=38)
                                      ggCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTa                                                                                                                        >  1:5988829/1‑143 (MQ=255)
                                          cTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                             <  2:5640673/150‑1 (MQ=255)
                                          cTGGGGCCGCTGTGGTTTGGGGGGGGTGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCATACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                             <  2:3513888/150‑1 (MQ=37)
                                              ggCCGCGGTGGTGGGGGGGGGGGGGGCGGCGTGGGGCCGGGGCGGGAGCTCGATAAGCGGTATGGCCGGGGCCGGCCTCTGCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCGGCTTTGTAGGAGCGTGTTTACCCGCGAAGAGGCCa                                                                                                         <  2:4615913/150‑1 (MQ=21)
                                                      gggTTTGGGGGGGGGGGGCGGCGTGGGGACGGGGCTGTCGCTCGATAATCGGTATGGCCGGGGCCGGCCTCTTCGCGGGTACACCCGCTCCTACACCGGGCCGCGCCGCTTGGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGGc                                                                                                 <  2:4166525/149‑1 (MQ=25)
                                                        gTTTGTGGTGGGGGGGGGGCCTGGGTACGGGGCTGACGCGCGATACCCGTTATTGCGGGGGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACACCGGGCCGCGCCGCTTTTTCGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCaa                                                                                               <  1:6502254/150‑1 (MQ=21)
                                                         tttGGGGGGGGGGGGGGGCCGGGGTACGGGGCTGTAGCTCTATAATGGTTATGGCCGGGGCCGGCCTCTTGGCGGTTAACCCCGCACCTACACCGGGCCGCGCCGCTTTGTCGGCGGGGTTTTACCCGCGAAGAGGCCAGCACAGGCAAt                                                                                              <  1:3552061/150‑1 (MQ=255)
                                                                 gggtgggCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCGCGCCAACACAGGGCCGCGCCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAAAACAAAGca                                                                                      >  2:3596796/1‑149 (MQ=38)
                                                                    ggggCGGCCTGGGGACGGGGCTGTAGCTCGATACTCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCGCCTACCCAGGGCCGCGCCGCTTTGTAGGTGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCa                                                                                   <  2:5931059/149‑1 (MQ=37)
                                                                    ggggCGGCCGGGGGACGGGGCTTTAGCTCTATAAGGGGTATTGCGGGGGCCGGCCTCTGCGGGGGTAAACCCGCTCCTACACTGGGCCGCGCCGCTTTGTTGGTGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCa                                                                                   <  1:728706/149‑1 (MQ=25)
                                                                     gggCGGCGGGGGGACGGGGCTGTCGCTCGATAATGGGTTTTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGGCCGCGCCGCTTTGTCGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAg                                                                                  <  2:3921352/150‑1 (MQ=25)
                                                                              ggggACGGGGCTGTAGCTCGATAATCGTTTTTGCGGGGGCCGGCCTCTTCGCGGGAAACCCCGCCCCTCCACAGGGCCGCGCCGCTTTTTCGGCGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCa                                                                         <  1:4245677/150‑1 (MQ=255)
                                                                                      cgttttAGCTCGATATGCGGTATGGCCGGGGCCGGCAACTTGGCGTGAAACCCCCCCCCTACACCGGGCCGCGCGGATTTGTTGGCGGGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGa                                                                 <  1:5326095/145‑1 (MQ=255)
                                                                                         tGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTg                                                              <  2:6446137/150‑1 (MQ=255)
                                                                                                            gggcgggggCCGCCCTCTGGGGGGCAAACCCCGCCCCTACACCGGGCCGCGCCGTTTTGTCGGGGGGGGATCCCCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTc                                           <  1:1743439/145‑1 (MQ=255)
                                                                                                                cGGGGCCGGCCTCTTCGCGGGGTAAACCCCGCCAAAAAAAGGGCCCGCCCCGTTTGTGGTGCCGGGTTTAACCCCAAGAGGCCACCACACGCAAAACAAATCATCAGGGGCCCACAGCCCCCCCCTGGTCGTCCGCCACCCGCCTccccc                                       >  2:7108283/1‑150 (MQ=255)
                                                                                                                   ggccggccTCTCCGCGGGTAAACCCGCTCCTACACCGGGCCGCGCCGCTTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgcc                                    <  2:3677092/150‑1 (MQ=255)
                                                                                                                       ggccTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCTCTTTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGc                                <  2:3324791/150‑1 (MQ=255)
                                                                                                                          ctctTCGCGGGTAAACCCGCTCCGACACAGGGCCGCGCCGCGTTGTAGGAGCGGGGTTACCCGCGAAGGGGCCACCACAGGCAATACCAAGCGTCAGGGCCGCTCAGCCCCCGCGGGGGCGGCCGCCCCCGGCGTCCCCCCCCCGGCTTc                             >  2:2387542/1‑150 (MQ=25)
                                                                                                                          ctctGCGCGGGTAAACCCGCTCCTACCCTGGGCCGCGCCGCTTTGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc                             <  1:4989289/150‑1 (MQ=255)
                                                                                                                             ttCGCGGGTAAACCCCCTCCGACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTAACCGCCAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTccgg                          >  2:1529439/1‑148 (MQ=21)
                                                                                                                               cgcgGGTACACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGAt                        <  2:2575640/150‑1 (MQ=255)
                                                                                                                                gcgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGCGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATg                       >  2:8194824/1‑150 (MQ=255)
                                                                                                                                gcgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCg                                                      >  1:2827217/1‑119 (MQ=255)
                                                                                                                                gcgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATg                       >  1:2177349/1‑150 (MQ=255)
                                                                                                                                 cgGGTAAACCCGCTCCTACCCCGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGc                      <  2:8240593/150‑1 (MQ=255)
                                                                                                                                    ggTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCa                    >  2:6891837/1‑150 (MQ=255)
                                                                                                                                     gTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCAc                   >  1:1660924/1‑150 (MQ=25)
                                                                                                                                      tAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACt                  >  1:1226222/1‑150 (MQ=255)
                                                                                                                                         aCCCGCCCCTACACCGGGCCGCGCCGTTTTGTCGGCGGGGGTTCCCCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCa               <  1:4868038/150‑1 (MQ=255)
                                                                                                                                          cccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg              >  2:1789736/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg              >  1:8428029/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg              >  1:5143683/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTACAAAGCGGCGCGGCCCTGTGGAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGccg              >  2:1320714/1‑148 (MQ=25)
                                                                                                                                           ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg             <  2:5106350/150‑1 (MQ=255)
                                                                                                                                           ccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg             >  1:3459427/1‑150 (MQ=25)
                                                                                                                                           ccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg             >  1:6895197/1‑150 (MQ=37)
                                                                                                                                           ccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg             >  2:1283694/1‑150 (MQ=37)
                                                                                                                                           ccGCTCCAACACCGGGCCGCGCCGCTTTGGCGGCGGGGTTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg             <  1:3768294/150‑1 (MQ=255)
                                                                                                                                             gCTCCTACCCCGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc           <  2:2270193/150‑1 (MQ=255)
                                                                                                                                             gCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc           >  1:6385906/1‑150 (MQ=255)
                                                                                                                                              cTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCg          >  2:5392988/1‑150 (MQ=37)
                                                                                                                                               tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga         >  2:6154398/1‑150 (MQ=255)
                                                                                                                                               tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga         >  1:6761799/1‑150 (MQ=37)
                                                                                                                                               tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga         >  1:7945548/1‑150 (MQ=37)
                                                                                                                                                 cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc       >  1:8229980/1‑150 (MQ=255)
                                                                                                                                                 cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc       >  2:6227979/1‑150 (MQ=255)
                                                                                                                                                 cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc       >  1:981095/1‑150 (MQ=25)
                                                                                                                                                 cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc       >  1:7436468/1‑150 (MQ=37)
                                                                                                                                                 cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGACCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc       >  2:8296253/1‑150 (MQ=37)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  2:8517923/1‑150 (MQ=37)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  1:4023230/1‑150 (MQ=37)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  2:5266552/1‑150 (MQ=37)
                                                                                                                                                     acaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcagga   >  2:4097905/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGAT  >  NC_002947/1549331‑1549629

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: