Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,115,849 C→G 9.3% S16S (TCG→TCC) ‡ wcaM ← colanic acid biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,115,8490CG 9.3% 339.1 / 19.8 107S16S (TCG→TCC) ‡wcaMcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (46/51);  new base G (5/5);  total (51/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 108 total aligned reads displayed.

CCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACC  >  NC_000913/2115701‑2115997
                                                                                                                                                     |                                                                                                                                                    
     ccGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGGATGGAGGAAGGCAAGCGCCGAGCtt                                                                                                                                                >  1:1222065/1‑150 (MQ=255)
      cGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTg                                                                                                                                               <  2:1244685/150‑1 (MQ=255)
       ggcaggcaCGACGACTGATTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTGGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGc                                                                                                                                              <  1:957618/150‑1 (MQ=255)
         caggccCGACGACTGTATGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGGGTGGGGTGTAGTCG‑CAAATGTTAACGCTATGGCGGGCGGGACGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCg                                                                                                                                            <  1:165817/150‑1 (MQ=255)
          aggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGcctt                                                                                                                                           <  2:1555091/150‑3 (MQ=255)
          aggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑GTAATGTTAACCCTTTGTCGGGCGGGAACCGCGCGGGCGAAACGGGTATTGAGGAAGGCAAGCGCCGAGCTTGCCGt                                                                                                                                           >  2:581034/1‑150 (MQ=255)
          aggcaCGACGACTGTTTGCCCTTCGCGGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCGTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGcctt                                                                                                                                           <  2:788963/150‑3 (MQ=255)
           ggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGCCTTc                                                                                                                                          <  2:873558/150‑1 (MQ=255)
             caCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGCCTTcct                                                                                                                                        <  2:1422048/150‑3 (MQ=255)
             caCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAg                                                                                                                                        >  1:990247/1‑150 (MQ=255)
             caCGACGACTGTTTGCCCTTCGCTGAAGGAGTGTTTAAATGAGGCGATCCAGTCGTGTTGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGCCTTcct                                                                                                                                        <  2:214845/150‑3 (MQ=255)
               cgacgaCTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTTGTCGTTTAATTTTAA‑GCGTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAACTCAACCGCCGAGCTTTCCGTCAGGa                                                                                                                                      >  2:1558494/1‑150 (MQ=255)
               cgacgaCTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGa                                                                                                                                      >  1:1556323/1‑150 (MQ=255)
                  cgacTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTTTTCG‑TTTATGTTAACGCTTTGTCGGGCGGGAATCGCGCGGGCGAAATGGGTATGGAGGAAGGCAATCGCCGAGCGTGCCGGCAGGAACg                                                                                                                                   >  1:475275/1‑150 (MQ=255)
                    acTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTg                                                                                                                                 <  2:1621286/150‑1 (MQ=255)
                    acTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTg                                                                                                                                 <  2:1109476/150‑1 (MQ=255)
                    acTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTg                                                                                                                                 <  1:1705454/150‑1 (MQ=255)
                     cTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGc                                                                                                                                <  1:757372/150‑1 (MQ=255)
                      tGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGCGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCg                                                                                                                               <  1:974208/150‑1 (MQ=255)
                      tGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCg                                                                                                                               >  1:10332/1‑150 (MQ=255)
                       gTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGt                                                                                                                              >  2:623158/1‑150 (MQ=255)
                       gTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGt                                                                                                                              >  1:1186395/1‑150 (MQ=255)
                       gTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCAAGTCGTGTGGGTTGTAGTCG‑TTAATTTTAACGCTTTGTCGGGCGGGAAGCGCCCGGGCGAAATGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGcct                                                                                                                              >  2:612636/1‑148 (MQ=255)
                          tGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGGTGTAGTCG‑TTAATTTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCgg                                                                                                                           >  1:1318952/1‑150 (MQ=255)
                           gCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTATTAC‑TTAATTTTTAAGCTTTTTCGGGCGGGAAGCGCGCGGGCGAAAAGGGTATTGATGAACGCAACCCCCGCGCATGCCGTCACGAACCTGCGTCggg                                                                                                                          >  2:663306/1‑150 (MQ=255)
                           gCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCggg                                                                                                                          <  2:1187051/150‑1 (MQ=255)
                           gCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCggg                                                                                                                          >  2:1150710/1‑150 (MQ=255)
                            cccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAAATTTAACGCTTTGTCGGGCCGGAAGCGCGCGGGCCAAACGGGTATGGAGGAAGGCAAGCCCCGAGCTTGACCTCAGGAACCTGCCTCCGga                                                                                                                         >  2:430850/1‑150 (MQ=255)
                               ttCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGagag                                                                                                                      <  2:1375903/150‑1 (MQ=255)
                                 cGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAAAGCTTTGTAGGGCGGGAAGCGCGCGGGCGAAAGGGGTTAGGAGGAAGGCAAGCGCCGAGCTTTCCGTCAGGAACGTGCGTCGGGAGAGtt                                                                                                                    >  1:181721/1‑150 (MQ=255)
                                   cTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGtttt                                                                                                                  >  2:1223528/1‑150 (MQ=255)
                                    tGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGttttt                                                                                                                 <  1:349355/150‑1 (MQ=255)
                                      aaGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTa                                                                                                               <  2:570474/150‑1 (MQ=255)
                                         gCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAAt                                                                                                            <  1:1288159/150‑1 (MQ=255)
                                         gCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAAt                                                                                                            <  1:581034/150‑1 (MQ=255)
                                                     gAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTtct                                                                                                <  1:1106404/150‑1 (MQ=255)
                                                     gAGCCGATCCAGACGTGTGGGTTGTAGTCG‑TTAAAGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTtct                                                                                                <  2:1248664/150‑1 (MQ=255)
                                                           ttCCTGTCGTGTGGGTTGTAGTCG‑TTAATGTTTACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCta                                                                                          <  1:1280654/149‑1 (MQ=255)
                                                             ccAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtata                                                                                        <  2:1556323/150‑1 (MQ=255)
                                                                  cGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGcc                                                                                   <  1:618064/150‑1 (MQ=255)
                                                                    tgtgGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTg                                                                                 >  2:14143/1‑150 (MQ=255)
                                                                     gtgGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgc                                                                                <  1:1558494/150‑1 (MQ=255)
                                                                     gtgGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGGATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgc                                                                                <  2:1162878/150‑1 (MQ=255)
                                                                      tgGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgca                                                                               <  2:1302368/150‑1 (MQ=255)
                                                                       gggTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAATGGGGTATGGAGGATGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgcag                                                                              >  1:235083/1‑150 (MQ=255)
                                                                       gggTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgcag                                                                              >  1:1237920/1‑150 (MQ=255)
                                                                       gggTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTATTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgcag                                                                              >  1:131037/1‑150 (MQ=255)
                                                                         gTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgcagca                                                                            <  2:548612/150‑1 (MQ=255)
                                                                            gTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGc                                                                         <  1:1199602/150‑1 (MQ=255)
                                                                            gTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGc                                                                         <  2:1207228/150‑1 (MQ=255)
                                                                             tAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCt                                                                        >  1:1434821/1‑150 (MQ=255)
                                                                              aGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTg                                                                       >  1:435211/1‑150 (MQ=255)
                                                                                 cG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCg                                                                    >  2:984809/1‑150 (MQ=255)
                                                                                 cG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCg                                                                    >  2:565656/1‑150 (MQ=255)
                                                                                    ttAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAg                                                                  <  2:823297/150‑1 (MQ=255)
                                                                                       aTGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtc                                                               <  1:409258/150‑1 (MQ=255)
                                                                                        tGTTATCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGTAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtct                                                              <  1:430850/150‑1 (MQ=255)
                                                                                        tGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtct                                                              >  1:351208/1‑150 (MQ=255)
                                                                                         gTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                             <  1:1150710/150‑1 (MQ=255)
                                                                                           tAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGa                                                           >  2:1449300/1‑150 (MQ=255)
                                                                                            aaCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGAt                                                          >  2:1602121/1‑150 (MQ=255)
                                                                                             aCGCTTTGTGGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGAtt                                                         <  2:351208/150‑1 (MQ=255)
                                                                                                cTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                      >  2:1158885/1‑150 (MQ=255)
                                                                                                cTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                      >  2:545737/1‑150 (MQ=255)
                                                                                                cTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                      >  1:1148291/1‑150 (MQ=255)
                                                                                                cTTTGTCGGGCGGAAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                      >  2:331216/1‑150 (MQ=255)
                                                                                                 tttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATc                                                     >  2:1517297/1‑150 (MQ=255)
                                                                                                  ttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTGTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCa                                                    >  1:817616/1‑150 (MQ=255)
                                                                                                  ttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCa                                                    <  2:1318952/150‑1 (MQ=255)
                                                                                                   tggcgggcgggAAGCGCGCGGGCGAAAGGGGTATGTAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTATATGGCATGCGTTCTACTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCAc                                                   <  2:1339271/147‑1 (MQ=255)
                                                                                                   tGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCGGGCGAGTTCTCGATTGATCAc                                                   >  2:152819/1‑150 (MQ=255)
                                                                                                     tCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACct                                                 >  2:592847/1‑150 (MQ=255)
                                                                                                     tCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATACGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACct                                                 >  2:563210/1‑150 (MQ=255)
                                                                                                      cgggcgggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                <  2:235083/150‑1 (MQ=255)
                                                                                                        ggcgggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTAAAATGGCATGCGTTCTCCTCTAAAAAGCCTGCAGCAAGCTGGCGAGATCTCGATTGATCACctgct                                              <  2:459764/150‑1 (MQ=255)
                                                                                                          cgggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTgg                                            >  1:590379/1‑150 (MQ=255)
                                                                                                          cgggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTgg                                            <  1:1165177/150‑1 (MQ=255)
                                                                                                            ggAAGCGCGCGGGCGAAAGTGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGtt                                          <  1:894805/150‑1 (MQ=255)
                                                                                                             gAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTa                                         <  1:1602121/150‑1 (MQ=255)
                                                                                                                 cgcgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAAt                                     <  2:952965/150‑1 (MQ=255)
                                                                                                                     ggggCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtg                                 <  1:518843/149‑1 (MQ=255)
                                                                                                                     cgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtg                                 <  1:1646806/150‑1 (MQ=255)
                                                                                                                       ggCGAAAAGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtgat                               >  1:1331055/1‑148 (MQ=255)
                                                                                                                        gCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTc                              <  2:435211/150‑1 (MQ=255)
                                                                                                                         cGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCa                             >  1:1187558/1‑150 (MQ=255)
                                                                                                                              ggggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttt                        >  2:1032922/1‑150 (MQ=255)
                                                                                                                              ggggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttt                        <  1:565656/150‑1 (MQ=255)
                                                                                                                              ggggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttt                        <  2:1508978/150‑1 (MQ=255)
                                                                                                                                 gTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTc                     <  1:331216/150‑1 (MQ=255)
                                                                                                                                 gTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTc                     <  2:1222065/150‑1 (MQ=255)
                                                                                                                                  tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcg                    <  2:475275/150‑1 (MQ=255)
                                                                                                                                    tGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcg                  >  1:1315751/1‑150 (MQ=255)
                                                                                                                                     ggaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgc                 <  1:984809/150‑1 (MQ=255)
                                                                                                                                      cctgacggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcg                >  1:1337787/7‑150 (MQ=255)
                                                                                                                                       ctgacggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcgc               >  1:33786/6‑150 (MQ=255)
                                                                                                                                         gacggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcgcgc             >  2:650736/4‑150 (MQ=255)
                                                                                                                                          acggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcgcgcg            >  1:1370043/3‑150 (MQ=255)
                                                                                                                                             gCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGttt         >  1:123224/1‑150 (MQ=255)
                                                                                                                                             gCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGAGCGCGCGttt         >  1:340641/1‑150 (MQ=255)
                                                                                                                                              cAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTg        <  2:590379/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCG‑TTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACC  >  NC_000913/2115701‑2115997

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