Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,747,025 G→C 8.0% intergenic (+20/‑106) 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,0250GC 8.0% 307.8 / 11.3 96intergenic (+20/‑106)ydhR/ydhSputative monooxygenase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (51/36);  new base C (6/2);  total (57/39)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

AAAAGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCT  >  NC_000913/1746877‑1747173
                                                                                                                                                    |                                                                                                                                                    
aaaaGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGt                                                                                                                                                     >  1:1328416/1‑150 (MQ=255)
 aaaGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGtt                                                                                                                                                    <  1:1464066/150‑1 (MQ=255)
   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                                                                                                                                        <  1:479078/150‑4 (MQ=255)
             ccTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttgat                                                                                                                                        <  2:1680139/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                                                                                                                              >  2:148093/1‑150 (MQ=255)
                       aaaaCATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCt                                                                                                                              >  1:940920/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:751596/1‑150 (MQ=255)
                                  cccGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCa                                                                                                                   >  2:837131/1‑150 (MQ=255)
                                     gACTGAAAAATCTCGTCGGTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAAt                                                                                                                <  2:1101016/150‑1 (MQ=255)
                                         gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                            >  2:1336627/1‑150 (MQ=255)
                                         gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                            >  2:1635669/1‑150 (MQ=255)
                                         gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                            >  1:1394784/1‑150 (MQ=255)
                                         gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                            >  2:657970/1‑150 (MQ=255)
                                         gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                            >  2:1469258/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:1635669/150‑1 (MQ=255)
                                                                     gCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAac                                                                                <  1:1618714/150‑1 (MQ=255)
                                                                     gCCAAAGTATTCGATGTCAATGAGCCACTTAGTCAAATCAATCCGGCAAAACTCGCCTGACAGAATTTAATCACGGGCTGTTAGCGCCCTTTTCAGCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTAGCTGATCac                                                                                >  2:1593680/1‑150 (MQ=255)
                                                                      ccAAAGTTTTCGATGTCAATGAGCCTCTTAGTCAATTCAATAAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGAGGTTTGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                               <  2:820333/150‑1 (MQ=255)
                                                                      ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                               <  2:1460104/150‑1 (MQ=255)
                                                                      ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                               >  2:445395/1‑150 (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                                                    >  1:485239/1‑148 (MQ=255)
                                                                                                   aGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgta                                                    <  2:485239/148‑1 (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:1455405/1‑150 (MQ=255)
                                                                                                      cAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagtaa                                               >  2:52848/1‑150 (MQ=255)
                                                                                                           caatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATg                                          <  1:1695645/150‑1 (MQ=255)
                                                                                                            aatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGa                                         <  2:1631924/150‑1 (MQ=255)
                                                                                                            aatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGa                                         >  2:1354723/1‑150 (MQ=255)
                                                                                                               caGGCAAAACTCGCCTGACAGAATTTAATCAAAGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaa                                      >  2:1152236/1‑150 (MQ=255)
                                                                                                                aGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaaa                                     >  1:1684857/1‑150 (MQ=255)
                                                                                                                   cAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAt                                  >  1:1124052/1‑150 (MQ=255)
                                                                                                                   cAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAt                                  >  1:1129876/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             >  2:898039/1‑150 (MQ=255)
                                                                                                                                        tAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGcc             >  2:1716927/1‑150 (MQ=255)
                                                                                                                                        tAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGcc             <  1:483013/150‑1 (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)
                                                                                                                                                   ggTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCt  <  1:519104/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAAAGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCT  >  NC_000913/1746877‑1747173

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: