Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 284,003 A→G 15.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,0030AG15.3% 160.8 / 2.4 73E268G (GAG→GGG) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base G (8/4);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGGGAAGACGCGGGCAAGGTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGTGCT  >  NC_000913/283781‑284240
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
tCGGGGAAGACGCGGGCAAGGTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGGGTGGCTGGGGGAcccccgcccggcg                                                                                                                                                                                                                                >  1:422805/1‑234 (MQ=255)
           gcgGGCAAGGTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGGGTCAGTTCCTCGGCACGGCGGGCACCCCGCCGGGGGGCGCGGAGGCGCTGGGGCCGGGGCCGCCGCCCTCCCCCCCGGCGCCGCCCGGGCCGGCGGTGTGGCCGGCGATGCGCCCCCACGCGGCGCCCCCgg                                                                                                                                                                                                                      >  2:546884/1‑238 (MQ=255)
              ggCAAGGTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGGGATCGCCCCCCCGTCGGCGCCCGCGGTcc                                                                                                                                                                                                                   >  1:1314796/1‑237 (MQ=255)
                   ggTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCCGGGGGTCGCGGAGGCGCTGGGGCTGGCGCCGCCGCACTCCGCGCCGGCGCCGTCCCGGCCGGCGGGGGGGGTGGGGAGCCCCCGCCCGTGGCCGCGCGCGGGGAgcgcgc                                                                                                                                                                                                             >  2:1396126/1‑239 (MQ=255)
                           ccATCGGCGCGCGTTTCGCCAACCACGAACTCCCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGc                                                                                                                                                                                                     <  1:483927/239‑1 (MQ=255)
                               cGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGGGGTCGCGGGGGCGCTGGGTCTGGCGCTGCCGCCCTCCGCGCCGGGGCCGTCCGGGCCGGCGGGGGGGCGGGAGAGTGGCCCCCCGGCCGCGGCCGGGGTCAGCCGGCCGGGTTGCcgcg                                                                                                                                                                                                 >  2:171358/1‑239 (MQ=255)
                               cGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGGGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACCCCGCGCTGGCGCCGTCCGGGCCGGCGGGGTGGCTGGGGGTCGCCCGCCCCGCCGCGCGCCCGGGCAGCCGGCCGGGTAGCcccg                                                                                                                                                                                                 >  2:1208479/1‑239 (MQ=255)
                                   gcgcgTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCAt                                                                                                                                                                                             <  1:171358/239‑1 (MQ=255)
                                      cgTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCCCCTCGCCGGTGGGCCCGGCGGGGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCCGGCGGTGTGGCTGGCGGTCTCCCCCCCCTCGCCCCGCGCGGTCCGCGGGCTGGGTAGCcgcgccgccc                                                                                                                                                                                           >  2:773026/1‑232 (MQ=255)
                                         tacGCCACCCGCGAGCTCTCCCGGCAGGAGGCCGCCGGACGGGGCGGTCGGGCCGGGCCCTGGCGGGGGGGGGGGGGTCGGTTCCTGGGCACGGCGGGCACCTCGCAGGGGGTCGGGGGGGCGCTGGGTCTGGGGCTGCCGCGCCCCGCGCGGGCGCCGTCCGGGCGGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATcaccac                                                                                                                                                                                       <  2:709824/237‑1 (MQ=255)
                                                     gAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGGGCTGGGTAGCCGCGGCATCACCCCGCGGGGGATcc                                                                                                                                                                            >  1:624608/1‑238 (MQ=255)
                                                       aCTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCctct                                                                                                                                                                         <  1:968830/239‑1 (MQ=255)
                                                            cccTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCCGGGGGGCGCGGAGGCGCTGGGGCTGGGGCTGCCGCACTCCCCGCTGGCGCCGTCCGGGCCGGCGGGGGGGGTGGGGGTAGGCCGCCCCGCGGCCGGCGCGGGCAGCGGGCGGGGCGGCCGCGGCCGCACCCCCCGGGGGATCCCCCCCggt                                                                                                                                                                    >  2:1039107/1‑237 (MQ=255)
                                                                  agggggCCGCCGAGGGGGGCTGCCGCGCGCGCGCCTGGCCGGGCGGCGGGTGTCGGTCCCCCGGCACGGCGGGCCCCCCGCAGGTGGTCGCGGGGGCGCGGGGTCTGGCGCTGCCGCGCTCCGCGCTGGCGCCGTCCGGGCAGGCTGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGc                                                                                                                                                               <  2:103982/238‑1 (MQ=255)
                                                                       gccgccGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGGGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGGGGGGCTGGAGGTCGCCCCCCAGGCGGGGCGCGCGGGCGGGCAGCCGGATAGCCGCGGCCCCCCCCCGCCGGGGATCCTCTCCGCTAAACCCCTCgc                                                                                                                                                         >  2:1236664/1‑238 (MQ=255)
                                                                            cgACCTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAAcg                                                                                                                                                    <  1:1236664/239‑1 (MQ=255)
                                                                                        tCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGGGATCCCCCCCCCGTCGGCGCGCGCGGGCAGCCGGCTGGGGAGCCGCGGCATCCCCACCCGGGGTATCCCCTCCCCTAAAGCCCTCCAAAAcccggaggggaac                                                                                                                                         >  2:798870/1‑229 (MQ=255)
                                                                                        tCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGGGGGGCTGGGGAACGCCCGCCCGGCGGCCCGCGCGGGCCGCGGGCTGGGGAGCCGCGGCCTCCCCCCCCGGGGGATCCTCTCCCATAAAAACATCGAAAAcccgggggggctc                                                                                                                                         >  2:1317052/1‑229 (MQ=255)
                                                                                              ctgcccctCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCGGGTGGTCGCGGGGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACgcgg                                                                                                                                  <  1:1317052/239‑1 (MQ=255)
                                                                                                         cGGGCGGCGGGTGTCAGTTCCCCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCGGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCggcggc                                                                                                                       <  1:1208479/239‑1 (MQ=255)
                                                                                                               gggggTGTCGGTTCCTCGGCACGGCGGGCCCCTCGCAGGTGGTCGGGGAGGCGCGGGGTCGGGCGCGGCCGCACTCCGGGCTGGCGCGGTCCGGGGAGGCGGGGGGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTccac                                                                                                                  <  1:1496020/236‑1 (MQ=255)
                                                                                                               gggggTGTCGGTTCCCGGGCACGGCGGCCCCCTCGCGGGGGGCCGCGGGGGCGCTGGGTCTGGCGCGGCCCCACCCCGCGCTGGCGCCGCCGGGGCAGGCGGTGTGGCTGGAGACCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTccac                                                                                                                  <  2:624608/236‑1 (MQ=255)
                                                                                                               gcggGTGTCAGTTCCTCGGCCCGGCGGGCCCCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCCCTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTccac                                                                                                                  <  1:79901/238‑1 (MQ=255)
                                                                                                                     gtCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGGGATCGCCCGCCCGTCGGCGCGCGCGGGCAGCGAGCTGGGTAGCCGCGGCATCCCCCCGCGGGGGATCCCCTCCCCGAAAACCCCCGAAAACACGAAGGGGGTCCCCGCGGGGGTCGGCGGGGCCACCAActt                                                                                                            >  2:1162061/1‑238 (MQ=255)
                                                                                                                                                        gtCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGTTATCCTCTCCGATAAAGCCATCGAAAACGAGATGGTGATCCACGCGGCGGTCGGCGGCTCCACCAATTTACTGCTGCGCATTCCGGCCCTCGCCCACGGGGCcgg                                                                        >  2:1019922/1‑239 (MQ=255)
                                                                                                                                                          cgcgGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGcggcgg                                                                         >  2:185828/1‑236 (MQ=255)
                                                                                                                                                          cgcgGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATAGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCCATAAAGCCATCGAAAACGCGGTGGTGATCCACGCGGCGTTCGGCGGCTCCCCCCATTTACTGCTGCACCTTCCGGCCATCGCCCCCGcggcgg                                                                         >  2:1468621/1‑236 (MQ=255)
                                                                                                                                                              ggggCGCTGGGTCGGGCCCGGCCGCACCCCGCGCTGGCGCCGTCCGGGCGGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCa                                                                   <  1:773026/236‑1 (MQ=255)
                                                                                                                                                              gAGGCGCTGGGTCTGGCGCTGCCGCCCTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGACCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCa                                                                   <  1:798870/238‑1 (MQ=255)
                                                                                                                                                              gAGGCGCTGGGTCTGGCGCGGCCGCACTCGCGGCTGGCGCCGGCCGGGCAGGCGGTGTGGTTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCa                                                                   <  1:1529696/238‑1 (MQ=255)
                                                                                                                                                                ggCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAAATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGGAAGCCGCGGCCTCACCACGCGGGGAATCCTCTCCGATAAAGCCATCGCAAACGCGGTGGTGGTCCACGCGGCGTTCGGCGGCGCCACCAATTTACAGCTGCACATTCCGGCCATCGCCCCACCGGCGGGCTTCagag                                                                >  2:1165975/1‑236 (MQ=255)
                                                                                                                                                                       ggtggggCGCCGCCGCCCCCCGCGCTGGCGCGGCCCGGGCGGGCGTGGTGGCGGGAGCTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCgg                                                         <  2:1382210/233‑1 (MQ=255)
                                                                                                                                                                           tGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCCCCACGCGGGATATCCTCTCCGATAAAGCCATCGCAAACGCGATGGTGATTCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACCCCGCGGGCAGCACGATCCCGGACGtt                                                    >  2:176331/1‑239 (MQ=255)
                                                                                                                                                                           tGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGtt                                                    >  1:1382210/1‑239 (MQ=255)
                                                                                                                                                                                 tGCCGCACTCCGGGCGGGCGCCGCCCGGGCGGCGGGGGGGGCGGGGCACGGCCCGCCCGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCa                                               <  2:999992/238‑1 (MQ=255)
                                                                                                                                                                                   ccGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGTGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCCGGATAGCCGCGGCATCACCACGCCGGATATCCTCTCCGATAAAACCATCGAAAACGCCATGGTGAACCCCGCGGCGTTCGGCGGCCCCACCAATTTACTGCTGCCCACTCCGTGCATCGCCCACGCCGCGGGCCTCATGAACCCCGGACGtgaccaatg                                             >  2:397003/1‑231 (MQ=255)
                                                                                                                                                                                   ccGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGCGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGCTAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGCCCCTCCCCCCCCCGGCGGGGCGCACCATCCCGGACGGTGAGcccac                                            >  2:1162077/1‑235 (MQ=255)
                                                                                                                                                                                   ccGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCCTCGCCCACGCGGCGGGCTGCCCGATCCCGGCCGTTTGGCAccg                                            >  2:1167506/1‑237 (MQ=255)
                                                                                                                                                                                   ccGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCCCGATCCCGGCCGTTGGGCACTg                                            >  2:1093720/1‑239 (MQ=255)
                                                                                                                                                                                   ccGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGGCGGTGAGCACtc                                            >  2:1181412/1‑238 (MQ=255)
                                                                                                                                                                                     gCACTCCGCGCTGGGGCCGTCCGGGCAGGCGGTGTGGCGGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGa                                          <  1:1468621/239‑1 (MQ=255)
                                                                                                                                                                                        cTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgc                                       >  2:1102353/1‑239 (MQ=255)
                                                                                                                                                                                           cgcgCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCCGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCCCCCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                                     >  2:1407248/1‑238 (MQ=255)
                                                                                                                                                                                           cgcgCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                                     >  2:352500/1‑238 (MQ=255)
                                                                                                                                                                                            gcgtgggCGCCGTCCGGGCGGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                                   <  1:1162061/234‑1 (MQ=255)
                                                                                                                                                                                            gcgcTGGCGCCGCCGGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                                   <  1:1158316/239‑1 (MQ=255)
                                                                                                                                                                                             cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCCCCAATTTACTGCTGCACATTCCGGCCATCGCCCCCCCGGCGGGCTGCCCGATCCCGGACGTTGAGCACTGGACGCGCATc                                  >  2:546239/1‑239 (MQ=255)
                                                                                                                                                                                             cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATc                                  >  1:1264359/1‑239 (MQ=255)
                                                                                                                                                                                             cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATc                                  >  2:664826/1‑239 (MQ=255)
                                                                                                                                                                                             cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAAATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGCGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGGTAAAGCCATCGGAAACGCGAACGGGATCCACGCGGCGATCGGCGGCCCCACCAATTTACCGCTGCACCTTCCGGCCATCCCCCACGCAGCCGGCTGCACGATCCCGGCAGGTGAGCACtcggagggcatc                                  >  2:1385438/1‑228 (MQ=255)
                                                                                                                                                                                             cgcTGGCGCCGGCCGGGCAGGGGGTGGGGCGGGGGACCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                                   <  1:763975/238‑1 (MQ=255)
                                                                                                                                                                                              gcgggCGCCGGCCGGGCGGGGGGGGTGGCGGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATc                                  <  1:841098/235‑1 (MQ=255)
                                                                                                                                                                                                  gcgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCCCGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCCCATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACCCGCATCAACCg                             >  1:309916/1‑239 (MQ=255)
                                                                                                                                                                                                  gcgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGCGCACTGGAcgccgcacaaccg                             >  2:1083339/1‑229 (MQ=255)
                                                                                                                                                                                                  gcgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCg                             >  1:1205885/1‑239 (MQ=255)
                                                                                                                                                                                                   cgcCGTCCGGGCAGGCGGTGTGGCTGGAGGTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGCCGTTGAGCCCTGGACGCGCATCACCCGt                            >  2:135664/1‑239 (MQ=255)
                                                                                                                                                                                                   cgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGGGGGCTGCACGCTCCCCGGCGGTGAGCACTGGAAGCGCATCAACCGt                            >  2:856101/1‑239 (MQ=255)
                                                                                                                                                                                                   cgcCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGt                            >  1:1173194/1‑239 (MQ=255)
                                                                                                                                                                                                    gcggTCGGGGCGGGCGGTGTGGCGGGAGTACGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                           <  1:1067960/236‑1 (MQ=255)
                                                                                                                                                                                                     ccGTCCGGGGGGGGGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                           <  1:1298599/238‑1 (MQ=255)
                                                                                                                                                                                                     ccGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGGGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGCCGCGCCTCAACCCTaa                          >  1:601513/1‑239 (MQ=255)
                                                                                                                                                                                                         ccGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGt                       <  1:185828/238‑1 (MQ=255)
                                                                                                                                                                                                          gggggcggCGGTGGGGCGGGGGACCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc                     <  1:712152/233‑1 (MQ=255)
                                                                                                                                                                                                          ggggcgggcGGTGGGGCGGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc                     <  1:546239/238‑1 (MQ=255)
                                                                                                                                                                                                          cgggcgggGGGGGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc                     <  1:1019922/239‑1 (MQ=255)
                                                                                                                                                                                                          cGGGCAGGCGGTGTGGCTGGGGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc                     <  1:550614/239‑1 (MQ=255)
                                                                                                                                                                                                          cGGGCAGGCGGGGGGGCTGGAGATCGCCCGCCAGGCGGCGCGCGCGGGCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc                     <  2:826936/239‑1 (MQ=255)
                                                                                                                                                                                                             gCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCcg                   <  1:1093720/238‑1 (MQ=255)
                                                                                                                                                                                                                  cGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                           >  1:513208/1‑225 (MQ=255)
                                                                                                                                                                                                                      gtgGCTGCAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGa         <  1:839297/239‑1 (MQ=255)
                                                                                                                                                                                                                      gtgGCGGGAGAGCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCGCCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGa         <  1:397003/239‑1 (MQ=255)
                                                                                                                                                                                                                       tgGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGGGCTGGATAGCCGCCGCATCACCCCGCGGGATATCCTCTCCCATAAAGCCATCGAAAACGCGGTGGTGATTCACGCGGCGTTCCGCGGACCCACCAATTTACCGCTGCACATTCCGGCCATCGCCCACCCGCCGGGCGGCACGAGCCCGCACGGTGGGCACTCGACGAGCAACAACCGGCAAGGTCCGCGCCCGGTGa         >  2:1231404/1‑238 (MQ=255)
                                                                                                                                                                                                                       tgGCTGGAGATCCCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGa         <  1:1113490/238‑1 (MQ=255)
                                                                                                                                                                                                                              agaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGtgct  <  1:1102353/238‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
TCGGGGAAGACGCGGGCAAGGTGCAGACCATCGGCGCGCGTTTCGCCAACCACGAACTCTCCCTGCAGGAGGCCGCCGAACTGGGCTGTCGCGCCTGCGCCTCGCCGGGCGGCGGGTGTCAGTTCCTCGGCACGGCGGGCACCTCGCAGGTGGTCGCGGAGGCGCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGTGCT  >  NC_000913/283781‑284240

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