Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,531 C→A 33.8% 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,5310CA33.8% 255.6 / 205.2 204intergenic (+48/+66)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base C (51/84);  new base A (32/37);  total (83/121)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

Only 100 of 216 total aligned reads displayed.

CTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGAT  >  NC_002947/1544382‑1544680
                                                                                                                                                      |                                                                                                                                                      
  gcATTTTGCCCTGGCGGGGGGGGTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAg                                                                                                                                                      <  2:260239/150‑1 (MQ=255)
  gcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAg                                                                                                                                                      <  2:3774990/150‑1 (MQ=255)
     ttttGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACACACGgc                                                                                                                                                   >  2:9010354/1‑150 (MQ=255)
             tGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGGGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCCCTCCTCCACAGGgccgcgccgc                                                                                                                                           >  2:7068778/1‑150 (MQ=255)
             gggcgggGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCggcgcggccct                                                                                                                                          <  2:469816/149‑6 (MQ=37)
                cgggggGGGTGGCGGGGGGGTGGGCGCTGTGGCCGGTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCttt                                                                                                                                        <  1:994249/150‑1 (MQ=255)
                cgggggGGGGGGGGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCttt                                                                                                                                        <  2:594592/150‑1 (MQ=255)
                 ggcggGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGGGGGGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                      <  1:2834314/150‑1 (MQ=255)
                  ggggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGGGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                      <  2:4796676/148‑1 (MQ=255)
                  gggggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                      <  2:6940163/148‑1 (MQ=255)
                  gggggggTGGCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGGGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                      <  1:1378706/148‑1 (MQ=255)
                   gggggtTGGCGGGGGCTTTGGCGCTGTGGCCGCTGTGGTTTGTGTTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                     <  2:4879015/149‑1 (MQ=255)
                   gggggggggCGGGGGGGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                     <  1:6634264/149‑1 (MQ=37)
                   gggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                     <  1:4161152/149‑1 (MQ=255)
                   ggggggTGGGGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                     <  2:645653/149‑1 (MQ=37)
                   ggggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                     <  1:8395874/149‑1 (MQ=37)
                   ggggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                     <  2:645602/149‑1 (MQ=37)
                   ggggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTa                                                                                                                                     <  2:8230474/149‑1 (MQ=37)
                      ggggggggggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                  <  1:8904736/150‑1 (MQ=255)
                      ggggggggggggCGTGGGCGCGGTGGCCGCTGTGGTTGGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                  <  1:3230433/150‑1 (MQ=39)
                      ggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGa                                                                                                                                  <  2:6084821/150‑1 (MQ=37)
                      ggggggCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGa                                                                                                                                  <  2:2804408/150‑1 (MQ=25)
                      gggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGa                                                                                                                                  <  2:8052105/150‑1 (MQ=37)
                      gggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGa                                                                                                                                  <  2:6087888/150‑1 (MQ=37)
                      gggTGGCGGGGGCGTGGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                                  <  1:8163947/150‑1 (MQ=255)
                        ggggCGGGGGCGGTGGCGCTGTGGCCGCTGTGGTGTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCGGTAGCTCGATCATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGc                                                                                                                                <  2:1826847/148‑1 (MQ=38)
                          ggCGGGGGCGTGGGCGCTGTGGCCGCGGTGGTTGGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCgg                                                                                                                              <  1:5214249/150‑1 (MQ=18)
                           gggggggCGTTGGCGCTGGGGCCGGTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg                                                                                                                             <  1:3039681/148‑1 (MQ=255)
                           gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg                                                                                                                             <  2:5784384/150‑1 (MQ=255)
                           gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCggg                                                                                                                             <  2:730847/150‑1 (MQ=37)
                             gggggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGtt                                                                                                                           <  1:9063759/150‑1 (MQ=255)
                                ggggTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTAc                                                                                                                       <  2:8710609/147‑1 (MQ=25)
                                  cGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTAccc                                                                                                                     <  2:2349997/150‑1 (MQ=255)
                                       gcgcTGTGGCCGGTGTGTTTTGTGGTGGGTGGGCGGCCGGGGGATGGGGCTGTAGCTCGATAATCGTTATTGCCGGTGCGGGCCTATTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGATTTGTAGGAGGGGGTTTACCCGCGaa                                                                                                                <  2:3949751/150‑1 (MQ=37)
                                       gcgcTGTGGCCGCTGGGTTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAAGCGGTATGGCCGGGGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTCGGAGCGGTTTTACCCGCGaa                                                                                                                <  2:831967/150‑1 (MQ=255)
                                         gcTGTGGCCGCTGGGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAaga                                                                                                              <  2:4460155/150‑1 (MQ=21)
                                            ggggCCGCTGGGGTTTGGGGTGGGGGGGCGGCGTGGGGACGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                           <  2:7491601/148‑1 (MQ=37)
                                                     gtgGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCCCGGGCCGCGCCGCTTTGTCGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAgg                                                                                                  <  2:5840056/150‑1 (MQ=255)
                                                          tgggggggggtgggCGGCGGGGGCCCGGGGCGTTCGCTCCATACGGGGTTGTGCGGGGGCGGGCCCCTGCGGGGTAAACCCCGCTCCCACACCGGGCCGCGCCGCTTGGTCGGCGGGGTTTACCCCGCGAAGAGGCCAGCACAGGCAATa                                                                                             <  2:7961/143‑1 (MQ=255)
                                                                       gcggcCTGGGTAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGGGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAggg                                                                                <  2:8042417/150‑1 (MQ=255)
                                                                              ggggAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCa                                                                         <  2:6094218/150‑1 (MQ=255)
                                                                                       gCTGTCGCTCGATAATCGTTATTGCGGGGGCCGGCCTCTTCGGGGTTAAACCCGCCCCTACACCGGGCCGCGCCGCTTTGTAGGCGGGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAg                                                                <  1:2183666/150‑1 (MQ=38)
                                                                                         tGTAGCTCGATAATCGGTCTGGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGCCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTg                                                              <  2:2923759/150‑1 (MQ=255)
                                                                                         ggTAGCGCGATAATCGGTATTGCCGGGGCCGGCCTGTTCGGGGGTAAACCCGCTCCTACACCGGGCCGCGCCGCTTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTg                                                              <  2:173815/149‑1 (MQ=37)
                                                                                          gttgcccGATAATCGGTATGGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg                                                             <  1:6843260/144‑1 (MQ=255)
                                                                                            aGCTCGATAATGGGTATTGCGGGGGCCGGCCTCTGCGGGGGAAAACCCGCTCCTACACCGGGCCGCGCCGCTTTGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGgtc                                                           <  1:4517343/150‑1 (MQ=37)
                                                                                                      tCGGTATTGCCGGGGCCGGCCGCTTAGCGGGGAAAAACCCGCCAACACACGGCCGCGCCCCGGTGTAGGGGCGGGGTTAAACCCGCAAAGGGCCCCACAGGCCAAAAAAAAAGAAAGGGCCGCGCAGCCCCCCAGGGGGCGGGCGCCAcc                                                 >  2:4285757/1‑150 (MQ=255)
                                                                                                       cGGTAGGGCGGGGGCCGGCCTCTTCGCGGTTAACCCCGCGCCTACACCGGGCGGCGCGGCTTTTTCGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCg                                                <  2:3081815/150‑1 (MQ=37)
                                                                                                            ttGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                              <  1:6995544/68‑1 (MQ=255)
                                                                                                                gggggCCGGCGTGTGGGCGGTAAACCCCGCCCCACCACGGGGCCGCGCGGTTTTGTGGGCGGGTGTTACCCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTccccc                                       <  1:111781/149‑1 (MQ=255)
                                                                                                                  gggcgggcCTCTGCGCGGGAACCCCCGCTCCCCCCCGGGGCCGCGCGGCTTTGTCGGCGGGGTTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgc                                     <  1:6899266/150‑1 (MQ=255)
                                                                                                                   ggccggccTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgcc                                    <  1:5282364/150‑1 (MQ=255)
                                                                                                                       ggccTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGc                                >  2:1496973/1‑150 (MQ=255)
                                                                                                                           tctTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTcc                            >  2:8788467/1‑150 (MQ=255)
                                                                                                                           tctTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCCGTGTAGGAGCGGGTTTACCCGCGAAGTGGCCAGCACCGGCCCTACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGGCGTCCGCCACCGGCGTCCCCCGCCCGGCTGcc                            >  2:6942133/1‑150 (MQ=18)
                                                                                                                            ctTCGCGGGTAAACCCGCTCCAACACAGGGCCGCGCCGCTTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCt                           <  1:6795525/150‑1 (MQ=255)
                                                                                                                                  gggTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGcc                     >  1:5651013/1‑150 (MQ=255)
                                                                                                                                  gggTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGcc                     >  1:4784505/1‑150 (MQ=255)
                                                                                                                                  gggTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGcc                     >  2:1385192/1‑150 (MQ=255)
                                                                                                                                  gggTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGcc                     >  1:707562/1‑150 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCAc                   >  2:3409579/1‑150 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCAc                   >  1:3674458/1‑150 (MQ=255)
                                                                                                                                     tAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTcc                                          <  2:5044782/126‑1 (MQ=255)
                                                                                                                                     tAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGGAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTcc                                          >  1:5044782/1‑126 (MQ=255)
                                                                                                                                     tAAACCCGCACCTAAACCGGGCCGGGCCGCTTTGTCGGCGGGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACt                  <  1:469265/150‑1 (MQ=255)
                                                                                                                                      tacccccGCCCCAACACGGGGCCGGGCCGTTTTGTCGGTGGGGTTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACt                  <  1:291510/146‑1 (MQ=255)
                                                                                                                                        aCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTCCa               >  2:5449810/1‑150 (MQ=255)
                                                                                                                                        aCCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGcc               >  2:4348079/1‑149 (MQ=25)
                                                                                                                                          cccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg              >  1:2356849/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg              >  2:2795872/1‑150 (MQ=37)
                                                                                                                                           ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTg                                                              <  2:7579414/101‑1 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTg                                                              >  1:7579414/1‑101 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg             >  1:4848397/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg             >  1:6284329/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg             >  2:2059560/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTcc                                                       <  2:5252009/108‑1 (MQ=18)
                                                                                                                                             gCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc           <  2:1384345/150‑1 (MQ=255)
                                                                                                                                              cTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATCGCCCTGCAgtggc          <  1:5245470/150‑5 (MQ=255)
                                                                                                                                               tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga         >  2:8084920/1‑150 (MQ=255)
                                                                                                                                               tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga         >  2:612276/1‑150 (MQ=37)
                                                                                                                                               tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga         >  1:3200583/1‑150 (MQ=37)
                                                                                                                                                ccTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCGGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat        >  2:1555879/1‑150 (MQ=37)
                                                                                                                                                 cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc       >  1:5570002/1‑150 (MQ=255)
                                                                                                                                                 cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc       >  1:4766651/1‑150 (MQ=25)
                                                                                                                                                 cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCGTCAGGGCCGCTCAGCCCCCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc       >  2:8159283/1‑150 (MQ=37)
                                                                                                                                                  tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATACAAAGCTTCAGGGCCGCTCa                                                                         <  1:4750031/83‑1 (MQ=255)
                                                                                                                                                  tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  1:4375802/1‑150 (MQ=255)
                                                                                                                                                  tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  2:1208209/1‑150 (MQ=255)
                                                                                                                                                  tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  1:1674071/1‑150 (MQ=255)
                                                                                                                                                  tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  1:1042127/1‑150 (MQ=255)
                                                                                                                                                  tACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  1:669422/1‑150 (MQ=255)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  1:3856175/1‑150 (MQ=37)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  2:1935159/1‑150 (MQ=37)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      >  1:7840892/1‑150 (MQ=37)
                                                                                                                                                  tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca      <  1:612276/150‑1 (MQ=37)
                                                                                                                                                   acacAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcag     >  1:8945344/1‑150 (MQ=255)
                                                                                                                                                   acacAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcag     >  2:3643556/1‑150 (MQ=255)
                                                                                                                                                    cacaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcagg    >  2:2681462/1‑150 (MQ=255)
                                                                                                                                                     ccaGGGCCGCGCCGCATTGCAGGAGCGGGATTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcagga   <  2:5297410/149‑1 (MQ=255)
                                                                                                                                                      caGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATCGCCCTGCAGTGGcatcaggaa  <  1:9010354/150‑2 (MQ=38)
                                                                                                                                                      |                                                                                                                                                      
CTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGAT  >  NC_002947/1544382‑1544680

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: