Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 282,417 T→C 29.8% L47P (CTG→CCG)  yagE → 2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913282,4170TC29.8% 125.4 / 73.2 97L47P (CTG→CCG) yagE2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (37/31);  new base C (15/14);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 8.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCT  >  NC_000913/282318‑282513
                                                                                                   |                                                                                                
tctcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCt                                                                                                  >  1:1686335/1‑100 (MQ=255)
 ctcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTg                                                                                                 >  2:2000704/1‑100 (MQ=255)
  tcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTgg                                                                                                <  1:392141/100‑1 (MQ=255)
  tcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTgg                                                                                                >  1:1689628/1‑100 (MQ=255)
  tcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTgg                                                                                                >  2:2105923/1‑100 (MQ=255)
   ccaccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTggg                                                                                               >  2:1482603/1‑100 (MQ=255)
   ccaccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCgg                                                                                                >  1:40010/1‑99 (MQ=255)
    caccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGc                                                                                              >  1:540190/1‑100 (MQ=255)
     accaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCa                                                                                             >  2:951071/1‑100 (MQ=255)
      ccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAg                                                                                            <  2:505280/100‑1 (MQ=255)
       caTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGc                                                                                           <  2:1393687/100‑1 (MQ=255)
       caTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGc                                                                                           <  2:1396368/100‑1 (MQ=255)
            tttccgccgccgGCCAGCTCTAGAAGCCGGGCACCACCGCGCTTCTAGACGATTTTATAAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGAGGTgg                                                                                      <  2:1841579/97‑1 (MQ=255)
             tACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGc                                                                                     <  2:1667548/100‑1 (MQ=255)
               ccgccgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGa                                                                                   >  2:563364/1‑100 (MQ=255)
               ccgccgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGa                                                                                   >  2:428207/1‑100 (MQ=255)
                cgccgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAg                                                                                  <  2:1662251/100‑1 (MQ=255)
                 gccgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGt                                                                                 <  1:1589050/100‑1 (MQ=255)
                 gccgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGt                                                                                 <  1:1484015/100‑1 (MQ=255)
                 gccgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGt                                                                                 >  1:1090693/1‑100 (MQ=255)
                 gccgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGt                                                                                 >  2:421631/1‑100 (MQ=255)
                  ccgACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGtt                                                                                >  2:873698/1‑100 (MQ=255)
                      cGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGATCTcc                                                                            >  1:1059479/1‑100 (MQ=255)
                       ggCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTccc                                                                           <  2:1050265/100‑1 (MQ=255)
                        gacAGCTCGATATGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTcc                                                                            <  1:1358038/96‑1 (MQ=255)
                        gCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCa                                                                          >  2:1979466/1‑100 (MQ=255)
                        gCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCa                                                                          >  2:389188/1‑100 (MQ=255)
                          cAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGc                                                                        >  2:995416/1‑100 (MQ=255)
                               cGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgc                                                                   >  2:675148/1‑100 (MQ=255)
                                gATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcg                                                                  <  2:1237393/100‑1 (MQ=255)
                                 aTAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcgc                                                                 >  1:1994944/1‑100 (MQ=255)
                                  tAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGcc                                                                <  2:760800/100‑1 (MQ=255)
                                   aaGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCg                                                               <  2:571558/100‑1 (MQ=255)
                                   aaGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCg                                                               >  1:1339449/1‑100 (MQ=255)
                                   aaGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCg                                                               >  2:989665/1‑100 (MQ=255)
                                    aGCCGGGCACCGTCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGa                                                              <  1:1377636/100‑1 (MQ=255)
                                     gCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa                                                             >  2:1252705/1‑100 (MQ=255)
                                     gCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa                                                             >  2:704940/1‑100 (MQ=255)
                                     gCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa                                                             >  1:1974433/1‑100 (MQ=255)
                                      ccGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa                                                             >  1:1846112/1‑99 (MQ=255)
                                      ccGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAag                                                            >  2:398295/1‑100 (MQ=255)
                                      ccGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAag                                                            <  1:675148/100‑1 (MQ=255)
                                        gggCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAagag                                                          >  2:1944614/1‑100 (MQ=255)
                                         ggCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGc                                                         <  1:24370/100‑1 (MQ=255)
                                         ggCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGc                                                         >  1:453147/1‑100 (MQ=255)
                                          gCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCg                                                        >  2:1960675/1‑100 (MQ=255)
                                           cACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGt                                                       <  2:205827/100‑1 (MQ=255)
                                            aCCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTa                                                      <  2:2055028/100‑1 (MQ=255)
                                             ccgccgCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGt                                                       <  1:59126/98‑1 (MQ=255)
                                             ccgccgCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTaa                                                     >  2:1550440/1‑100 (MQ=255)
                                             ccgccgCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTaa                                                     >  1:728773/1‑100 (MQ=255)
                                                ccgcGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGc                                                  <  2:1038250/100‑1 (MQ=255)
                                                     cTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTg                                             >  2:1792800/1‑100 (MQ=255)
                                                      tGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGc                                            >  1:1220538/1‑100 (MQ=255)
                                                       gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGcc                                           <  1:951073/100‑1 (MQ=255)
                                                       gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGcc                                           >  2:908231/1‑100 (MQ=255)
                                                       gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGc                                            >  1:1056983/1‑99 (MQ=255)
                                                         tCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGccc                                          >  1:748867/1‑99 (MQ=255)
                                                            acgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCTGCGCCGAAGAGCGTAAAGCCATTGCCCGCtt                                      >  2:1101914/1‑100 (MQ=255)
                                                             cgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCttt                                     >  1:429547/1‑100 (MQ=255)
                                                             cgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCttt                                     <  2:1524608/100‑1 (MQ=255)
                                                              gaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCttt                                     >  1:1218216/1‑99 (MQ=255)
                                                              gaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTg                                    <  1:1387483/100‑1 (MQ=255)
                                                              gaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTg                                    >  1:946878/1‑100 (MQ=255)
                                                               aTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGc                                   >  1:1678744/1‑100 (MQ=255)
                                                               aTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGc                                   <  1:473314/100‑1 (MQ=255)
                                                                tCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCt                                  <  1:1065995/100‑1 (MQ=255)
                                                                  tGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTAt                                >  1:1716386/1‑100 (MQ=255)
                                                                   gATCAAAGCAGGCGTTAACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATc                               >  1:2112037/1‑100 (MQ=255)
                                                                    aTCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCg                              <  2:668441/100‑1 (MQ=255)
                                                                     tCAAAGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGa                             <  1:566284/100‑1 (MQ=255)
                                                                      cAAAGCATGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGTGCGTAAAGCCATTGCCCGCTTTGCTATCGAt                            <  1:282997/100‑1 (MQ=255)
                                                                      cAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGTCCGCTTTGCTATCGa                             <  2:1148638/99‑1 (MQ=255)
                                                                      cAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGAt                            <  1:1696144/100‑1 (MQ=255)
                                                                      cAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGAt                            <  2:1703329/100‑1 (MQ=255)
                                                                        aaGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCa                          <  2:86240/100‑1 (MQ=255)
                                                                        aaGCAGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCa                          <  2:40010/100‑1 (MQ=255)
                                                                        aaGCAGGCGTTGACGGCCCGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCa                          <  2:759955/100‑1 (MQ=255)
                                                                         aGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCAt                         >  2:1215940/1‑100 (MQ=255)
                                                                            tggCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTc                      <  2:594777/99‑1 (MQ=255)
                                                                            aGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTc                      <  1:49935/100‑1 (MQ=255)
                                                                            aGGCGTTGACGGCCTGTTCTTCCCGGGCAGCGGTGGCGAGTCCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTc                      <  1:1285769/100‑1 (MQ=255)
                                                                             ggCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTca                     >  1:594777/1‑99 (MQ=255)
                                                                               cGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGAt                   <  2:1177724/100‑1 (MQ=255)
                                                                                gTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGAt                   >  1:146211/1‑99 (MQ=255)
                                                                                  tGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCACCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGAtcgt                >  2:1192514/1‑100 (MQ=255)
                                                                                   gACGGCCTGTTCTTCACGGGCAGCAGTGGAGAGTTCTCCCAGCTCTGCGCCGAAGAGAGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTcttacttc               >  2:561279/1‑93 (MQ=255)
                                                                                     cGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTcgc             <  1:507917/100‑1 (MQ=255)
                                                                                        ccTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTg          <  2:94837/100‑1 (MQ=255)
                                                                                        ccTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTcgc             >  1:1358434/1‑97 (MQ=255)
                                                                                         cTGTTCTTCCCGGGCAGCGGTGGCGAGTTCTCCCAGGTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGc         <  1:986830/100‑1 (MQ=255)
                                                                                          tGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGcc        >  1:916583/1‑100 (MQ=255)
                                                                                           gTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCg       >  2:1447117/1‑100 (MQ=255)
                                                                                             tcttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCgg      >  1:2179412/1‑99 (MQ=255)
                                                                                             tcttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGt     >  1:2000649/1‑100 (MQ=255)
                                                                                               tcacTGGGCAGCGGTCGCCAGTCCTCTCAGCTCGGCGCCGAATAGCGTATATCCATTGCCCGTTTTGCTATCTATCATGTCTATCGTCGCGTGCCg       <  2:1870984/93‑1 (MQ=255)
                                                                                                tcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCt  <  1:503703/100‑1 (MQ=255)
                                                                                                tcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCt  <  2:2074876/100‑1 (MQ=255)
                                                                                                   |                                                                                                
TCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCT  >  NC_000913/282318‑282513

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: