Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,747,027 T→A 5.1% intergenic (+22/‑104) ydhR → / → ydhS putative monooxygenase YdhR/FAD/NAD(P) binding domain‑containing protein YdhS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,747,0270TA 5.1% 300.5 / 3.1 93intergenic (+22/‑104)ydhR/ydhSputative monooxygenase YdhR/FAD/NAD(P) binding domain‑containing protein YdhS
Reads supporting (aligned to +/- strand):  ref base T (48/40);  new base A (2/3);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACA  >  NC_000913/1746880‑1747176
                                                                                                                                                   |                                                                                                                                                     
aGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAg                                                                                                                                                     <  1:937562/150‑1 (MQ=255)
  cgcCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAgcg                                                                                                                                                   >  2:1243304/1‑150 (MQ=255)
   gcCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAgcgc                                                                                                                                                  >  1:910368/1‑150 (MQ=255)
       ttGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCtt                                                                                                                                              <  1:957203/150‑1 (MQ=255)
         gCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCtttt                                                                                                                                            <  2:896738/150‑1 (MQ=255)
         gCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCtttt                                                                                                                                            <  2:104662/150‑1 (MQ=255)
          ccTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCATAGTTTTCGATGTCAATGAGCCACTTATTCAAATCAATCAGGCATAACTCGCCTGATATAATTTAATCTAGGGCGCTAACCGCCCttgat                                                                                                                                           <  2:937410/150‑4 (MQ=255)
          ccTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTc                                                                                                                                           <  1:1151160/150‑1 (MQ=255)
          ccTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTc                                                                                                                                           >  1:532343/1‑150 (MQ=255)
          ccTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGCTAACCGCCCttgat                                                                                                                                           <  1:786946/150‑4 (MQ=255)
            tATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCAt                                                                                                                                         >  2:579104/1‑150 (MQ=255)
                 ggAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTg                                                                                                                                    >  2:909649/1‑150 (MQ=255)
                 ggAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTg                                                                                                                                    <  2:1239393/150‑1 (MQ=255)
                   aaaaaCATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAAACAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTc                                                                                                                                  >  1:127817/1‑150 (MQ=255)
                       aCATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGaa                                                                                                                              >  1:688773/1‑150 (MQ=255)
                              gCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTtctc                                                                                                                       >  1:696441/1‑150 (MQ=255)
                               cccGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCa                                                                                                                      <  2:1097681/150‑1 (MQ=255)
                                ccGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCATGACAGAATTTAATCAAGGGCGGTTAGCGGCCTTTTCATCGCTGGCTGAAATTTCTGaa                                                                                                                     >  1:609284/1‑150 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTtct                                                                                                                        <  2:1183577/141‑1 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTtct                                                                                                                        >  1:1183577/1‑141 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                               >  1:1277581/1‑150 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                               >  1:246285/1‑150 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                               >  1:9202/1‑150 (MQ=255)
                                      gAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATAAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTAAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTa                                                                                                               >  1:1175706/1‑150 (MQ=255)
                                          aaTCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTaaaaa                                                                                                           <  2:293213/150‑1 (MQ=255)
                                            tctcGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAAtc                                                                                                         <  2:1020401/150‑1 (MQ=255)
                                            tctcGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAAtc                                                                                                         <  1:1182532/150‑1 (MQ=255)
                                              tcGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAAtctc                                                                                                       <  2:532343/150‑1 (MQ=255)
                                              tcGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAAtctc                                                                                                       >  1:776431/1‑150 (MQ=255)
                                                   gTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCa                                                                                                  >  1:1056825/1‑150 (MQ=255)
                                                       aGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACt                                                                                              >  1:685746/1‑150 (MQ=255)
                                                       aGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACt                                                                                              <  2:210022/150‑1 (MQ=255)
                                                         gAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTa                                                                                            <  1:28149/150‑1 (MQ=255)
                                                           aGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATc                                                                                          >  1:100062/1‑150 (MQ=255)
                                                             ttGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTg                                                                                        <  2:1056825/150‑1 (MQ=255)
                                                             ttGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTg                                                                                        >  2:527694/1‑150 (MQ=255)
                                                               gTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGAt                                                                                      <  2:688773/150‑1 (MQ=255)
                                                               gTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGAt                                                                                      <  1:1282667/150‑1 (MQ=255)
                                                                  gCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAac                                                                                   <  1:909649/150‑1 (MQ=255)
                                                                  gCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAac                                                                                   <  1:1263445/150‑1 (MQ=255)
                                                                  gCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCCCAAATTCTAAAAATCTCAACCAAACTTATCTGATAac                                                                                   >  2:1115700/1‑150 (MQ=255)
                                                                   ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                                  >  2:166598/1‑150 (MQ=255)
                                                                   ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                                  >  2:567268/1‑150 (MQ=255)
                                                                   ccAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAaca                                                                                  >  2:1201189/1‑150 (MQ=255)
                                                                    cAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTATTTTCTTTTTTTTTACTGTGTGAAATTTCTTAAATTCTATATTTCTCTATTAAACTTATCTTATAacac                                                                                 <  2:1141105/150‑1 (MQ=255)
                                                                    cAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAacac                                                                                 <  2:1329441/150‑1 (MQ=255)
                                                                     aaaGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACt                                                                                <  2:1215469/150‑1 (MQ=255)
                                                                          tttCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAAtt                                                                           >  2:1377477/1‑150 (MQ=255)
                                                                               aTGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGaaa                                                                      <  2:58108/150‑1 (MQ=255)
                                                                                tGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAg                                                                     <  2:127817/150‑1 (MQ=255)
                                                                                      tGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCg                                                                         >  2:986501/1‑140 (MQ=255)
                                                                                      tGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCg                                                                         <  1:986501/140‑1 (MQ=255)
                                                                                      tGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCg                                                               <  2:696441/150‑1 (MQ=255)
                                                                                           cACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAg                                                          >  2:61362/1‑150 (MQ=255)
                                                                                               tAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaa                                                      >  2:824014/1‑150 (MQ=255)
                                                                                               tAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTCTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTCATAACACTAAATTCGAAAGAATGCGTACATgtaa                                                      >  2:205713/1‑150 (MQ=255)
                                                                                                 ctCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGGTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagt                                                    <  2:1340455/149‑1 (MQ=255)
                                                                                                   cAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagtaa                                                  >  2:821412/1‑150 (MQ=255)
                                                                                                      atcaatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagt                                                    <  2:719903/145‑1 (MQ=255)
                                                                                                      atcaatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGgtaagt                                                    >  1:719903/1‑145 (MQ=255)
                                                                                                         aatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGa                                            <  1:64516/150‑1 (MQ=255)
                                                                                                         aatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGa                                            >  1:338925/1‑150 (MQ=255)
                                                                                                         aatcaGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCACAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGa                                            >  2:1155581/1‑150 (MQ=255)
                                                                                                            caGGCAAAACTCGCCTGACAGAATTTAATCAATGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaa                                         <  1:826407/150‑1 (MQ=255)
                                                                                                            caGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaa                                         <  2:9202/150‑1 (MQ=255)
                                                                                                             aGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTc                                                                          >  1:95386/1‑116 (MQ=255)
                                                                                                             aGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTc                                                                          <  2:95386/116‑1 (MQ=255)
                                                                                                             aGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaaa                                        >  1:1136452/1‑150 (MQ=255)
                                                                                                              ggCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGGATGCGTACAGGTAAGTAACAATGaaaaaa                                       >  1:421931/1‑150 (MQ=255)
                                                                                                              ggCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaaaa                                       >  2:1292004/1‑150 (MQ=255)
                                                                                                              ggCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaaaa                                       >  2:897582/1‑150 (MQ=255)
                                                                                                              ggCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaaaa                                       <  2:1277581/150‑1 (MQ=255)
                                                                                                              ggCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCACTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGaaaaaa                                       >  1:1209657/1‑150 (MQ=255)
                                                                                                                cAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAt                                     >  1:881499/1‑150 (MQ=255)
                                                                                                                 aaaaCTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAtt                                    <  1:174197/150‑1 (MQ=255)
                                                                                                                  aaaCTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAttt                                   >  1:422923/1‑149 (MQ=255)
                                                                                                                   aaCTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTCTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATATGATAACACTAAATTCGAAAGAATGCGTACACGTAAGTAACAATGAAAAAAATTGc                                  >  2:921908/1‑150 (MQ=255)
                                                                                                                    aCTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCt                                 <  1:166598/150‑1 (MQ=255)
                                                                                                                    aCTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCt                                 <  1:1201189/150‑1 (MQ=255)
                                                                                                                           tGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTgg                          >  2:598560/1‑150 (MQ=255)
                                                                                                                            gACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTggg                         >  2:1386165/1‑150 (MQ=255)
                                                                                                                             aCAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGt                        >  2:321116/1‑150 (MQ=255)
                                                                                                                                gAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGcc                     >  2:765619/1‑150 (MQ=255)
                                                                                                                                     tAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGcc                <  2:910368/150‑1 (MQ=255)
                                                                                                                                       aTCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTa              >  1:129930/1‑150 (MQ=255)
                                                                                                                                        tCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCACAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTAGAAAGAATGCGTACACGTAAGTAACAATGAAAAAAATTGATATTGTGGGTGCCGGGCCTAc             >  1:1408262/1‑150 (MQ=255)
                                                                                                                                        aaaaGGGCGCTAACCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTAc             >  1:598982/3‑150 (MQ=255)
                                                                                                                                          aaGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAATTTGCTATTGTGGGTGCCGGGCCTACgg           >  1:816144/1‑150 (MQ=255)
                                                                                                                                           aGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACggg          <  1:1243304/150‑1 (MQ=255)
                                                                                                                                              gCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGAt       >  2:908343/1‑150 (MQ=255)
                                                                                                                                               cGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATc      <  1:824014/150‑1 (MQ=255)
                                                                                                                                               cGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATc      <  1:541735/150‑1 (MQ=255)
                                                                                                                                                ggTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCt     <  1:908343/150‑1 (MQ=255)
                                                                                                                                                   tAGCGCCCTTTTCATCCCTGTCTGAAATTTCCCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTaca  <  1:1115700/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCGCCCTTGCCTATCTGGAAAAACATACTGCCCGACTGAAAAATCTCGGCGTTGAGGAAGTTGTCGCCAAAGTTTTCGATGTCAATGAGCCACTTAGTCAAATCAATCAGGCAAAACTCGCCTGACAGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACA  >  NC_000913/1746880‑1747176

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