Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,747,027 T→A 6.3% intergenic (+22/‑104) 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,0270TA 6.3% 274.2 / 6.7 95intergenic (+22/‑104)ydhR/ydhSputative monooxygenase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (52/37);  new base A (4/2);  total (56/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

AAAGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCT  >  NC_000913/1746878‑1747173
                                                                                                                                                     |                                                                                                                                                  
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                                                                                                                                        <  2:1680139/150‑4 (MQ=255)
            ccTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttgat                                                                                                                                        <  1:479078/150‑4 (MQ=255)
             cTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCa                                                                                                                                       <  1:971147/150‑1 (MQ=255)
             cTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCa                                                                                                                                       >  2:1618714/1‑150 (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:751596/1‑150 (MQ=255)
                                 cccGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCa                                                                                                                   >  2:837131/1‑150 (MQ=255)
                                    gACTGAAAAATCTCGTCGGTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAAt                                                                                                                <  2:1101016/150‑1 (MQ=255)
                                        gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                            >  2:1469258/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:1336627/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: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:1455405/1‑150 (MQ=255)
                                                                                                     cAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagtaa                                               >  2:52848/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             >  2:1716927/1‑150 (MQ=255)
                                                                                                                                       tAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGcc             <  1:483013/150‑1 (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   <  2:830593/150‑1 (MQ=255)
                                                                                                                                                 cGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATc   <  1:240957/150‑1 (MQ=255)
                                                                                                                                                  gtttAGCGACATTTTCATCCATGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACATTGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGTATCt  <  2:288138/148‑1 (MQ=255)
                                                                                                                                                  ggTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCt  <  1:519104/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAAGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCT  >  NC_000913/1746878‑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: ‑