Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 284,003 A→G 23.3% E268G (GAG→GGG)  yagF → CP4‑6 prophage, dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913284,0030AG23.3% 100.4 / 4.0 56E268G (GAG→GGG) yagFCP4‑6 prophage, dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base A (18/23);  new base G (12/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGCCGCCGACCGTCGGGGAAGACGCGGGCAAGGTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGTG  >  NC_000913/283766‑284238
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
cGCTGCCGCCGACCGTCGGGGAAGACGCGGGCAAGGTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGga                                                                                                                                                                                                                                             <  2:1567641/238‑1 (MQ=255)
          gACCGTCGGGGAAGACGCGGGCAAGGTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGACGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCCCCCCGGGGGTCGCCGTGCGGCTCGGTCCGGCGCCCCCCGCCCCCCCGGTCGCGCCGGTCGGGCACGCGCGGTGGCCGGACGTcgcccgc                                                                                                                                                                                                                                   >  1:740885/1‑238 (MQ=255)
                      aGACGCGGGCAAGGTGCAGCCCATCGGCGGGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCGGGGCGGGGGGTGTCAGTTCCTCGGCACGGCGGGCCCCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGgcgcgcg                                                                                                                                                                                                                      <  2:554099/239‑1 (MQ=255)
                        ggcgggcatgctggcgcccaccgggCGCGGTTGCGCCACCCCCGCACCCCCCCGGCGGGGGCCCGCGGACCTGGGCGGCGCCGCCTGCGCCTCGCCGGGGGGCGGGTGTCAGTTCCTGGGCCGGGGGGGCACCTCGCAGGTGGTGGCGGAGGCGCTGGGTCGGGCGCGGCCGCACTCCGCGCTGGCTCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGcgcgcgcg                                                                                                                                                                                                                    <  1:685941/215‑1 (MQ=255)
                             ggCAAGGTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGGGTCAGTTCCTCGGCACGGCGGGCCCCCCGCCGGGGGGCTCGGAGGCGCTGGGGCTGGCGCGGCGGCCCCCCGCGCCGGCGCCGGCCGGGCCGGCGGGGGGGCTGGTGGTACCCCGCCCGTGGCCGCCCGCGGGCa                                                                                                                                                                                                                >  1:424812/1‑238 (MQ=255)
                             ggCAAGGTGCAGACCACCGGCGCGCGTTTCGCCAACCACGAACTCCCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCCCCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCa                                                                                                                                                                                                                <  2:1705084/238‑1 (MQ=255)
                                  ggTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCCCGCTGGCGCCGCCCGGGCAGGCGGGGGGGCGGGGGATCGCCCGCCAGTCCGCGCGCCCGGTCAGCggcg                                                                                                                                                                                                          >  1:1657219/1‑236 (MQ=255)
                                        gACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCCCGCTGGCGCCCGCCGGGCAGGCGGGGTGGCTGGACATCGCCCCCCAGTCGGCGCGCGGCGCCAGCGCGCCgggga                                                                                                                                                                                                    >  1:729783/1‑236 (MQ=255)
                                              cGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGGGGCGGGGGATCGCCCGCCCGACGGCGCGCGCGGGCAGCGGGCTGGATAAGcgcg                                                                                                                                                                                              >  2:1216859/1‑239 (MQ=255)
                                              cGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGGGTGGCTGGGGATCGCCCGCCAGTCGGCGCGCGCGGTCCGCGAGGTGGGTAGCcccg                                                                                                                                                                                              >  2:879426/1‑239 (MQ=255)
                                              cGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGGCCGGGCAGGCGGGGTGGCTGGAGATCGCCCCCCAGTCGGCGCGCGCGGTCAGCGAGCTGGGTAGCcgcg                                                                                                                                                                                              >  2:1237858/1‑239 (MQ=255)
                                              cGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCCGGCGCCGGCCGGGCAGGCGGTGGGGCTGGAGATCCCCCGCCCGCCGGCCGCCGCCGGCAGCCCGCTGGGTAGcccgg                                                                                                                                                                                              >  2:405306/1‑236 (MQ=255)
                                                   cgcgTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGGGGTGGTGGGGGTCGCCCGCCCGTCGGCGCGCGCGGGCCGCGGGGCGGGTAGGCGCgggct                                                                                                                                                                                          >  1:816552/1‑235 (MQ=255)
                                                               aCCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGGGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGGGGGGCTGGGGGACCCCCCCCCGTCGGCGCGCGCGGGCCGCCAGCTGGGAAGCCGCgggctccccccccggggg                                                                                                                                                                             >  1:496076/1‑223 (MQ=255)
                                                                ccACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGGGTGGCTGGGGATCGCCCGCCCGGCCGCGCGCGCGGTCAGCGAGCTGGGTAGCCGCGGCATCACCCCGCCGGgta                                                                                                                                                                            >  2:837739/1‑239 (MQ=255)
                                                                          tcCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGGGATCGCCCGCCCGTCGGCGCGCGCGGTCAGCGGGCCGGATAGCCCCGGCCTCACCACGCGGGGTATCCTCCCCgg                                                                                                                                                                  >  2:815902/1‑238 (MQ=255)
                                                                           cccTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGAt                                                                                                                                                                 <  2:161103/239‑1 (MQ=255)
                                                                              tGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGGGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGc                                                                                                                                                                                                  <  2:452098/203‑1 (MQ=255)
                                                                                cAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCCCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCCGGCAGGCGGGGTGGCTGGGGATCGCCCGCCCGTCGGCGCGCGCGGGCCGCGGGCTGGATAGCCGCGGGCCCACCCCCCGGGGTATCCCCTCCCAGAAAg                                                                                                                                                             >  1:1325493/1‑238 (MQ=255)
                                                                                  ggaggCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCCCAGGGGGTCGCGGAGGCGCGGGGTCTCGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGGAGGCGGCGTGGCTGGGGATCGCCCGCCCGTCGGCGCGCGCGGTCAGCGAGCTGGGGCGCCGCGGCATCACCCCCCGGGAGATACTCTCCGACACAGcc                                                                                                                                                           >  1:219494/1‑238 (MQ=255)
                                                                                      gccgccGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCCGGTGGTCGCGGAGGCGCTGGGTCTGGCCCTGCCGCCCTCCCCGCCGGCGCCGTCCGGGCCGGCGGGGGGGCCGGGGGTCGCCCCCCCGGCCGCGCGCGCGGGCCGCGCGCCGGGTAGCCGCGGCACCCCCCCCCCGGATATCCTCTcccgcaacgacaccca                                                                                                                                                      >  1:69253/1‑225 (MQ=255)
                                                                                          ccgAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGGGGTCGCGGGGGCGCTGGGTCTGGGGCTGCCGCACTCCGCGCTGGCGCCGGCCGGGCAGGCGGGGGGGGTGGGGGACGCCCGCCCGTCGGCGCGCGCGGGGAGCGGGCCGGGGAGCCGCGGGACCACCCCGCGGGGGCTCCTCCCCGGTAAAACGCTCGGAAAc                                                                                                                                                  >  1:846197/1‑239 (MQ=255)
                                                                                                                              gggggTGGCGGTTCCCCGGCGGGGCGGGCCCCTCGGGGGTGGGCGGGGAGGGGCGGGGGCGGGCGCTGCCGCCCTCCGCGCTGGCGCCGCCCGGGCGGGCGGTGGGGCTGGAGACCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTccac                                                                                                               <  1:837739/236‑1 (MQ=255)
                                                                                                                                   tgtCAGTTCCTCGGCCCGGCGGGCCCCTCGCGGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCCCTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAAtt                                                                                                          <  2:496076/238‑1 (MQ=255)
                                                                                                                                      cggTTCCGCGGCACGGCGGGCACCTCGCAGGTGGTGGCGGAGGCGCGGGGTCTGGCGCTGCCGCACCCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTAct                                                                                                      <  2:1543192/237‑1 (MQ=255)
                                                                                                                                           ccTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGGGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACtgctgc                                                                                                 <  2:1065334/239‑1 (MQ=255)
                                                                                                                                                 cACGGCGGGCCCCTCGCAGGTGGTCGCGGGGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTc                                                                                                                  <  2:1667307/216‑1 (MQ=255)
                                                                                                                                                      cggGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGGGGCTGGAGATCGCCCCCCAGTCGGCCCGCGCGGTCAGCGGGCTGGATAGCCGCGGCCTCACCACCCGGGATATCCTCTCCGGTAAAGCCATCGAAAACGCGATGGGTATCCACGCGGCCTTCGGGGGCCCCCCCCATTTACTCCTGCCCATTCCCGc                                                                                       >  1:172721/1‑238 (MQ=255)
                                                                                                                                                      cggGCACCTCGCAGGGGGGCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGGGTGGCTGGGGATCGCCCGCCAGTCCGCGCGCGCGGTCAGGGGGCTGGATAGCCGCCGGCGCACCACGCCGGGAATCCCCTCCCATAAAGCCATCCAAAACGGCCGGGGGATGCCCGCCGCGCTCGGCGCATCCACCCAATTTCCGCTTCACACTCCCGcc                                                                                      >  1:1628791/1‑239 (MQ=255)
                                                                                                                                                               cGCAGGTGGGCGCGGGGGCGCTGGGCCGGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGGGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCAc                                                                             <  2:1478266/239‑1 (MQ=255)
                                                                                                                                                                   ggtgtggccgGGGGCGGTGGGTGGGGGGCGGCCGCCCCCCGCGCAGGCCCCGCCCGGGCAGGGGGTGGGGCTGGGGACCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGTGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCCCCAATTTACTGCTGCACATTCCGGCCATCGCCCAcgcg                                                                          <  1:1216859/230‑1 (MQ=255)
                                                                                                                                                                   gggggtCGGGGGGGCGCTGGGTCGGGCGCTGCCGCCCTCCGCGCTGGCGCCGCCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCAcgcg                                                                          <  2:1525068/235‑1 (MQ=255)
                                                                                                                                                                       gtCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGGGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCCCCACCAATTTACTGCTGCACATTCCAGCCCTCGCCCACGcggcgg                                                                      >  1:1659234/1‑238 (MQ=255)
                                                                                                                                                                       gtCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCCCGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGGGATCCACGCGGGGTTCGGCGGCGCCACCAATTTACTGCTGACCATTCCGGCCATCGCCCACGcggcgg                                                                      >  1:1708315/1‑238 (MQ=255)
                                                                                                                                                                       gtCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGcggcgg                                                                      >  1:1651976/1‑238 (MQ=255)
                                                                                                                                                                       gtCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGcggcgg                                                                      >  2:1388543/1‑238 (MQ=255)
                                                                                                                                                                             ggggggcggggTCGGCCGCGCCCGCACGCCGCGCGGGCGCCGCCCGGGCGGGCGGTGGGGCTGGAGATCCCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCa                                                                <  1:879426/230‑1 (MQ=255)
                                                                                                                                                                             ggggCGCTGGGTCGGGCGCGGCCGCGCGCCGCGCGGGCGCCGTCCGGGCAGGCGGTGTGGCGGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCa                                                                <  2:1420679/236‑1 (MQ=255)
                                                                                                                                                                                gcgcTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATACCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCAATCGCAAACGCGATGGTGAACCACGCGGCCGTCGGCGGCTCCACCCATTTACGGCTGCACATTCCGGCCATCGCCCCCGCCACGGGCCGCAccg                                                             >  1:373520/1‑236 (MQ=255)
                                                                                                                                                                                    ggggTCGGGCGCGGCCGCACTCCGCGCGGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATccc                                                         <  2:24811/237‑1 (MQ=255)
                                                                                                                                                                                      ggTCTGGCGCTGCCGCACCCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCgg                                                       <  2:1614233/238‑1 (MQ=255)
                                                                                                                                                                                           ggCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGGGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                                 >  2:672866/1‑239 (MQ=255)
                                                                                                                                                                                                tGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCCGCGGGCGGCCCGCTCCCGGGCGTTGGGcc                                             >  1:1251469/1‑237 (MQ=255)
                                                                                                                                                                                                tGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCCCGCTGGCTAGCCGCGGCATCACCACGCGGCATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGGTCCGCGGATCCACCCATTTCTTGCTGCAGATTCCCGCCATCCCCCACGCGGCGGGCCGCACGATCCCGGACGCTGAGcc                                             >  1:24811/1‑237 (MQ=255)
                                                                                                                                                                                                    gCACCCCGCGCTGGCGCCTTCCGGGCGGGCGGTGTGGCTGGGGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGa                                        <  2:200062/239‑1 (MQ=255)
                                                                                                                                                                                                      aCTCCGCGCTGGCGCCGCCCGGGCGGGGGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                                      <  2:1456588/239‑1 (MQ=255)
                                                                                                                                                                                                      aCTCCGCGCTGGCGCCGCCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                                      <  2:816552/239‑1 (MQ=255)
                                                                                                                                                                                                           aaaagggacagggcggggCGGGGGGTGGGGCGGGAGACCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                                 <  1:1569567/224‑1 (MQ=255)
                                                                                                                                                                                                            cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCCGGCTGCACGAACCCGGACGTTGAGCACTGGGCGCGCATc                                >  1:1665881/1‑239 (MQ=255)
                                                                                                                                                                                                                         cgggcgggcgggGTGGCGGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc                   <  2:1050506/239‑1 (MQ=255)
                                                                                                                                                                                                                         cGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc                   <  2:846197/239‑1 (MQ=255)
                                                                                                                                                                                                                               ggCGGTGTGGTTGGAGACCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGt              <  2:92520/238‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                       >  2:1569567/1‑227 (MQ=255)
                                                                                                                                                                                                                                 cGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTg           >  1:132916/1‑239 (MQ=255)
                                                                                                                                                                                                                                      tgGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCGACCGTTAACTGCCGCGTCTGGTGa       >  1:1410166/1‑238 (MQ=255)
                                                                                                                                                                                                                                      tgGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGa       >  2:436722/1‑238 (MQ=255)
                                                                                                                                                                                                                                           ggAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGtg  <  2:1251469/238‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
CGCTGCCGCCGACCGTCGGGGAAGACGCGGGCAAGGTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGTG  >  NC_000913/283766‑284238

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