Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,747,029 G→C 6.0% intergenic (+24/‑102) ydhR → / → ydhS putative monooxygenase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,747,0290GC 6.0% 298.7 / 7.5 96intergenic (+24/‑102)ydhR/ydhSputative monooxygenase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (52/36);  new base C (4/2);  total (56/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

AGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACC  >  NC_000913/1746880‑1747178
                                                                                                                                                     |                                                                                                                                                     
aGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAg                                                                                                                                                       <  2:1413794/150‑1 (MQ=255)
    cccTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGcc                                                                                                                                                   <  1:95549/150‑1 (MQ=255)
    cccTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGcc                                                                                                                                                   >  2:519104/1‑150 (MQ=255)
        tGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCttt                                                                                                                                               >  1:472868/1‑150 (MQ=255)
          ccTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttgat                                                                                                                                             <  2:1680139/150‑4 (MQ=255)
          ccTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttgat                                                                                                                                             <  1:479078/150‑4 (MQ=255)
           cTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCa                                                                                                                                            >  2:1618714/1‑150 (MQ=255)
           cTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCa                                                                                                                                            <  1:971147/150‑1 (MQ=255)
               cTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATccc                                                                                                                                        <  2:1466249/150‑1 (MQ=255)
                 ggAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTg                                                                                                                                      >  2:240957/1‑150 (MQ=255)
                    aaaaCATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCt                                                                                                                                   >  1:940920/1‑150 (MQ=255)
                    aaaaCATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCt                                                                                                                                   >  2:148093/1‑150 (MQ=255)
                    aaaaCATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGAGCCCTTTTCATCCCTGtat                                                                                                                                   >  1:130347/1‑148 (MQ=255)
                       aCATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGAAAGATTTTAATCAATGGCGGTTAGCGCCCTTTTCCTCCGTGTCTGaa                                                                                                                                >  1:859631/1‑150 (MQ=255)
                             tGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTAAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTtct                                                                                                                          <  2:260601/150‑1 (MQ=255)
                              gCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTTAGCGCCCTTTTCATCCCTGTCTGAAATTtctc                                                                                                                         >  1:1610423/1‑150 (MQ=255)
                               cccGACTGAAAAATCTCTGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCa                                                                                                                        <  2:531709/150‑1 (MQ=255)
                               cccGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCa                                                                                                                        >  2:837131/1‑150 (MQ=255)
                               cccGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCa                                                                                                                        >  2:751596/1‑150 (MQ=255)
                                  gACTGAAAAATCTCGTCGGTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAAt                                                                                                                     <  2:1101016/150‑1 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                                 >  2:1336627/1‑150 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                                 >  2:1469258/1‑150 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                                 >  2:1635669/1‑150 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                                 >  2:657970/1‑150 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                                 >  1:1394784/1‑150 (MQ=255)
                                           aTCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCGGGACAGAATTTAATCAAGGGCGGTTAGGGCCCTTTTCATCCCTGACTGATATTTCTCTAATGCTaataat                                                                                                            >  2:1604188/1‑150 (MQ=255)
                                            tctcGGCGTTGAGGAAGTTGTCGGCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAAtc                                                                                                           <  1:20052/150‑1 (MQ=255)
                                             ctcGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAAtct                                                                                                          >  2:228646/1‑150 (MQ=255)
                                             ctcGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAAtct                                                                                                          <  1:972155/150‑1 (MQ=255)
                                                    ttGAGGAAGTTGGCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCaa                                                                                                   <  2:842273/150‑1 (MQ=255)
                                                     tGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTTACAGTATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCaaa                                                                                                  <  1:1533428/150‑1 (MQ=255)
                                                     tGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCaaa                                                                                                  >  2:224370/1‑150 (MQ=255)
                                                      gAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAAACAAACAGGCAAAACGCGCATCACAGATTTTAAACAAGGGCGGTTAGCGCCCTTTAAATCCAGGTATGAAATTTATCAAATACTAAAAAAATCAACCaaaa                                                                                                 >  2:483013/1‑149 (MQ=255)
                                                       aGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACt                                                                                                <  2:1377317/150‑1 (MQ=255)
                                                        ggAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACtt                                                                                               >  2:1695645/1‑150 (MQ=255)
                                                              tGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGa                                                                                         <  1:751596/150‑1 (MQ=255)
                                                               gTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGAt                                                                                        <  1:228646/150‑1 (MQ=255)
                                                                tCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATa                                                                                       <  2:1394784/150‑1 (MQ=255)
                                                                  gCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAac                                                                                     <  1:1618714/150‑1 (MQ=255)
                                                                  gCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAac                                                                                     <  1:1635669/150‑1 (MQ=255)
                                                                  gCCAAAGTATTCGATGTCAATGAGCCACTTAGTCAAATCAATCCGGCAAAACTCGCCTGACAGAATTTAATCACGGGCTGTTAGCGCCCTTTTCAGCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTAGCTGATCac                                                                                     >  2:1593680/1‑150 (MQ=255)
                                                                   ccAAAGTTTTCGATGTCAATGAGCCTCTTAGTCAATTCAATAAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGAGGTTTGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                                    <  2:820333/150‑1 (MQ=255)
                                                                   ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                                    >  2:445395/1‑150 (MQ=255)
                                                                   ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                                    <  2:1460104/150‑1 (MQ=255)
                                                                    cAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAacac                                                                                   <  2:1549300/150‑1 (MQ=255)
                                                                     aaaGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACt                                                                                  >  1:115631/1‑150 (MQ=255)
                                                                       aGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTaa                                                                                >  2:1514821/1‑150 (MQ=255)
                                                                             cGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGa                                                                          >  2:631146/1‑150 (MQ=255)
                                                                               aTGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCATTTTCATCCCTGTCTGAAATTTCTAAAATTCTAAAAATCTCAACCTAAATTATCTGATAACACTAAATTCTaaa                                                                        >  1:288138/1‑150 (MQ=255)
                                                                                   cAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAAt                                                                    <  2:940920/150‑1 (MQ=255)
                                                                                   cAAAGAGCCACTTAGACAAAACAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTTAGCGCCCTTATCATCCCTGTCTGAAATTTCTCAAATTCTAAAACTCTCAACCAAACTTATCTGATAACAATAAATTCGAAAGaaa                                                                    >  1:1708360/1‑149 (MQ=255)
                                                                                     aTGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGc                                                                  >  1:229244/1‑150 (MQ=255)
                                                                                      tGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCg                                                                 <  2:115631/150‑1 (MQ=255)
                                                                                              ttAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgta                                                         <  1:148093/150‑1 (MQ=255)
                                                                                               tAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaa                                                        >  1:1114858/1‑150 (MQ=255)
                                                                                                aGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaag                                                       <  2:229244/150‑1 (MQ=255)
                                                                                                aGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgta                                                         <  2:485239/148‑1 (MQ=255)
                                                                                                aGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgta                                                         >  1:485239/1‑148 (MQ=255)
                                                                                                 gTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagt                                                      >  1:830593/1‑150 (MQ=255)
                                                                                                 ctCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAAATCTAAAAATCTCAACCAAACTAATCTGATAACACTAAAATCGAAAGAATGCGTACAGgtaagt                                                      >  2:471629/2‑150 (MQ=255)
                                                                                                  tCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCATGTCTGAAATTTATCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACTCTAAATTCGAAAGAATGCGTTCAGgtaagtt                                                     >  1:1087449/1‑149 (MQ=255)
                                                                                                   cAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGCATGCGTACAGgtaagtaa                                                    >  1:1394056/1‑150 (MQ=255)
                                                                                                   cAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagtaa                                                    >  2:52848/1‑150 (MQ=255)
                                                                                                   cAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagtaa                                                    >  2:1455405/1‑150 (MQ=255)
                                                                                                        caatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATg                                               <  1:1695645/150‑1 (MQ=255)
                                                                                                         aatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGa                                              >  2:1354723/1‑150 (MQ=255)
                                                                                                         aatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGa                                              <  2:1631924/150‑1 (MQ=255)
                                                                                                            caGGCAAAACTCGCCTGACAGAATTTAATCAAAGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaa                                           >  2:1152236/1‑150 (MQ=255)
                                                                                                             aGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaaa                                          >  1:1684857/1‑150 (MQ=255)
                                                                                                                cAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAt                                       >  1:1129876/1‑150 (MQ=255)
                                                                                                                cAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAt                                       >  1:1124052/1‑150 (MQ=255)
                                                                                                                   aaCTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGc                                    <  2:1114858/150‑1 (MQ=255)
                                                                                                                    aCTCGCCTGACAGAATTTAATCAAGTGCGTTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCt                                   <  2:1610423/150‑1 (MQ=255)
                                                                                                                            gACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTggg                           >  1:451007/1‑150 (MQ=255)
                                                                                                                               aGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGc                        <  1:1336627/150‑1 (MQ=255)
                                                                                                                                gAATTTAATCACGGGCGGTTATCGCCCTTTTCATCCCTGTCTGAAATTTCTCTAATTTTAAATATCTCATCCAATCTTATCTGATAACACTAAATTAGAAAGAATGCGTACAGGTAAGTAACAATGAAATAAATTGCTATTGTGTGTGcc                       <  2:1708360/150‑1 (MQ=255)
                                                                                                                                gAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGcc                       >  2:449545/1‑150 (MQ=255)
                                                                                                                                gAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGATATTGTGGGTGcc                       >  2:1041099/1‑150 (MQ=255)
                                                                                                                                   tttAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCggg                    >  2:1694755/1‑150 (MQ=255)
                                                                                                                                    gatgaaaaGGGCGCTAACCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAACCTTATCTGATAACACTAAATTCGTAAGAATGCGTACAGGTAAGTTACAATGAAAAAAATTGCTATTGTGGGTGCCGTGc                   >  2:1078141/7‑150 (MQ=255)
                                                                                                                                    gatgaaaaGGGCGCTAACCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGATATTGTGTGTTCCTGGc                   >  2:326873/7‑150 (MQ=255)
                                                                                                                                     tAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGcc                  <  1:483013/150‑1 (MQ=255)
                                                                                                                                     tAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGcc                  >  2:1716927/1‑150 (MQ=255)
                                                                                                                                     tAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGcc                  >  2:898039/1‑150 (MQ=255)
                                                                                                                                       aTCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTa                >  2:597051/1‑150 (MQ=255)
                                                                                                                                        aaaaGGGCGCTAACCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATTAAAAAAATTGCTATTGTGGGTCAAGGGCCTaa               >  1:145204/3‑149 (MQ=255)
                                                                                                                                        aaaaGGGCGCTAACCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTAc               >  1:104808/3‑150 (MQ=255)
                                                                                                                                          aaGGGCGGTTATCGCCCTTTTCCTCCCTGTCTGAAATTACTCAAATTCTAAAAATCTCAACCAACCTTATCTGATTAGACTAAATTCGTATTAATGCGGACAGGTAAGTAACAAAGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACgg             <  1:1593680/150‑1 (MQ=255)
                                                                                                                                          aaGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACgg             >  1:135699/1‑150 (MQ=255)
                                                                                                                                          aaGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACgg             >  2:1324898/1‑150 (MQ=255)
                                                                                                                                               cGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAATATCTCAACCAAACTTTTCTGATAACACTATATTCGAATGAATGCGTTCTGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATc        <  1:471629/150‑1 (MQ=255)
                                                                                                                                               cGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATc        >  2:709197/1‑150 (MQ=255)
                                                                                                                                               cGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATc        <  1:240957/150‑1 (MQ=255)
                                                                                                                                               cGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATc        <  2:830593/150‑1 (MQ=255)
                                                                                                                                                gtttAGCGACATTTTCATCCATGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACATTGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGTATCt       <  2:288138/148‑1 (MQ=255)
                                                                                                                                                ggTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCt       <  1:519104/150‑1 (MQ=255)
                                                                                                                                                    aGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTacac   <  1:1324898/150‑1 (MQ=255)
                                                                                                                                                     gcgcCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAATATTTTTTTTTTTTTGTTGCGGGGACTACGGTGATCTACAcc  <  2:1394056/150‑1 (MQ=255)
                                                                                                                                                     gcgcCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACAcc  <  1:837131/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACC  >  NC_000913/1746880‑1747178

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