Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,157,636 G→C 10.1% intergenic (+62/‑233) 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,6360GC10.1% 297.2 / 23.6 109intergenic (+62/‑233)ycfH/ptsGputative DNase/fused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base G (44/52);  new base C (5/6);  total (49/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-02

Only 100 of 110 total aligned reads displayed.

TGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAGGGC  >  NC_000913/1157487‑1157778
                                                                                                                                                     |                                                                                                                                              
tGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTg                                                                                                                                                <  2:420817/150‑1 (MQ=255)
   cGTTGAAGAACTGGCTCAGTTAACCAACGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTTAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGc                                                                                                                                             <  1:905252/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)
       gaagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTg                                                                                                                                         <  2:642585/150‑1 (MQ=255)
       gaagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTg                                                                                                                                         <  2:568888/150‑1 (MQ=255)
       gaagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTg                                                                                                                                         <  2:828678/150‑1 (MQ=255)
       gaagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGACGAAAATTGGGCGGTg                                                                                                                                         >  2:1335430/1‑150 (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)
                  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)
                             ccGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtg                                                                                                                   <  2:863269/150‑1 (MQ=255)
                               gATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCTGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtgtg                                                                                                                 <  1:1413315/150‑1 (MQ=255)
                               gATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtgtg                                                                                                                 <  1:216658/150‑1 (MQ=255)
                               gATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtgtg                                                                                                                 <  2:1655579/150‑1 (MQ=255)
                                 tAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGAc                                                                                                               <  1:1436475/150‑1 (MQ=255)
                                 tAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGAc                                                                                                               <  1:1090152/150‑1 (MQ=255)
                                 tAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGAc                                                                                                               <  2:934893/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)
                                                   tCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtg                                                                                             >  1:887795/1‑150 (MQ=255)
                                                     acaTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtgtg                                                                                           <  2:593992/150‑1 (MQ=255)
                                                     acaTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtgtg                                                                                           >  1:1069407/1‑150 (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                                                                               >  2:549971/1‑150 (MQ=255)
                                                                 cccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCaa                                                                               >  1:1695498/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)
                                                                          aaTCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGaa                                                                      <  2:825057/150‑1 (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                                                    <  1:493827/150‑1 (MQ=255)
                                                                                            gTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCg                                                    <  2:880401/150‑1 (MQ=255)
                                                                                            gTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCg                                                    >  2:1411373/1‑150 (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)
                                                                                                          cGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTTGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGc                                      <  1:1513422/150‑1 (MQ=255)
                                                                                                           gTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTCAAATATTGTGACAAATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGTCa                                     >  1:72343/1‑150 (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:1144241/1‑150 (MQ=255)
                                                                                                                      cTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTAttt                          <  2:1069407/150‑1 (MQ=255)
                                                                                                                        aaaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTa                        >  1:1429297/1‑150 (MQ=255)
                                                                                                                        aaaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTa                        >  1:766467/1‑150 (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)
                                                                                                                                       cACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCATTGATCTGAAGTTGAAACGTGATAGCCGTCAAACTAATTGGCACTGTATTATTTTAGTCTGTGTaataaat         >  2:1123394/1‑150 (MQ=255)
                                                                                                                                       cACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTaataaat         >  1:1382506/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     <  2:228882/150‑1 (MQ=255)
                                                                                                                                           gccgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAg     >  1:645396/1‑150 (MQ=255)
                                                                                                                                            ccgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAgg    <  2:674382/150‑1 (MQ=255)
                                                                                                                                             cgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAggg   <  1:1244725/150‑1 (MQ=255)
                                                                                                                                              gAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAATTGTGCAACTTCTCCAATGATCTCAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAGGgc  <  2:172806/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAGGGC  >  NC_000913/1157487‑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: