Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,127,799 C→A 9.6% intergenic (+70/‑40) yceM → / → murJ putative oxidoreductase YceM/putative lipid II flippase MurJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,127,7990CA 9.6% 336.5 / 16.5 104intergenic (+70/‑40)yceM/murJputative oxidoreductase YceM/putative lipid II flippase MurJ
Reads supporting (aligned to +/- strand):  ref base C (47/47);  new base A (6/4);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 107 total aligned reads displayed.

TCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAGGGATG  >  NC_000913/1127651‑1127946
                                                                                                                                                    |                                                                                                                                                   
tCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCgcctgcct                                                                                                                                                    <  2:1582084/150‑1 (MQ=255)
 ccGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCtt                                                                                                                                                   >  2:787095/1‑150 (MQ=255)
  cGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGGCTCTTTCAGCGCCTGCCTTg                                                                                                                                                  >  2:1007675/1‑150 (MQ=255)
  cGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTg                                                                                                                                                  >  2:343367/1‑150 (MQ=255)
     aaaCCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAg                                                                                                                                               <  1:860615/150‑1 (MQ=255)
          gccgGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGtt                                                                                                                                          <  1:301245/150‑1 (MQ=255)
           ccgGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGttt                                                                                                                                         >  1:594194/1‑150 (MQ=255)
           ccgGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCAAGGCAGGCgctg                                                                                                                                         <  2:1194643/150‑4 (MQ=255)
             ggCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCAAGGCAGGCgctgaa                                                                                                                                       <  1:991599/150‑6 (MQ=255)
             ggCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAACGCCTGCAAGGCAGGCgctgaa                                                                                                                                       <  2:1272325/150‑6 (MQ=255)
              gCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGc                                                                                                                                      <  1:84060/150‑1 (MQ=255)
               cGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCAAGGCAGGCgctgaaag                                                                                                                                     <  1:491067/150‑8 (MQ=255)
                gAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGccc                                                                                                                                    >  1:363278/1‑150 (MQ=255)
                  aCAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGt                                                                                                                                  >  1:1279557/1‑150 (MQ=255)
                         gTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTg                                                                                                                           <  2:134362/150‑1 (MQ=255)
                            cTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGAt                                                                                                                        >  1:1334058/1‑150 (MQ=255)
                             tGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGAtt                                                                                                                       <  2:153722/150‑1 (MQ=255)
                             tGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGAtt                                                                                                                       <  1:345008/150‑1 (MQ=255)
                             tGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGAtt                                                                                                                       >  1:1415202/1‑150 (MQ=255)
                             tGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGAtt                                                                                                                       <  1:925147/150‑1 (MQ=255)
                              ggCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTa                                                                                                                      <  1:1160530/150‑1 (MQ=255)
                              ggCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTa                                                                                                                      <  2:563004/150‑1 (MQ=255)
                               gcggAACGTATCGTTGACAAGATCTGGCGGGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAg                                                                                                                     <  2:1627281/150‑1 (MQ=255)
                                 gcAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGaa                                                                                                                   <  2:132297/150‑1 (MQ=255)
                                  cAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAac                                                                                                                  >  1:1058196/1‑150 (MQ=255)
                                     cGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTTGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACAcc                                                                                                               <  1:1558251/150‑1 (MQ=255)
                                     cGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACAcc                                                                                                               <  2:565571/150‑1 (MQ=255)
                                         tCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATg                                                                                                           >  2:889303/1‑150 (MQ=255)
                                         tCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATg                                                                                                           <  2:813259/150‑1 (MQ=255)
                                          cGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGa                                                                                                          >  2:32919/1‑150 (MQ=255)
                                           gTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGaa                                                                                                         >  2:1214136/1‑150 (MQ=255)
                                             ggACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAAtt                                                                                                       <  2:1510949/149‑1 (MQ=255)
                                              gACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAAttt                                                                                                      >  2:1610544/1‑150 (MQ=255)
                                                cAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATttat                                                                                                    >  1:1176930/1‑150 (MQ=255)
                                                  aGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGCATTAGAACACCGATGAATttatta                                                                                                  >  1:430000/1‑150 (MQ=255)
                                                   gATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTaa                                                                                                 <  1:881180/150‑1 (MQ=255)
                                                   gATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTaa                                                                                                 <  1:1280062/150‑1 (MQ=255)
                                                    aTCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTaaa                                                                                                >  1:124990/1‑150 (MQ=255)
                                                     tCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTaaaa                                                                                               >  1:423475/1‑150 (MQ=255)
                                                      cTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAAt                                                                                              >  1:1579051/1‑150 (MQ=255)
                                                          cgcgATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCt                                                                                          >  1:1441953/1‑150 (MQ=255)
                                                            cgatgcgatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTgg                                                                                        >  2:1272499/1‑150 (MQ=255)
                                                                 cgatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGgccgcc                                                                                   >  2:366876/1‑150 (MQ=255)
                                                                  gatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGccgccg                                                                                  <  1:830503/150‑1 (MQ=255)
                                                                    tgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTc                                                                                <  2:1391604/150‑1 (MQ=255)
                                                                      aGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAg                                                                              >  1:20084/1‑150 (MQ=255)
                                                                         taaTAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTc                                                                           <  2:846786/149‑1 (MQ=255)
                                                                         gAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTc                                                                           >  2:958054/1‑150 (MQ=255)
                                                                         gAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTc                                                                           <  1:343367/150‑1 (MQ=255)
                                                                         gAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGGCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTc                                                                           >  1:151563/1‑150 (MQ=255)
                                                                           aTAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGa                                                                         >  2:1332364/1‑150 (MQ=255)
                                                                           aTAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGa                                                                         >  1:1529610/1‑150 (MQ=255)
                                                                             aaaCCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATg                                                                       <  1:127402/150‑1 (MQ=255)
                                                                               aCCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGAc                                                                     <  2:1071628/150‑1 (MQ=255)
                                                                                      aaCATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCa                                                                               >  2:363513/1‑133 (MQ=255)
                                                                                      aaCATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCa                                                                               <  1:363513/133‑1 (MQ=255)
                                                                                         aTCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTcg                                                           <  1:787095/150‑1 (MQ=255)
                                                                                          tCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTcgc                                                          >  2:794360/1‑150 (MQ=255)
                                                                                          tCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTcgc                                                          <  1:162414/150‑1 (MQ=255)
                                                                                           cTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTcgcg                                                         >  2:1497051/1‑150 (MQ=255)
                                                                                           cTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTcgcg                                                         >  2:199033/1‑150 (MQ=255)
                                                                                             ggCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGGCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCgtg                                                       <  2:151563/150‑1 (MQ=255)
                                                                                                ggTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGc                                                    >  1:1630861/1‑150 (MQ=255)
                                                                                                  tAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCtt                                                  >  1:1184258/1‑150 (MQ=255)
                                                                                                      aaTTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCt                                              >  2:730290/1‑150 (MQ=255)
                                                                                                       aTTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGTCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCtt                                             >  1:1625435/1‑150 (MQ=255)
                                                                                                       aTTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCtt                                             >  2:1542263/1‑150 (MQ=255)
                                                                                                        ttCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTc                                            <  1:1272499/150‑1 (MQ=255)
                                                                                                           aCCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCa                                         <  2:1058196/150‑1 (MQ=255)
                                                                                                           aCCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCa                                         <  1:366876/150‑1 (MQ=255)
                                                                                                             cATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACg                                       <  1:337581/150‑1 (MQ=255)
                                                                                                             cATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACg                                       <  1:331669/150‑1 (MQ=255)
                                                                                                             cATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACg                                       <  2:1238494/150‑1 (MQ=255)
                                                                                                               tAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACgag                                     >  2:1338239/1‑150 (MQ=255)
                                                                                                               tAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACgag                                     <  2:1409243/150‑1 (MQ=255)
                                                                                                                aaTCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACgaga                                    <  2:1579051/150‑1 (MQ=255)
                                                                                                                  tCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGCCCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACg                                  >  1:1488625/1‑150 (MQ=255)
                                                                                                                  tCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACg                                  >  1:205286/1‑150 (MQ=255)
                                                                                                                   ccAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGc                                 <  2:124990/150‑1 (MQ=255)
                                                                                                                      ggTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAAt                              >  2:1434938/1‑150 (MQ=255)
                                                                                                                       gTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAAtt                             >  1:972350/1‑150 (MQ=255)
                                                                                                                         aGACTATTCGGCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGt                           <  1:1007675/150‑1 (MQ=255)
                                                                                                                           aCTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCg                         >  1:417122/1‑150 (MQ=255)
                                                                                                                             tATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGcc                       <  1:1338239/150‑1 (MQ=255)
                                                                                                                             tATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGcc                       <  1:1471291/150‑1 (MQ=255)
                                                                                                                             tATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGGCTGGAGTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGcc                       <  1:893803/150‑1 (MQ=255)
                                                                                                                                 cGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGaa                   >  1:739549/1‑150 (MQ=255)
                                                                                                                                    ctctTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATCAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCt                >  2:1552361/1‑150 (MQ=255)
                                                                                                                                      ggcaaaacgcCTGCAAGGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCttt              >  2:1507174/8‑150 (MQ=255)
                                                                                                                                       tttCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTg             >  1:1060635/1‑150 (MQ=255)
                                                                                                                                       tttCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTg             <  2:1279557/150‑1 (MQ=255)
                                                                                                                                       gcaaaacgcCTGCAAGGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAGCACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTg             >  1:524680/7‑150 (MQ=255)
                                                                                                                                        ttCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTgg            <  2:1626008/150‑1 (MQ=255)
                                                                                                                                        ttCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTgg            <  2:1529610/150‑1 (MQ=255)
                                                                                                                                        caaaacgcCTGCAAGGCAGGCGTTTTGCCCGTGGGTCTGGATTAGACCACCGATGACTTTATTAAAATCGCTGGCCGCCGTCAGCTGGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCGCGCGACGCAATTGTCGCCAGAATCTTTgg            >  2:1114298/6‑150 (MQ=255)
                                                                                                                                        caaaacgcCTGCAAGGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTgg            >  2:1038338/6‑150 (MQ=255)
                                                                                                                                        caaaacgcCTGCAAGGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTgg            >  2:85916/6‑150 (MQ=255)
                                                                                                                                        caaaacgcCTGCAAGGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTga            >  1:1406666/6‑149 (MQ=255)
                                                                                                                                                ctgcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAGGGa    >  1:453958/1‑150 (MQ=255)
                                                                                                                                                  gcctTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAGGGATg  >  1:32096/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGCCCGTGGGTCTGGATTAGAACACCGATGAATTTATTAAAATCGCTGGCCGCCGTCAGCTCGATGACCATGTTTTCGCGTGTGCTTGGCTTCGCACGAGACGCAATTGTCGCCAGAATCTTTGGCGCAGGGATG  >  NC_000913/1127651‑1127946

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