Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 313,852:1 +A 5.4% intergenic (+23/+28) ygaU → / ← yrfG murein cross‑linking regulator/purine nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947313,8521.A 5.4% 463.2 / 18.2 111intergenic (+23/+28)ygaU/yrfGmurein cross‑linking regulator/purine nucleotidase
Reads supporting (aligned to +/- strand):  ref base . (52/53);  new base A (3/3);  total (55/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 112 total aligned reads displayed.

AGCCATTGCCATGGTCGTGTATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGT  >  NC_002947/313706‑314000
                                                                                                                                                   |                                                                                                                                                    
aGCCATTGCCATGGTCGTGTATGGGAATGCAAATCATTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACAACCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATAGGTGAAGTCTGTACAGGCCTCTTc                                                                                                                                                   >  1:467281/1‑150 (MQ=255)
    aTTGCCATGGTCGTGTATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCgg                                                                                                                                               <  1:1156170/150‑1 (MQ=255)
             gTCGTGTATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCcg                                                                                                                                      >  2:208978/1‑150 (MQ=255)
                  gtATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCcgc                                                                                                                                     <  1:617390/146‑1 (MQ=255)
                  gtATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCcgc                                                                                                                                     >  2:617390/1‑146 (MQ=255)
                  gtATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGaag                                                                                                                                 >  2:134065/1‑149 (MQ=255)
                  gtATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATACTAAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTATACAGGCCTCTTCGCAAGCGCGCCCACAAAt                                                                                                                                 >  2:348780/1‑150 (MQ=255)
                   tATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGaaga                                                                                                                                <  1:134065/150‑3 (MQ=255)
                   tATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATg                                                                                                                                >  1:505101/1‑150 (MQ=255)
                   tATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACCGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGaaga                                                                                                                                <  1:1665034/150‑3 (MQ=255)
                     tGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGGGTATTCCTGACTGATCGCTAAAGTCTGTACAGGCCTCTTCGCGTGCTCCCCCGGGAATggg                                                                                                                              <  2:349483/150‑1 (MQ=255)
                         aaTGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCt                                                                                                                          >  2:2037127/1‑150 (MQ=255)
                             ccAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACa                                                                                                                      <  1:2037127/150‑1 (MQ=255)
                             ccAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTAGCGGGCGCGCCAGCGAATGGGTCCTTACa                                                                                                                      <  2:1354525/150‑1 (MQ=255)
                              cAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGgcgc                                                                                                                                           >  1:2162096/1‑128 (MQ=255)
                              cAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAg                                                                                                                     >  1:51366/1‑150 (MQ=255)
                              cAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAg                                                                                                                     <  1:2202627/150‑1 (MQ=255)
                              cAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAg                                                                                                                     <  2:51366/150‑1 (MQ=255)
                              cAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAg                                                                                                                     >  2:2265058/1‑150 (MQ=255)
                               aaTCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAgg                                                                                                                    <  1:943116/150‑1 (MQ=255)
                               aaTCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCCCATTCGCGGGCGCGCCCGCGaagag                                                                                                                               <  1:1755358/140‑4 (MQ=255)
                               aaTCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCCCATTCGCGGGCGCGCCCGCGaagag                                                                                                                               >  2:1755358/1‑137 (MQ=255)
                                aTCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCggg                                                                                                                                              >  2:490340/1‑123 (MQ=255)
                                aTCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCggg                                                                                                                                              <  1:490340/123‑1 (MQ=255)
                                aTCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGc                                                                                                                   <  2:301019/150‑1 (MQ=255)
                                aTCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGc                                                                                                                   <  2:2306618/150‑1 (MQ=255)
                                aTCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCATATTCCTGACTGAGCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGc                                                                                                                   <  1:1691146/150‑1 (MQ=255)
                                    gTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATgg                                                                                                                               <  1:241699/134‑1 (MQ=255)
                                    gTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATgg                                                                                                                               >  2:241699/1‑134 (MQ=255)
                                    gTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTAc                                                                                                               <  2:1566977/150‑1 (MQ=255)
                                        aacaaGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGaa                                                                                                                                  <  2:394552/127‑1 (MQ=255)
                                        aacaaGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGaa                                                                                                                                  >  1:394552/1‑127 (MQ=255)
                                        aacaaGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGt                                                                                                           >  1:312580/1‑150 (MQ=255)
                                          caaTATTTTCTCGGCCAACAATGCTATGCTCAAACGCCCCGACAAGCTCTATCCGGGTCATGTGCTGCGTATTCCTGACTGATCTCTGAAGTCTTTACAGGCCTCTTCACGTGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTaa                                                                                                         <  2:467281/150‑1 (MQ=255)
                                                tttCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTTACTCCCg                                                                                                   <  1:1461003/150‑1 (MQ=255)
                                                tttCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCg                                                                                                   <  1:208978/150‑1 (MQ=255)
                                                 ttCGAGGCAAACAAGCCGATGCTCAAACACCTGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCgg                                                                                                  <  1:348780/150‑1 (MQ=255)
                                                    gAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTaa                                                                                               <  1:1235600/150‑1 (MQ=255)
                                                       gCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATc                                                                                             >  2:739566/1‑149 (MQ=255)
                                                         gaacaaGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTc                                                                                          >  2:1649904/2‑150 (MQ=255)
                                                         caacaaGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTc                                                                                          <  1:1384897/150‑1 (MQ=255)
                                                          aacaaGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCgg                                                                                                  >  2:162463/1‑141 (MQ=255)
                                                          aacaaGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCgg                                                                                                  <  1:162463/141‑1 (MQ=255)
                                                           acaaGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGa                                                                                        <  1:2330711/150‑1 (MQ=255)
                                                             aaGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATgg                                                                                                                               >  1:1559720/1‑109 (MQ=255)
                                                                  gATGGTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGCATGGGGCCTTATTGGCCTACCCGTAACTCCCGGTAATCCTCGACcgccgc                                                                                 >  2:535860/1‑150 (MQ=255)
                                                                         aaaCACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTAc                                                                                                               >  1:248640/1‑113 (MQ=255)
                                                                         aaaCACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCt                                                                                           <  1:1426816/133‑1 (MQ=255)
                                                                         aaaCACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCt                                                                                           >  2:1426816/1‑133 (MQ=255)
                                                                         aaaCACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTc                                                                                          >  2:514279/1‑134 (MQ=255)
                                                                         aaaCACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTc                                                                                          <  1:514279/134‑1 (MQ=255)
                                                                         aaaCACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAActct                                                                          >  2:1093137/1‑150 (MQ=255)
                                                                          aaCACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCtt                                                                         <  1:520072/150‑1 (MQ=255)
                                                                          aaCACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTAAGGCCCCATTCGCGGGCGCGCCcgc                                                                                                                                     >  1:1098821/1‑90 (MQ=37)
                                                                          aaCACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTAAGGCCCCATTCGCGGGCGCGCCcgc                                                                                                                                     <  2:1098821/90‑1 (MQ=37)
                                                                              cccGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGt                                                                     >  1:1673237/1‑150 (MQ=255)
                                                                                     aaGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGTAAACTCTTCAGTATCCcgc                                                              >  2:1581221/1‑150 (MQ=255)
                                                                                       gATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATggg                                                                                                                              <  1:2165295/84‑1 (MQ=255)
                                                                                       gATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATggg                                                                                                                              >  2:2165295/1‑84 (MQ=255)
                                                                                       gATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAAc                                                                                                        <  1:1849379/106‑1 (MQ=255)
                                                                                       gATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAAc                                                                                                        >  2:1849379/1‑106 (MQ=255)
                                                                                        aTCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTAc                                                                                                               >  1:2210141/1‑98 (MQ=255)
                                                                                        aTCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTAc                                                                                                               <  2:2210141/98‑1 (MQ=255)
                                                                                        aTCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGc                                                           <  2:1673237/150‑1 (MQ=255)
                                                                                          cTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTAAGGCCCCATTCGCGGGCGCGCCCGCGaaga                                                                                                                                <  1:1798644/79‑3 (MQ=21)
                                                                                          cTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTAAGGCCCCATTCGCGGGCGCGCCCGCGaaga                                                                                                                                >  2:1798644/1‑77 (MQ=21)
                                                                                            aTCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGGGCGCCCGCGAATGGGGCCTTACAGGCCTCCCAGTAACTCACGGTAATCCTCCACCGCCGGAAACTCTGCAGTTCCGAGCGGCCCTg                                                       >  2:655753/1‑150 (MQ=255)
                                                                                             tCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGc                                                      >  1:1961975/1‑150 (MQ=255)
                                                                                               cGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGt                                                                                                           <  2:2293859/95‑1 (MQ=255)
                                                                                               cGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGt                                                                                                           >  1:2293859/1‑95 (MQ=255)
                                                                                               cGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCtt                                                    <  2:425599/150‑1 (MQ=255)
                                                                                                ggggCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCttt                                                   <  2:1961975/150‑1 (MQ=255)
                                                                                                           tgcGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCtt                                                                         <  1:16802/117‑1 (MQ=255)
                                                                                                           tgcGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCtt                                                                         >  2:16802/1‑117 (MQ=255)
                                                                                                            gcGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCTGCCCTGCTTTGCTATCCGGCTg                                       >  2:1758839/1‑150 (MQ=255)
                                                                                                            gcGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTg                                       <  2:312580/150‑1 (MQ=255)
                                                                                                             cGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTgg                                      >  1:1537252/1‑150 (MQ=255)
                                                                                                               tATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTcgcg                                                                                                                                                >  1:1298944/1‑42 (MQ=255)
                                                                                                               tATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTcgcg                                                                                                                                                <  2:1298944/42‑1 (MQ=255)
                                                                                                                 ttCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCa                                  >  1:1027614/1‑150 (MQ=255)
                                                                                                                  tCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATACTCGACAGGGGCACTTCCTCCGGTATCGCGCGGCCCTGCTTTGCTATCCGGCTGGCGCAc                                 >  1:2341750/1‑150 (MQ=255)
                                                                                                                   ccTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCAcc                                <  1:1155226/150‑1 (MQ=255)
                                                                                                                     tgactgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTg                                                  <  2:769263/130‑1 (MQ=255)
                                                                                                                     tgactgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTg                                                  >  1:769263/1‑130 (MQ=255)
                                                                                                                        ctgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATggg                                                                                                                              <  2:1691492/51‑1 (MQ=255)
                                                                                                                        ctgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATggg                                                                                                                              >  1:1691492/1‑51 (MQ=255)
                                                                                                                           aTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCaa                        <  2:1628077/150‑1 (MQ=255)
                                                                                                                           aTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCaa                        <  2:2283864/150‑1 (MQ=255)
                                                                                                                            tCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAg                       >  2:1967040/1‑150 (MQ=255)
                                                                                                                                    gTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATACGGCTGGCGCACCGCCAGCAAGTGCCCcac               >  1:2127462/1‑150 (MQ=255)
                                                                                                                                     tCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGccc                  <  2:1449107/146‑1 (MQ=255)
                                                                                                                                     tCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGccc                  >  1:1449107/1‑146 (MQ=255)
                                                                                                                                     tCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCcaca              >  2:2222911/1‑150 (MQ=255)
                                                                                                                                         tACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCaa                        >  2:2172497/1‑136 (MQ=255)
                                                                                                                                         tACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCaa                        <  1:2172497/136‑1 (MQ=255)
                                                                                                                                           cAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAAc        >  1:1483915/1‑150 (MQ=255)
                                                                                                                                            aGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCgg                                                                                                  >  1:648992/1‑59 (MQ=255)
                                                                                                                                            aGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCgg                                                                                                  <  2:648992/59‑1 (MQ=255)
                                                                                                                                              gCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAAccgc     >  1:1249456/1‑150 (MQ=255)
                                                                                                                                                 tctTCGCGGGCGCGCCCGCGAATGGGGCCTTACAAGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCTGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCTGt  >  1:1180294/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCCATTGCCATGGTCGTGTATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGT  >  NC_002947/313706‑314000

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