Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,770,203 T→A 5.1% intergenic (+17/‑169) ydiK → / → rprA UPF0118 family inner membrane protein/sRNA antisense activator for RpoS translation, Hfq‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,770,2030TA 5.1% 559.8 / 14.9 158intergenic (+17/‑169)ydiK/rprAUPF0118 family inner membrane protein/sRNA antisense activator for RpoS translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base T (80/70);  new base A (4/4);  total (84/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

Only 100 of 166 total aligned reads displayed.

TTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTC  >  NC_000913/1770055‑1770351
                                                                                                                                                    |                                                                                                                                                    
tttATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGct                                                                                                                                                      >  3:1267528/1‑149 (MQ=255)
  tATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGctct                                                                                                                                                    <  2:6835/149‑1 (MQ=255)
       gTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgc                                                                                                                                               >  4:1459435/1‑149 (MQ=255)
        tCCGGTGCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                                                                              <  3:1437349/149‑1 (MQ=255)
         ccGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCt                                                                                                                                             >  1:100440/1‑149 (MQ=255)
         ccGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCt                                                                                                                                             >  2:1171612/1‑149 (MQ=255)
           ggTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGa                                                                                                                                           <  1:1702197/149‑1 (MQ=255)
            gctCTGTTAGCCGTTTCCTGGCGTCTGTTTGCTGCGTGGGTGGAAGAAGTGCCGCCGCCTACGGACTAACTGGAAGGAGGTTGGGAAGAACTTGGCGAACGCGAGAAACAGAATAAGTAATTTCTCATCAGGCGGCTCAGACGCCTGaa                                                                                                                                          >  4:1059555/3‑148 (MQ=255)
             ttCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAtt                                                                                                                                         <  1:475112/149‑1 (MQ=255)
             ttCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAtt                                                                                                                                         <  3:786489/149‑1 (MQ=255)
             ttCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAtt                                                                                                                                         <  4:1106350/149‑1 (MQ=255)
              tCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtga                                                                                                                                        <  1:824569/149‑3 (MQ=255)
                 gTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTa                                                                                                                                     >  2:1723960/1‑149 (MQ=255)
                  ttAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTaa                                                                                                                                    <  1:424403/149‑1 (MQ=255)
                   tAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAac                                                                                                                                   <  2:1459305/149‑1 (MQ=255)
                      ccGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccac                                                                                                                                >  1:265410/1‑149 (MQ=255)
                             cTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAt                                                                                                                         >  2:836381/1‑149 (MQ=255)
                             cTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAt                                                                                                                         <  1:709585/149‑1 (MQ=255)
                              tGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAtt                                                                                                                        <  4:553033/149‑1 (MQ=255)
                              tGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAtt                                                                                                                        >  1:332237/1‑149 (MQ=255)
                                gCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                                                                      >  1:1291645/1‑149 (MQ=255)
                                   tCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGatat                                                                                                                   >  2:161546/1‑149 (MQ=255)
                                      gTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                                                                              <  3:1209134/119‑1 (MQ=255)
                                      gTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                                                                              >  4:1209134/1‑119 (MQ=255)
                                       tttGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgc                                                                                                               <  1:1617595/149‑1 (MQ=255)
                                           ccGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttgc                                                                                                           <  3:1059555/149‑1 (MQ=255)
                                                tgggtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTc                                                                                                      <  2:1172351/149‑1 (MQ=255)
                                                   gtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTcccc                                                                                                   >  1:1325724/1‑149 (MQ=255)
                                                   gtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTcccc                                                                                                   >  3:1627904/1‑149 (MQ=255)
                                                    tggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccac                                                                                                                                <  3:359647/119‑1 (MQ=255)
                                                    tggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccac                                                                                                                                >  4:359647/1‑119 (MQ=255)
                                                    tggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCa                                                                                                  >  1:372115/1‑149 (MQ=255)
                                                       aagaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCatat                                                                                               >  1:855873/1‑149 (MQ=255)
                                                       aagaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCatat                                                                                               >  3:1627431/1‑149 (MQ=255)
                                                           agTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTc                                                                                           <  3:778024/149‑1 (MQ=255)
                                                               ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGa                                                                                                                                           >  3:1156878/1‑97 (MQ=255)
                                                               ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGACTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGc                                                                                       >  2:452495/1‑149 (MQ=255)
                                                                  ccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTa                                                                                    <  4:577418/149‑1 (MQ=255)
                                                                       gACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAac                                                                                                                                   >  2:1066227/1‑97 (MQ=255)
                                                                       gACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAac                                                                                                                                   <  1:1066227/97‑1 (MQ=255)
                                                                       gACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGcc                                                                               >  3:210651/1‑149 (MQ=255)
                                                                       gACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGcc                                                                               >  2:814185/1‑149 (MQ=255)
                                                                            aCCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTc                                                                                           <  2:451506/132‑1 (MQ=255)
                                                                               aaCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCt                                                                       >  3:199859/1‑149 (MQ=255)
                                                                                aCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTg                                                                      >  4:25445/1‑149 (MQ=255)
                                                                                 ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTgg                                                                     <  4:760145/149‑1 (MQ=255)
                                                                                       gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc                                                               <  1:1183515/149‑1 (MQ=255)
                                                                                       gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc                                                               <  1:814185/149‑1 (MQ=255)
                                                                                        aaaTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc                                                              <  1:872014/149‑1 (MQ=255)
                                                                                           ttCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGt                                                           >  1:737668/1‑149 (MQ=255)
                                                                                            tctcGAAGAACTTGGCGAAATCCAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTc                                                          <  3:1034372/149‑1 (MQ=255)
                                                                                                gaagaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaa                                                      >  2:505138/1‑149 (MQ=255)
                                                                                                gaagaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaa                                                      >  2:621419/1‑149 (MQ=255)
                                                                                                 aagaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaaa                                                     >  4:189206/1‑149 (MQ=255)
                                                                                                 aagaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaaa                                                     <  2:855873/149‑1 (MQ=255)
                                                                                                 aagaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTACCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaaa                                                     >  1:1104660/1‑149 (MQ=255)
                                                                                                    aaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTc                                                  <  3:1096743/149‑1 (MQ=255)
                                                                                                     aCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCt                                                 <  2:267714/149‑1 (MQ=255)
                                                                                                      cTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCtt                                                <  1:45517/149‑1 (MQ=255)
                                                                                                       ttGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTcc                                                                                                     >  2:719288/1‑95 (MQ=255)
                                                                                                        tGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTa                                              <  3:25445/149‑1 (MQ=255)
                                                                                                            gAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTGCCTGTCGCAAATTCTTTACAGt                                          >  1:1476081/1‑149 (MQ=255)
                                                                                                            gAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGt                                          >  4:1644366/1‑149 (MQ=255)
                                                                                                             aaaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGtt                                         <  1:529153/149‑1 (MQ=255)
                                                                                                              aaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGttt                                        <  2:112894/149‑1 (MQ=255)
                                                                                                                tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGttttt                                      >  2:1025344/1‑149 (MQ=255)
                                                                                                                  gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTaa                                    >  3:1247658/1‑149 (MQ=255)
                                                                                                                   agaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTaaa                                   >  1:740064/1‑149 (MQ=255)
                                                                                                                   agaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTaaa                                   <  2:1325724/149‑1 (MQ=255)
                                                                                                                    gaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAAc                                  >  1:771764/1‑149 (MQ=255)
                                                                                                                     aaaCCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACt                                 >  1:910472/1‑149 (MQ=255)
                                                                                                                        ccGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCCCCGCTAATTGATATTTGCTTGCACTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAAt                              <  1:1171612/149‑1 (MQ=255)
                                                                                                                          gAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTGCTGCTTACCGCCATTCAGCTGATAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATgt                            >  1:584617/1‑148 (MQ=255)
                                                                                                                          gAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATga                            <  2:99819/149‑1 (MQ=255)
                                                                                                                            aTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATgaga                          <  3:1425756/149‑1 (MQ=255)
                                                                                                                            aTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATgaga                          >  1:403107/1‑149 (MQ=255)
                                                                                                                            aTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATgaga                          <  2:372115/149‑1 (MQ=255)
                                                                                                                             tAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGAc                         <  1:1025344/149‑1 (MQ=255)
                                                                                                                              aaGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTAAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGATACg                        >  2:1424959/1‑149 (MQ=255)
                                                                                                                                gTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGaa                      <  2:1104660/149‑1 (MQ=255)
                                                                                                                                    tttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATgaga                          <  2:1425952/141‑1 (MQ=255)
                                                                                                                                     ttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGa                 >  3:1681458/1‑149 (MQ=255)
                                                                                                                                       ctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATc               >  2:26804/1‑149 (MQ=255)
                                                                                                                                       ctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATc               <  1:1443328/149‑1 (MQ=255)
                                                                                                                                       ctcATCAGGCGGATCTGCCGCCTGCTTGTTAACCACCGCTCATTGATATTTGTTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTGGCAAATTCTTTACAGTTTTTAAACTAATGCGACGAATCTGAtt               <  2:584617/149‑2 (MQ=255)
                                                                                                                                        tcatcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCg              >  3:824532/1‑149 (MQ=255)
                                                                                                                                         catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGa             >  3:567952/1‑149 (MQ=255)
                                                                                                                                          atcaGGCGGCTCTGCCGCCTGATTGTTAACCGCCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGAACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGAc            >  4:1035941/1‑149 (MQ=255)
                                                                                                                                          atcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTc                                                                                           >  4:331322/1‑70 (MQ=255)
                                                                                                                                          atcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGAc            >  4:468390/1‑149 (MQ=255)
                                                                                                                                             aGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTc                                                                                           <  1:1580217/67‑1 (MQ=255)
                                                                                                                                             aGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTc                                                                                           >  2:1580217/1‑67 (MQ=255)
                                                                                                                                             aGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCa         >  3:354749/1‑149 (MQ=255)
                                                                                                                                               gcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAAt                              >  4:1318279/1‑126 (MQ=255)
                                                                                                                                               gcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAAt                              <  3:1318279/126‑1 (MQ=255)
                                                                                                                                                 ggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaaa                                                     >  1:578722/1‑101 (MQ=255)
                                                                                                                                                 ggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCaaaaa     <  1:1424959/149‑1 (MQ=255)
                                                                                                                                                 ggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGCTTTTAAACTAATGAGACGAATCTGATCGACGCaaaaa     <  1:1198016/149‑1 (MQ=255)
                                                                                                                                                   ctctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTa                                              <  2:562150/106‑1 (MQ=255)
                                                                                                                                                    tctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTc  >  3:1471077/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTATTGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTC  >  NC_000913/1770055‑1770351

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: