Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,157,629 G→C 10.3% intergenic (+55/‑240) ycfH → / → ptsG putative DNase/fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,157,6290GC10.3% 321.6 / 24.3 107intergenic (+55/‑240)ycfH/ptsGputative DNase/fused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base G (44/51);  new base C (5/6);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

Only 100 of 113 total aligned reads displayed.

GGTGTTGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAGGGC  >  NC_000913/1157482‑1157778
                                                                                                                                                   |                                                                                                                                                     
ggTGTTGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGaa                                                                                                                                                     >  2:818383/1‑150 (MQ=255)
 gtgtTGCCGTTGAAGAACTGGCGCAGGTTACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCGATCCGTGGAATGAGTTTTTTTAAAGCTCGTAAGTAATGGCTAAAACGAGTAAAGTTCACCGCCGaaa                                                                                                                                                    >  1:140205/1‑150 (MQ=255)
   gtTGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAAt                                                                                                                                                  <  1:54457/150‑1 (MQ=255)
        cGTTGAAGAACTGGCGCAGGAAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCATCCCGCCTTCAATCCATCCGTTGACTGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGc                                                                                                                                             <  2:1111576/150‑1 (MQ=255)
         gTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCg                                                                                                                                            >  2:521586/1‑150 (MQ=255)
          ttGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCgg                                                                                                                                           <  1:594778/150‑1 (MQ=255)
           tGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGt                                                                                                                                          <  1:1076387/150‑1 (MQ=255)
           tGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGt                                                                                                                                          <  1:1177581/150‑1 (MQ=255)
           tGAAGAACTGGCGCAGGTAACCACCGATAACTTAGCAAGTCTGTTTCACATCGACGCATCCCGCATTCAATCCATCCGGTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAAACAGTAAAGTTCACCGACGATAATTGGGCTGt                                                                                                                                          >  1:281869/1‑150 (MQ=255)
            gaagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTg                                                                                                                                         >  1:985948/1‑150 (MQ=255)
            gaagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTg                                                                                                                                         <  2:828678/150‑1 (MQ=255)
            gaagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTg                                                                                                                                         <  2:568888/150‑1 (MQ=255)
            gaagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACAGCGGAAAATTGGGCGGTg                                                                                                                                         >  1:172806/1‑150 (MQ=255)
              agaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGaa                                                                                                                                       >  1:330978/1‑150 (MQ=255)
               gaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAAt                                                                                                                                      >  1:825057/1‑150 (MQ=255)
               gaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTACCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGaac                                                                                                                                      <  1:1490045/150‑2 (MQ=255)
                aaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATa                                                                                                                                     <  2:607177/150‑1 (MQ=255)
                aaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGaact                                                                                                                                     <  1:1541875/150‑3 (MQ=255)
                    ggCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCa                                                                                                                                 >  1:444108/1‑150 (MQ=255)
                       gcAGGTAACCACCGATAACTTCGCCCGGCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGt                                                                                                                              <  2:1556922/150‑1 (MQ=255)
                        cAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGtt                                                                                                                             <  2:705023/150‑1 (MQ=255)
                        cAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGtt                                                                                                                             >  2:732146/1‑150 (MQ=255)
                              accaccGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAat                                                                                                                       >  2:1350083/1‑150 (MQ=255)
                              accaccGAAAACTTCGCCCGTCTGTTTCACTTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAat                                                                                                                       <  1:731551/150‑1 (MQ=255)
                               ccaccGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAata                                                                                                                      >  1:916697/1‑150 (MQ=255)
                                 accGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATAtt                                                                                                                    <  1:732146/150‑1 (MQ=255)
                                  ccGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtg                                                                                                                   <  2:863269/150‑1 (MQ=255)
                                    gATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCTGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtgtg                                                                                                                 <  1:1413315/150‑1 (MQ=255)
                                    gATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtgtg                                                                                                                 <  2:1655579/150‑1 (MQ=255)
                                      tAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGAc                                                                                                               <  1:1436475/150‑1 (MQ=255)
                                      tAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGAc                                                                                                               <  2:934893/150‑1 (MQ=255)
                                      tAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGAc                                                                                                               <  1:1090152/150‑1 (MQ=255)
                                           tCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATg                                                                                                          <  1:1514173/150‑1 (MQ=255)
                                                cGTCTGTTTCACATCGACGCTACCCGCCATCAATCAATCCGTTGAATGATTTTTTTTAAAGCTCTTAATTAATTGCTAAAACGAGTAAAGTTAAACGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTg                                                                                                     <  1:865685/150‑1 (MQ=255)
                                                 gTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGt                                                                                                    <  2:852374/150‑1 (MQ=255)
                                                       ttCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAt                                                                                              >  1:354933/1‑150 (MQ=255)
                                                        tCACATCGCCGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtg                                                                                             <  1:1410252/150‑1 (MQ=255)
                                                        tCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtg                                                                                             >  1:887795/1‑150 (MQ=255)
                                                          acaTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtgtg                                                                                           >  1:1069407/1‑150 (MQ=255)
                                                          acaTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtgtg                                                                                           <  2:593992/150‑1 (MQ=255)
                                                           caTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAATGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGc                                                                                          <  1:425340/150‑1 (MQ=255)
                                                           caTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGc                                                                                          <  1:632317/150‑1 (MQ=255)
                                                            ccccccGCTCCCCGCCCTCAAGCCATCCCGAGAAAGACTTTTTATAAAGCTCGTACTTCCTCCCTCAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAACTATTGTGACATATGTTTTGTCAAAATGTGCa                                                                                         <  2:1133417/145‑1 (MQ=255)
                                                             tCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCaa                                                                                        >  2:937373/1‑150 (MQ=255)
                                                               gACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACt                                                                                      >  1:541024/1‑150 (MQ=255)
                                                                aCGCTTCCCGCCTTCTATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAACTTAACCGCCGATATTTGGGCGGTGAATAACCACGTTTGAAATATTGAGAAATATGTTTTGTGAAAATGTTCAACtt                                                                                     >  1:278033/1‑150 (MQ=255)
                                                                 cGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACtta                                                                                    >  2:1636111/1‑149 (MQ=255)
                                                                 cGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTtc                                                                                    <  1:417604/150‑1 (MQ=255)
                                                                  gCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTtct                                                                                   <  1:1525335/150‑1 (MQ=255)
                                                                      cccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCaa                                                                               >  1:1695498/1‑150 (MQ=255)
                                                                      cccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCaa                                                                               >  2:549971/1‑150 (MQ=255)
                                                                      cccGCCTTCAATCCATACGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCTCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCaa                                                                               >  2:658448/1‑150 (MQ=255)
                                                                           cTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATc                                                                          <  2:1218525/150‑1 (MQ=255)
                                                                             tCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTg                                                                        >  2:1697361/1‑150 (MQ=255)
                                                                              cAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGa                                                                       >  2:1244725/1‑150 (MQ=255)
                                                                                atccatccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAg                                                                     <  1:1376558/150‑1 (MQ=255)
                                                                                  ccatccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGtt                                                                   <  1:529868/150‑1 (MQ=255)
                                                                                      ccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTATCCAATGATCTGAAGTTGaaa                                                               >  1:138654/1‑150 (MQ=255)
                                                                                        gTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACg                                                             <  1:767355/150‑1 (MQ=255)
                                                                                           gAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTAGTCAAAATGTGCAACATCTCCAATGATCTGAAGTTGAAACTTGa                                                          >  2:401845/1‑150 (MQ=255)
                                                                                                aGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGcc                                                     <  2:541024/150‑1 (MQ=255)
                                                                                                aGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTCGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGcc                                                     >  1:690005/1‑150 (MQ=255)
                                                                                                 gTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCg                                                    >  2:1411373/1‑150 (MQ=255)
                                                                                                 gTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCg                                                    <  2:880401/150‑1 (MQ=255)
                                                                                                  tttttttAAATCTCGTAATTAATGGATAAAACGAGTAAAGTTCACCGCCGAAAATTGGGAGGTTAATAACCTCGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGt                                                   <  2:813273/150‑1 (MQ=255)
                                                                                                   ttttttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTc                                                  <  1:937489/150‑1 (MQ=255)
                                                                                                     ttttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTcaa                                                >  1:1065714/1‑150 (MQ=255)
                                                                                                            gCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAAtt                                         <  1:1350083/150‑1 (MQ=255)
                                                                                                             cTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAAAGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTg                                        >  2:1239936/1‑150 (MQ=255)
                                                                                                              tCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTgg                                       >  1:1421140/1‑150 (MQ=255)
                                                                                                              tCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTgg                                       <  2:488656/150‑1 (MQ=255)
                                                                                                                gTAATTAATGGCTAAAACGAGTAAAGTTCACAGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCa                                     <  2:1288716/150‑1 (MQ=255)
                                                                                                                  aattaatGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACt                                   <  2:126377/150‑1 (MQ=255)
                                                                                                                       atGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATTGCCGTCAAACTAATTGGCACTGAatt                              >  1:161043/1‑150 (MQ=255)
                                                                                                                          gCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAattatt                           >  1:1733292/1‑150 (MQ=255)
                                                                                                                           cTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTAttt                          <  2:1069407/150‑1 (MQ=255)
                                                                                                                           cTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTAttt                          >  1:1661870/1‑150 (MQ=255)
                                                                                                                           cTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTAttt                          >  2:1144241/1‑150 (MQ=255)
                                                                                                                             aaaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTa                        >  1:1429297/1‑150 (MQ=255)
                                                                                                                             aaaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTa                        >  1:766467/1‑150 (MQ=255)
                                                                                                                              aaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTAc                       <  1:960934/150‑1 (MQ=255)
                                                                                                                                aCGAGTACAGTTCACCGCCTAAAATAGGGCGGGGAATAACCCCGTTTGAAAAATTGTGTCATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTActc                     <  2:262561/150‑1 (MQ=255)
                                                                                                                                aCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTActc                     <  2:766467/150‑1 (MQ=255)
                                                                                                                                 cGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTActct                    >  1:790973/1‑150 (MQ=255)
                                                                                                                                    gTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCtgtg                 >  1:674382/1‑150 (MQ=255)
                                                                                                                                     tAAAGTTCACGGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGGCAAACAAATTGGCACTGAATTATTTTACTCTgtgt                <  2:904299/150‑1 (MQ=255)
                                                                                                                                       aatttCACCTCCGAAATTTTTTCGGTGAATAACCACTTTTGAAATATTGTGACATATGTTTTGTCAAAATGTTCAACTTCTCTAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGtaa              <  2:3577/147‑1 (MQ=255)
                                                                                                                                       aaGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGTCACTGAATTATTTTACTCTGTGtaa              <  2:343502/150‑1 (MQ=255)
                                                                                                                                        tattCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGtaat             >  2:173150/3‑150 (MQ=255)
                                                                                                                                          ttCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGtaataa           >  2:1609929/1‑150 (MQ=255)
                                                                                                                                           tCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAATCGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTaataaa          >  2:734387/1‑150 (MQ=255)
                                                                                                                                            cACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTaataaat         >  2:540242/1‑150 (MQ=255)
                                                                                                                                            cACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGAACTGAATTATTTTACTCTGTGTaataaat         >  2:1325822/1‑150 (MQ=255)
                                                                                                                                            cACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTaataaat         >  2:1675949/1‑150 (MQ=255)
                                                                                                                                              ccgccgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTATGTGTAataaataa       >  1:936742/1‑150 (MQ=255)
                                                                                                                                                gccgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAg     >  1:645396/1‑150 (MQ=255)
                                                                                                                                                gccgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAg     <  2:228882/150‑1 (MQ=255)
                                                                                                                                                 ccgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAgg    <  2:674382/150‑1 (MQ=255)
                                                                                                                                                  cgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAggg   <  1:1244725/150‑1 (MQ=255)
                                                                                                                                                   gAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAATTGTGCAACTTCTCCAATGATCTCAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAGGgc  <  2:172806/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGTGTTGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAGGGC  >  NC_000913/1157482‑1157778

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