Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,153,792 T→G 5.9% intergenic (+13/‑224) 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,7920TG 5.9% 288.1 / 8.3 84intergenic (+13/‑224)sibB/mdtAsRNA antisense regulator of toxic IbsB protein/multidrug efflux system, subunit A
Reads supporting (aligned to +/- strand):  ref base T (42/37);  new base G (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 = 8.28e-01

GTGAGGGTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAATTCC  >  NC_000913/2153644‑2153940
                                                                                                                                                    |                                                                                                                                                    
gTGAGGGTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGt                                                                                                                                                      >  1:528148/1‑149 (MQ=255)
 tGAGGGTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTc                                                                                                                                                     >  3:186346/1‑149 (MQ=255)
       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                                                                                                                                 >  3:105596/1‑149 (MQ=255)
                     ccccGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTa                                                                                                                                 >  2:77501/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                                                                                                                   <  3:122403/130‑1 (MQ=255)
                                                      ttATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCAt                                                                                                                   >  4:122403/1‑130 (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                                                                                   <  4:284414/149‑1 (MQ=255)
                                                                   caccacGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATaaaa                                                                                   <  2:582254/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                                                                >  2:165032/1‑144 (MQ=255)
                                                                                           catGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTAc                                                                <  1:165032/144‑1 (MQ=255)
                                                                                             tGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGc                                                         >  2:86908/1‑149 (MQ=255)
                                                                                             tGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGc                                                         >  2:371712/1‑149 (MQ=255)
                                                                                               aCCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGt                                                       >  4:650784/1‑149 (MQ=255)
                                                                                                 ccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGc                                                     >  2:131759/1‑149 (MQ=255)
                                                                                                 ccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGc                                                     >  1:234150/1‑149 (MQ=255)
                                                                                                 ccTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGc                                                     >  3:283670/1‑149 (MQ=255)
                                                                                                   tttCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTc                                                   <  4:734495/149‑1 (MQ=255)
                                                                                                   tttCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTc                                                   >  1:625404/1‑149 (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             >  3:450439/1‑149 (MQ=255)
                                                                                                                                         ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAgtgt             >  4:20230/1‑149 (MQ=255)
                                                                                                                                         ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAgtgt             >  3:447582/1‑149 (MQ=255)
                                                                                                                                         ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAgtgt             >  2:51886/1‑149 (MQ=255)
                                                                                                                                         ccGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAgtgt             >  3:578561/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)
                                                                                                                                                    |                                                                                                                                                    
GTGAGGGTAGAGCGGGGTTTCCCCCGCCCTGGTAGTCTTAGTAAGCGGGGAAGCTTATGACTAAGAGCACCACGATGATGAGTAGCTTCATCATGACCCTTTCCTTATTTATGGCCCCTTCCTCGGGAGGGGCTTTCCCGTTTCAGCGTCCCGCTGAAATCGTCGGCTTACCTCCTTTCGCCATGCAAGCAGTCTATCGCTAACGCGTAGATAAAATAGTTTCCTGTGTTATTACTGGATGCGTGCTCGCAAATGTGCCCGTCATTCAGACGATTCCAGACAGTGTTTCATAATTCC  >  NC_000913/2153644‑2153940

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: