Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,115,851 A→C 7.4% S16A (TCG→GCG) ‡ wcaM ← colanic acid biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,115,8510AC 7.4% 364.1 / 16.6 132S16A (TCG→GCG) ‡wcaMcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (63/59);  new base C (4/6);  total (67/65)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

Only 100 of 135 total aligned reads displayed.

CAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGA  >  NC_000913/2115702‑2116000
                                                                                                                                                     |                                                                                                                                                     
cAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGCAGCGCGCGGGGGAAAGGGGTATGGAGGAAGGCAAGCGCCGa                                                                                                                                                       >  1:1101635/1‑150 (MQ=255)
  aTCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGc                                                                                                                                                     <  1:1143241/150‑1 (MQ=255)
      ggcaggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGGATGGAGGAAGGCAAGCGCCGAGCTTGc                                                                                                                                                 >  2:1185859/1‑150 (MQ=255)
        caggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGAGCCAGTCGTGTGGGTGGTAGTCGCTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCg                                                                                                                                               <  1:100232/150‑1 (MQ=255)
         aggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGcctt                                                                                                                                              <  2:1283563/150‑3 (MQ=255)
         aggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGcctt                                                                                                                                              <  1:1576312/150‑3 (MQ=255)
         aggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGt                                                                                                                                              <  1:40077/150‑1 (MQ=255)
          ggcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTc                                                                                                                                             >  1:196537/1‑150 (MQ=255)
           gcaCGACGAGTGTTTGCCCTTCGCTGATGGATTGTTTAACTGAGGCTATCCAGTCGTGTGGGTTGTAGTCGTTTATCTTAAGGCTTTGTCGGGCGGGAAGCGCGCGGGCGTAAGGGGTATGGAGGAAGGCATGCTCGGCGCTTGCCTTcc                                                                                                                                            <  2:753768/150‑2 (MQ=255)
           gcaCGACGACTGTTTGCCCTTCTCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGCCTTcc                                                                                                                                            <  1:1255540/150‑2 (MQ=255)
           gcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGCCTTcc                                                                                                                                            <  1:1395393/150‑2 (MQ=255)
           gcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCa                                                                                                                                            <  2:200151/150‑1 (MQ=255)
            caCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAg                                                                                                                                           <  1:971314/150‑1 (MQ=255)
            caCGACGACTGGTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGCCTTcct                                                                                                                                           <  2:1503880/150‑3 (MQ=255)
                acgacTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTTATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATTGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAAc                                                                                                                                       >  1:1147317/1‑150 (MQ=255)
                       tttGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACATGCGTc                                                                                                                                <  2:551462/150‑1 (MQ=255)
                         tGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGGTGTAGTCGTTAATGTTAACGCATTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTAgg                                                                                                                              >  1:162247/1‑150 (MQ=255)
                          gCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCggg                                                                                                                             <  2:207975/150‑1 (MQ=255)
                          gCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGGGTGGGTTTTAGGCGGTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCggg                                                                                                                             <  1:1463269/150‑1 (MQ=255)
                                 gCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGttt                                                                                                                      <  1:253368/150‑1 (MQ=255)
                                   tgcaGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGGTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGttttt                                                                                                                    <  1:370424/147‑1 (MQ=255)
                                    gAAGGCTTGTTTAAATGAGGCGATCCAGGCGTGTGGGGTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGtttttt                                                                                                                   <  1:318346/150‑1 (MQ=255)
                                     caGGGTTGTTTAAATGAGGCGAGCCAGGCGTGTGGGGGGTTGTGGTTATGGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTa                                                                                                                  <  1:295584/149‑1 (MQ=255)
                                         cTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGATATGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATg                                                                                                              >  2:2728/1‑150 (MQ=255)
                                         cTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATg                                                                                                              >  1:35079/1‑150 (MQ=255)
                                           tGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGc                                                                                                            <  2:532560/150‑1 (MQ=255)
                                            gTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTAGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGTAAGGAAAGCGCCGAGCTTGCCGTAAGGAACGTGCGTCGGTAGAGTTTTTTAAATGGCa                                                                                                           >  1:1281289/1‑150 (MQ=255)
                                            gTTTAAATGAGGCGATCCAGTCGTGTGGGTAGTAGTCGTTAATGTTTACGCTTTGTCGGGCGTGAATCGCGCGGGCGTAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCa                                                                                                           <  2:348076/150‑1 (MQ=255)
                                             tttAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTGGGTTATGTTAACGCGTTGTGGGGCGGGCAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAATCGCCGAGCTTGCCGTGAGGGGCGGGCGGGGGGGGAGTTTTTTAAATGGcct                                                                                                          >  2:1581685/1‑148 (MQ=255)
                                                aaaTGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCg                                                                                                       >  2:713855/1‑150 (MQ=255)
                                                      ggCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGTTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTctcc                                                                                                 <  1:935894/150‑1 (MQ=255)
                                                      ggCGATCCAGTCGTGGGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTctcc                                                                                                 <  2:1532904/150‑1 (MQ=255)
                                                       gCGAACCAGTCGTGGGGGGTGTAGTCGGTAATGTTTACGCGGTGTCGGGGGGGAAGCGCGCGGGCGAAAGGGGGATGGAGGAAGGCCAGCGCCGCGCTTGGCGGCCGGAACGGGCGTGGGGTGAGTTTTTTAAATGGCAAGGGTTgtccg                                                                                                >  2:1749418/1‑149 (MQ=255)
                                                         gATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCctct                                                                                              >  1:1010680/1‑150 (MQ=255)
                                                          aTCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCta                                                                                             >  2:41487/1‑150 (MQ=255)
                                                          aTCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACCCTTTGTAGGGCGGGAAGCCCGCGGGCGAAAAGGGTATTGAGGAAGGCAAACGCCGAGCATGCCGTAAAGAACGTGCGTCGGGAGAGTTTTTTAAAATGCATGCGGTCTACTCta                                                                                             >  1:1375224/1‑150 (MQ=255)
                                                          aTCCAGTCGTGGGGGTTGTAGTCGGTAATGTTCCCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCta                                                                                             <  2:249588/150‑1 (MQ=255)
                                                           tCCAGTCGTGTGGGTTGTATTCGTTACTGTTAACGATTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtat                                                                                            <  1:341591/150‑1 (MQ=255)
                                                           tCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtat                                                                                            <  1:1713211/150‑1 (MQ=255)
                                                           tCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtat                                                                                            <  1:640157/150‑1 (MQ=255)
                                                                 cGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGcc                                                                                      <  1:1276872/150‑1 (MQ=255)
                                                                 cGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGcc                                                                                      >  1:479997/1‑150 (MQ=255)
                                                                 cGTGGGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGAGCGCGGGCGAAAGGGGGATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGcc                                                                                      <  2:870621/150‑1 (MQ=255)
                                                                    gtgGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgc                                                                                   >  1:312875/1‑150 (MQ=255)
                                                                     tgGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgca                                                                                  <  2:1070032/150‑1 (MQ=255)
                                                                     tgGGTTGTAGTCGATAATGTTAACGCTGTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgca                                                                                  <  2:1114057/150‑1 (MQ=255)
                                                                      gggTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGACCGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgcag                                                                                 >  2:1169048/1‑150 (MQ=255)
                                                                       ggTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTgcagc                                                                                <  2:1466926/150‑1 (MQ=255)
                                                                           gTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTTTGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGGCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGc                                                                            >  2:1131910/1‑150 (MQ=255)
                                                                           gTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGc                                                                            >  2:835718/1‑150 (MQ=255)
                                                                            tAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCt                                                                           <  2:196537/150‑1 (MQ=255)
                                                                                cGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCg                                                                       <  1:645238/150‑1 (MQ=255)
                                                                                cGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCg                                                                       <  1:1411464/150‑1 (MQ=255)
                                                                                 gTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATTGAGGAAGGCAAGCGCCGAGCTTGCAGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGa                                                                      >  1:900695/1‑150 (MQ=255)
                                                                                  ttAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAg                                                                     <  1:1016296/150‑1 (MQ=255)
                                                                                      tGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtct                                                                 >  1:353834/1‑150 (MQ=255)
                                                                                       gTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGGCTGCAGCAAGCTGGAGAGTtcta                                                                >  1:110733/1‑149 (MQ=255)
                                                                                       gTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                <  1:1298117/150‑1 (MQ=255)
                                                                                        ttAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCg                                                               >  1:1778320/1‑150 (MQ=255)
                                                                                        ttAACGCGTTGTCGGGCGGTAATCGCGCGGGCGGAAGGGGGATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGTAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCg                                                               <  2:1676485/150‑1 (MQ=255)
                                                                                             gCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGa                                                          >  2:1277043/1‑150 (MQ=255)
                                                                                              cTTTGTCGGGCGGGAAGCTCGCGTGCGAAAGGGGTTTGGTCGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGGTTTTTAAATGGCATGCGGTCTCCGCTATAAATCCTGCAGCATGCTGGCGAGTTCTCGATTGAt                                                         <  1:120566/150‑1 (MQ=255)
                                                                                               tttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTGTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAta                                                        >  1:288956/1‑149 (MQ=255)
                                                                                               tttGTCGGGAGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATc                                                        >  2:1066369/1‑150 (MQ=255)
                                                                                                ttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCa                                                       <  2:162247/150‑1 (MQ=255)
                                                                                                  gTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCCGAGCTTGCAGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCAcc                                                     >  2:1199985/1‑150 (MQ=255)
                                                                                                      ggcgggAAGCGTGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGAGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgcg                                                 >  2:347838/1‑149 (MQ=255)
                                                                                                      ggcgggAAGCGTGCGGGCGAAAGGGGTATCGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGAGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgct                                                 >  2:485766/1‑150 (MQ=255)
                                                                                                      ggcgggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgct                                                 >  2:431329/1‑150 (MQ=255)
                                                                                                        cgggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTgg                                               <  2:100553/150‑1 (MQ=255)
                                                                                                         gggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGt                                              >  2:800948/1‑150 (MQ=255)
                                                                                                         gggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGt                                              >  1:413419/1‑150 (MQ=255)
                                                                                                          ggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGtt                                             >  2:29449/1‑150 (MQ=255)
                                                                                                          ggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGAAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGtt                                             >  1:1078378/1‑150 (MQ=255)
                                                                                                           gAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTa                                            >  2:1774120/1‑150 (MQ=255)
                                                                                                             aGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACATGCTGGTTaaa                                          >  2:860430/1‑150 (MQ=255)
                                                                                                               cgcgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAAt                                        <  2:646462/150‑1 (MQ=255)
                                                                                                                   cgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtg                                    >  2:344041/1‑150 (MQ=255)
                                                                                                                       cGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCa                                >  2:26390/1‑150 (MQ=255)
                                                                                                                       cGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCa                                >  2:296913/1‑150 (MQ=255)
                                                                                                                            ggggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttt                           >  1:538822/1‑150 (MQ=255)
                                                                                                                            ggggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttt                           >  1:29190/1‑150 (MQ=255)
                                                                                                                            ggggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttt                           >  1:531888/1‑150 (MQ=255)
                                                                                                                              ggTATGGAGGAAGGCCAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACttttt                         <  1:181594/150‑1 (MQ=255)
                                                                                                                                 aTGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgc                      >  1:1458480/1‑150 (MQ=255)
                                                                                                                                   ggaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgc                    <  2:46379/150‑1 (MQ=255)
                                                                                                                                    gaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCATCTTTTTcgcgcg                   >  2:1751235/1‑150 (MQ=255)
                                                                                                                                    cctgacggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCAGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTTTTCAAATTTTTcgcgcg                   >  2:1490789/7‑150 (MQ=255)
                                                                                                                                     ctgacggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcgc                  >  2:1120790/6‑150 (MQ=255)
                                                                                                                                         aGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGt              >  2:572264/1‑150 (MQ=255)
                                                                                                                                         aGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGt              >  2:1635668/1‑150 (MQ=255)
                                                                                                                                            cAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTg           >  1:1753535/1‑150 (MQ=255)
                                                                                                                                            cAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTg           <  1:765515/150‑1 (MQ=255)
                                                                                                                                              aGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTgac         >  2:47569/1‑150 (MQ=255)
                                                                                                                                                 gcCGAGCTTGCCGTCATGAAAGTGCGTCGGGAGAGTTTATATATTGGGATGCTTTCTCCTCTATAAAGCCTGAAGCACGCTGGCGAGTTCTCGATTGATCTCCTTCTGGTTAAAAACGTGTTCTACTTTTTCGCGCGCGCGTTTgacgac      <  2:1664110/150‑1 (MQ=255)
                                                                                                                                                 gcCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTgacgac      >  2:607511/1‑150 (MQ=255)
                                                                                                                                                  ccGATCTTTCCGTCAGGAACGTGCGTCGGGAGAGTTTTTGAAATGGCATGCGTTCTCCTCTATAAACCCTGCAGCAAGATGGCGTTTGATCGATTGATCACCTGCTGGTTAAAATCGTGTTCAAATTTTACGCGCGCGCGTTTGACGAcc     >  2:428930/1‑150 (MQ=255)
                                                                                                                                                  ccGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGAcc     >  1:1536892/1‑150 (MQ=255)
                                                                                                                                                    gAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCgg   >  1:113340/1‑150 (MQ=255)
                                                                                                                                                     aGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGa  <  1:1319482/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGA  >  NC_000913/2115702‑2116000

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: