Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,153,794 C→A 5.7% intergenic (+15/‑222) sibB → / → mdtA sRNA antisense regulator of toxic IbsB protein/multidrug efflux system, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,153,7940CA 5.7% 285.7 / 5.1 84intergenic (+15/‑222)sibB/mdtAsRNA antisense regulator of toxic IbsB protein/multidrug efflux system, subunit A
Reads supporting (aligned to +/- strand):  ref base C (42/37);  new base A (2/3);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

TAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAATTCCTC  >  NC_000913/2153651‑2153942
                                                                                                                                               |                                                                                                                                                    
tAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTg                                                                                                                                                 <  3:170229/149‑1 (MQ=255)
 agagCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGa                                                                                                                                                <  3:371064/149‑1 (MQ=255)
    ggcgggTTTCCTCCGCCCTGCTAGTCTTAGTAGGGGAGCAAGCCTATGACTAACAGCACCACGATGATGCGTAGGTTCATCATGACCCGTTCCTTATTCATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAACGGGACGCTGaaac                                                                                                                                             <  2:751250/146‑2 (MQ=21)
    gCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGGGACGCTGaaac                                                                                                                                             <  3:163010/149‑2 (MQ=255)
              ccccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTa                                                                                                                                   >  2:77501/1‑149 (MQ=255)
              ccccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTa                                                                                                                                   >  3:105596/1‑149 (MQ=255)
               cccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTAc                                                                                                                                  <  4:192091/149‑1 (MQ=255)
               cccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTACGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTAc                                                                                                                                  <  4:494540/149‑1 (MQ=255)
                   cccTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTAcctcc                                                                                                                              <  1:756326/149‑1 (MQ=255)
                       ggTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTc                                                                                                                          >  4:93740/1‑149 (MQ=255)
                         tAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGc                                                                                                                        <  1:167221/149‑1 (MQ=255)
                          aGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGcc                                                                                                                       <  1:669756/149‑1 (MQ=255)
                          aGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGcc                                                                                                                       <  4:447477/149‑1 (MQ=255)
                            tCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCAt                                                                                                                     <  3:622700/149‑1 (MQ=255)
                            tCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCAt                                                                                                                     >  2:508503/1‑149 (MQ=255)
                              ttAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGc                                                                                                                   <  2:155912/149‑1 (MQ=255)
                              ttAGTAAGCGGGGAAACTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCACTTCCTCGGGAGGGGCTTTCCCGTCTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGc                                                                                                                   >  4:513305/1‑149 (MQ=255)
                                aGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCaa                                                                                                                 <  3:639502/149‑1 (MQ=255)
                                   aaGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCa                                                                                                              <  2:528148/149‑1 (MQ=255)
                                          gAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATc                                                                                                       <  4:534694/149‑1 (MQ=255)
                                            aGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCGTGACCCTTTGCTTATTTATTGCTCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCTGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGc                                                                                                     <  3:10487/149‑1 (MQ=39)
                                               ttATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCAt                                                                                                                     >  4:122403/1‑130 (MQ=255)
                                               ttATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCAt                                                                                                                     <  3:122403/130‑1 (MQ=255)
                                               ttATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTaa                                                                                                  >  4:668346/1‑149 (MQ=255)
                                               ttATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGGGACGCTGaaacg                                                                                                                                            >  3:180933/1‑105 (MQ=25)
                                               ttATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGGGACGCTGaaacg                                                                                                                                            <  4:180933/107‑3 (MQ=25)
                                                   gACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAAcgcg                                                                                              >  2:344005/1‑149 (MQ=255)
                                                    aCTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGt                                                                                             <  4:186346/149‑1 (MQ=255)
                                                            caccacGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATaaaa                                                                                     <  2:582254/149‑1 (MQ=255)
                                                            caccacGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATaaaa                                                                                     <  4:284414/149‑1 (MQ=255)
                                                              ccacGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATa                                                                                   <  4:347452/149‑1 (MQ=255)
                                                                   atgatgAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTc                                                                              <  4:18970/149‑1 (MQ=255)
                                                                    tgatgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTcc                                                                             >  4:93167/1‑149 (MQ=255)
                                                                    tgatgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTcc                                                                             >  3:477923/1‑149 (MQ=255)
                                                                     gatgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCt                                                                            >  2:255994/1‑149 (MQ=255)
                                                                       tgaGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCtgt                                                                          <  3:252608/149‑1 (MQ=255)
                                                                          gTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGtt                                                                       >  3:208192/1‑149 (MQ=255)
                                                                           tAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGtta                                                                      >  1:145428/1‑149 (MQ=255)
                                                                            aGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGttat                                                                     <  1:508503/149‑1 (MQ=255)
                                                                             gCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGttatt                                                                    <  4:105596/149‑1 (MQ=255)
                                                                              cTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGttatta                                                                   >  2:678655/1‑149 (MQ=255)
                                                                                tcatcaTGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACt                                                                 <  4:585678/149‑1 (MQ=255)
                                                                                   tcatGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGa                                                              >  2:500451/1‑149 (MQ=255)
                                                                                    catGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTAc                                                                  <  1:165032/144‑1 (MQ=255)
                                                                                    catGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTAc                                                                  >  2:165032/1‑144 (MQ=255)
                                                                                      tGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGc                                                           >  2:86908/1‑149 (MQ=255)
                                                                                      tGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGc                                                           >  2:371712/1‑149 (MQ=255)
                                                                                        aCCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGt                                                         >  4:650784/1‑149 (MQ=255)
                                                                                          ccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGc                                                       >  1:234150/1‑149 (MQ=255)
                                                                                          ccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGc                                                       >  2:131759/1‑149 (MQ=255)
                                                                                          ccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGc                                                       >  3:283670/1‑149 (MQ=255)
                                                                                            tttCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTc                                                     >  1:625404/1‑149 (MQ=255)
                                                                                            tttCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTc                                                     <  4:734495/149‑1 (MQ=255)
                                                                                                 ttatttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCaaa                                                >  1:606070/1‑149 (MQ=255)
                                                                                                    tttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAAtgt                                             >  2:665013/1‑149 (MQ=255)
                                                                                                     ttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAAtgtg                                            <  1:183898/149‑1 (MQ=255)
                                                                                                     ttatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAAcgtg                                            >  2:533177/1‑149 (MQ=255)
                                                                                                      tatGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGc                                           <  3:93167/149‑1 (MQ=255)
                                                                                                        tGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGccc                                         >  1:567942/1‑149 (MQ=255)
                                                                                                        tGGCCCCTTCCTCGGGAGGGGCCTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGccc                                         <  1:665013/149‑1 (MQ=255)
                                                                                                         ggCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCg                                        <  3:279329/149‑1 (MQ=255)
                                                                                                          gCCCCTCCCGAGGAAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTGGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGt                                       >  1:53562/1‑149 (MQ=21)
                                                                                                           ccccTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTc                                      >  2:384345/1‑149 (MQ=255)
                                                                                                            cccTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCa                                     >  2:543982/1‑149 (MQ=255)
                                                                                                              cTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACt                                                                 >  2:539500/1‑119 (MQ=255)
                                                                                                              cTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACt                                                                 <  1:539500/119‑1 (MQ=255)
                                                                                                                tCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCa                                 >  4:689967/1‑149 (MQ=255)
                                                                                                                  cTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGa                               <  1:77501/149‑1 (MQ=255)
                                                                                                                   tCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGGCATTCAGAc                              >  3:639238/1‑149 (MQ=255)
                                                                                                                   tCGGGAGGGGCTTTCCCGTATCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGAc                              >  3:626373/1‑149 (MQ=255)
                                                                                                                    cGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACg                             <  1:371712/149‑1 (MQ=255)
                                                                                                                      ggAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGAt                           <  3:93740/149‑1 (MQ=255)
                                                                                                                          gggCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCa                       >  3:61494/1‑149 (MQ=255)
                                                                                                                           ggCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTAc                                                                                                                                  <  2:501593/41‑1 (MQ=40)
                                                                                                                            gCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGa                     <  2:606070/149‑1 (MQ=255)
                                                                                                                              tttCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACa                   >  4:179905/1‑149 (MQ=255)
                                                                                                                                  ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAgtgt               >  4:20230/1‑149 (MQ=255)
                                                                                                                                  ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAgtgt               >  2:51886/1‑149 (MQ=255)
                                                                                                                                  ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAgtgt               >  3:578561/1‑149 (MQ=255)
                                                                                                                                  ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAgtgt               >  3:450439/1‑149 (MQ=255)
                                                                                                                                  ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAgtgt               >  3:447582/1‑149 (MQ=255)
                                                                                                                                   gatttCAGCGGGACGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGCTGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGtt              >  3:382800/3‑149 (MQ=255)
                                                                                                                                         aGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATaa        <  3:20230/149‑1 (MQ=255)
                                                                                                                                             tCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAATtcc    >  2:712935/1‑149 (MQ=255)
                                                                                                                                               ccGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAATtcctc  <  2:234150/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAATTCCTC  >  NC_000913/2153651‑2153942

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