Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 602,200 A→G 17.7% E81G (GAG→GGG)  pheP → phenylalanine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913602,2000AG17.7% 136.5 / 5.6 71E81G (GAG→GGG) phePphenylalanine transporter
Reads supporting (aligned to +/- strand):  ref base A (22/32);  new base G (6/8);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTGGGCTGCCGCCTTCTTTATTATCATCAACGCCGTTAACCTGGTGA  >  NC_000913/601985‑602415
                                                                                                                                                                                                                       |                                                                                                                                                                                                                       
aGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCGGAATTGCGGGGTTACGGCGTCGCCGGGCTCCTCGCTTTCCTGCTTATGCGCCAGCTTGGCGAAATGGTGGTTGGGGGGCCGGAATCCGGt                                                                                                                                                                                                            >  2:568215/1‑229 (MQ=255)
    aCTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCAGCGCTTTCCTGATTTTGCGCCAGCTTGGCGAAATGGGGGTTCAGGGGCCGGGGTCCGGTTCAt                                                                                                                                                                                                        >  2:742598/1‑229 (MQ=255)
     cTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCCTTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTCCGCGGTCGCCGGGATCATCGCTTTCCTGATTTTGGGCCCGCTGGGGGAAATGGGGGTGGGGGGGCCGGTATCCGGTTCAt                                                                                                                                                                                                        >  2:184687/1‑228 (MQ=255)
      tGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTTCTGGGATACGGCGTCCCCGGGGTCATCGCCTTCCTGATTTTGCGCCAGCTTGGCGAAATGGGGGTTGAGGAGCCGGTGCCCGGGTCCttt                                                                                                                                                                                                      >  2:521506/1‑229 (MQ=255)
       cggTCGAATCAAGGGCCGACGTTTCATCGCGGATTACATAACCGTCATTTTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTg                                                                                                                                                                                                     <  2:211180/227‑1 (MQ=255)
        cGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTTTTGCTGGGCTTCGGCGGCGCCGGGATCATCGCTTTCCTTATTATGCGCCAGCTTGGGGGACTGGGGGTGGTGGAGCCGGTGTCCCGTTCCTTTGc                                                                                                                                                                                                    >  2:610858/1‑229 (MQ=255)
        cGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGGTCATCGCCTTCCTGATTATGCGCCCGCTTGGCGCATAGGGGGTTGAGGACCCGGTCTCCGGTTCAttt                                                                                                                                                                                                      >  2:319664/1‑227 (MQ=255)
        cGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGGATTGCTGGGCTACGGCGTCGCCGGGATCCTCGCTTTCCTGCTTTTGGGCCAGCTTGGCGAAATGGTGGTTGAGGGGCCGGTATCCGGTCCAttt                                                                                                                                                                                                      >  1:333742/1‑227 (MQ=255)
             aaTCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGGGGTTGAGGAGCCGGTATCCGGTTCCTTTGCCCCCt                                                                                                                                                                                               >  1:565980/1‑229 (MQ=255)
              aTCAAGGGCCGACCCTTTCTCGGGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACtt                                                                                                                                                                                              <  1:495900/229‑1 (MQ=255)
                ctaggtcggatgCTTTATCGGGTATTAGATAACCGTCTGATTCAGCTGATTGCGGTGGGTGGCGCAATGGGTTCTGGGCTGTTTCTTGGCATTGGCCCGGCGATTCGGGTGGCGGGTCCGGCTGTAGTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCGTGATTATGCGCCATCTTGGCGAAATTGTGGTTGAGGGGCCGGTATCCGGTTCATTTGCCCCCTTTg                                                                                                                                                                                            <  2:600008/218‑1 (MQ=255)
                cAAGAGCCGACTCTTCATCGCGGATTATATAACCGTTATATTCAACTGATTGCGTTGGGTGGCGGAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTg                                                                                                                                                                                            <  2:210738/229‑1 (MQ=255)
                     gccgagctgtCATCGCGGTTTACATAACCTTCATATTCAGCTGATTGCGTTGGGTGGCGCAATTGGTGCTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCtat                                                                                                                                                                                       <  2:42039/220‑1 (MQ=255)
                     gCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCtat                                                                                                                                                                                       <  1:95514/229‑1 (MQ=255)
                                             tagccGTCATTTTCAACGGGTTGGGTTGGGTGGCGGAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGc                                                                                                                                                               <  1:332893/226‑1 (MQ=255)
                                                     atatTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGGGCCGCCTTGGCGAAAAGGGGGTTGAGGGGCCGGTACCCGGGGACTTTGCCCCCTTTGCCTTATAATACGTGGGGCCGTTTGGGGGGTTctctcc                                                                                                                                                       >  2:710754/1‑224 (MQ=255)
                                                     atatTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCCGCTTGGCGCAATGGGGGTTGAGGGGCCGGTTTCCGGTTCATTTTGCCCCTTTGCCTGTAAATACTTGGGCGCGTTTGCGGGCTTCctctc                                                                                                                                                       >  2:4380/1‑229 (MQ=255)
                                                            agtggtttggTTGGGGGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGGGATTTAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGGGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTctggctgg                                                                                                                                                <  2:237328/220‑1 (MQ=255)
                                                             cTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGTGCCGGTATCCGGTTCATTTGCCCGCTTTGCCTATAAATACTGGGGACCGGTTCCGGGCTTCCTCTCTGCTTGGa                                                                                                                                               >  1:678023/1‑229 (MQ=255)
                                                             cTGATTGCGTTGGGCAGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTCCTGGGCTACGCCATCACCGGGATCATCGCTTTCCTGACTATGCGCCAGCTTGTCGACATGGGGGTTTCGGTGCGGGTCTCCGGTTAGTTTGCCCATTTTTCCTTTCAACTCTGCGGACCGTTTTCGGGCTTCctctctcggcgg                                                                                                                                                >  2:103520/1‑222 (MQ=255)
                                                                aTTGGGTGGGGTGGCGCAATTGGTTCTTGTTCTTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCATCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAact                                                                                                                                            <  1:231276/229‑1 (MQ=255)
                                                                 ttGCGTTGGGTGGCGCAATTGTTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTAGTGCTGGGCTACGGCGTCGGCGGGATCACCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGTGACCGTTTGCGGGCTTCCTCTCTGGCTGGAacta                                                                                                                                           <  1:739655/229‑1 (MQ=255)
                                                                       tGGGTGGCGCAAGTGGTTCTGGTCTGTTTCTTGGCATTGGCCCGGAGATTCAGATGGCGGGTCCGGCTGTATTGCGGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGt                                                                                                                                     <  1:371686/229‑1 (MQ=255)
                                                                         ggTGGCTCAATTGGTTCTGGTCGGTTTCTTGGCATTGGCCCGGCGATTGAGATGGGGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGt                                                                                                                                     <  1:319664/227‑1 (MQ=255)
                                                                              cgcAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGGTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGGAATGTTc                                                                                                                              >  1:681357/1‑229 (MQ=255)
                                                                                 attgggTTCTGGTGTGTTTTTTGGCTTTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTg                                                                                                                           <  2:663844/225‑1 (MQ=255)
                                                                                 aattgggatgtgTCTGTTTCTTGGCATTTGGCCGGCGTATCAGATGGCGGGTCCGGCTGTAATGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTg                                                                                                                           <  1:205253/218‑1 (MQ=255)
                                                                                     gggttgggTCTTTTTTTTGGCCTTGGCCCGGCGATTGAGATGGCGGGTCCGGCGGTATTGCGTGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCtg                                                                                                                        <  2:678022/222‑1 (MQ=255)
                                                                                            gTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCTGGCTGGGGGGaa                                                                                                                 >  1:531400/1‑228 (MQ=255)
                                                                                               tGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATgga                                                                                                             >  1:367518/1‑228 (MQ=255)
                                                                                                 tttCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGATTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGGGGGAAGGGccg                                                                                                           >  1:725453/1‑227 (MQ=255)
                                                                                                   tCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCtat                                                                                                                                                                                       <  2:100150/151‑1 (MQ=255)
                                                                                                   tCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCtat                                                                                                                                                                                       >  1:100150/1‑151 (MQ=255)
                                                                                                      tGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGGGCTGGGGGGAAGGGCAGAGCtt                                                                                                      >  1:21316/1‑228 (MQ=255)
                                                                                                      tGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGGGCTGGCGGGAATGGCCGCGCTg                                                                                                      >  1:181925/1‑229 (MQ=255)
                                                                                                      tGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACCTTGCCTCTCAATAATGGGGACCGTTTGCGGGCTTCGTCCTCGGCGGGAACTACTGGGTAATGTTGCGGCTGGGGGGAATGGccgggccg                                                                                                      >  2:356875/1‑222 (MQ=255)
                                                                                                      tGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCCTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGGTTCCGGGCTTCCTCCCTGGCTTGAACTTCTGGGCAATGTTCCTGTTGGTGGGAAGGGCAGGGCTg                                                                                                      >  2:484664/1‑229 (MQ=255)
                                                                                                      tGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGCCGCCGGGATCATCGCTTTCCTGTTTATGCGCCAGCTTGGCGAAATGGGGGTTGATGAGCCGGTTTCAGGTTCATTTGCCCACTTTGCCTCTAACAAAGGGGGCCCGCTTGCGGGGTTTATCTCCGGGTGGAATTACTGGGGCATGATCAGGCGGGGGGGCATGGCaaggcta                                                                                                      >  2:24131/1‑223 (MQ=255)
                                                                                                          cttGGCCCGGCGATTCGGATGGCGGTTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTATTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCg                                                                                                  <  2:333742/228‑1 (MQ=255)
                                                                                                             ggCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCAt                                                                                                                                                                                                        <  2:506610/124‑1 (MQ=255)
                                                                                                             ggCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCAt                                                                                                                                                                                                        >  1:506610/1‑124 (MQ=255)
                                                                                                             ggCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGGGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGGGGGAATGGCAGAGCGGACCCCtc                                                                                               >  2:660333/1‑228 (MQ=255)
                                                                                                               cccgggagttCAGGTGGCGGGTCCGGTTGTTTTGGTGGGGTAGGGCGTCGGCGGGATCATCGTTTTCCAGATTATGCGGCAGCTTGGCGAAATGGTGGTTGAGGGGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCg                                                                                             <  2:234090/221‑1 (MQ=255)
                                                                                                               cccGGCGGTTCAGTTGGCGGGTCCGGTTGTAGTGGCGGGGTTGGGCTTCGCCGGGTTCATCGCTTTCCTGGTTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGGGCCGGTATCCGGTTCATTTGCCCACGTTGCCTATAAATACTGGGGGCCGTTTGCGGGCGTACTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCg                                                                                             <  2:42628/229‑1 (MQ=255)
                                                                                                                 cggcggTTCGGATGGCGGGTCCGGCTGTTTTGCTGGGTTACGGCGTCGCCGGGATCATCTCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCggg                                                                                           <  2:725452/229‑1 (MQ=255)
                                                                                                                 ccggcaTTCTGTTGGGGGGGCCGGCTGTATTGCTGGGCTCCGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCTATTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCggg                                                                                           <  2:608870/224‑1 (MQ=255)
                                                                                                                 cGGTGATTCAGGTGGCGGGTCCGGGTGTATTGGTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACGTTGCCTATAAATTCTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCggg                                                                                           <  1:719485/229‑1 (MQ=255)
                                                                                                                 cGGCTATTCAGAGGGCGGGTCCGGTTTTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATGATGCGCCAGCTTGGCGAAATGGTGGTTGAGGGGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCggg                                                                                           <  1:742599/229‑1 (MQ=255)
                                                                                                                 cGGCGATTCAGGTGGCGGGTCGGGCTGTATTGCTGGGCTCTGGCGGCGCCGGGATCATCGCTTGCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCggg                                                                                           <  1:98711/229‑1 (MQ=255)
                                                                                                                 cGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCggg                                                                                           <  1:440836/229‑1 (MQ=255)
                                                                                                                      tttCAGGTGGCGGTTCCGGTTTTTTTGTGGGGCTCCGGTGGCGCCGGTATCCGCGCTTTCCGGTTTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCt                                                                                      <  2:307186/228‑1 (MQ=255)
                                                                                                                      aTTCAGATGGCGGGTCCGGTTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGGTGAGGATCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCt                                                                                      <  1:367047/229‑1 (MQ=255)
                                                                                                                           ggttgCGGGTCCGGTTTTTTTGGTGGGGTGCGGCGTCGCCGGGATTATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATg                                                                                 <  2:340970/225‑1 (MQ=255)
                                                                                                                           gATGGCGGGTTCGGCTGTTTTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATg                                                                                 <  1:224806/229‑1 (MQ=255)
                                                                                                                              ggCTGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCa                                                                               <  1:610859/228‑1 (MQ=255)
                                                                                                                                 gggTCCGGGTTTTTTGGCGGGGGACGGCGTTGCCGGGATATGCGCTTTCCTGATTATCCGCCAGCTTGGCAAAATGGTGGTTGAGGAGCCGGTAGCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTAc                                                                           <  2:531400/229‑1 (MQ=255)
                                                                                                                                 gggTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGCACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTTTATGCAGTaa                                                                           >  1:341529/1‑228 (MQ=255)
                                                                                                                                                      cTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAATTCCTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTGTTCCTGGTGGTAATGGCGGACCTGACCCCTGCGGGACTTTATATGCAGACCGGGTTCCCGGATGTCCCAAc                                                       >  2:47576/1‑228 (MQ=255)
                                                                                                                                                      cTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGTAATGGCAGAGCTGAACGCTGCGGGCGTCTATCTTCTGTACTGTTTCCCGGCTGTTCCAAc                                                       >  1:423587/1‑228 (MQ=255)
                                                                                                                                                       tACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACg                                                      <  1:47576/228‑1 (MQ=255)
                                                                                                                                                              tCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTGGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGCACCTATATGCAGTATTGGTTCCCGGATGTTCCAACGAGGATTAgg                                             >  1:529310/1‑229 (MQ=255)
                                                                                                                                                                 ccgggtttcttgCTTTCCTGTTTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGGGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTggg                                            <  2:21316/217‑1 (MQ=255)
                                                                                                                                                                   gggATCCTCGCTTTCCTGATTTAGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTGGGCTg                                         <  1:448596/229‑1 (MQ=255)
                                                                                                                                                                   gggATCATCGCTTTCCGTATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGTGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTGGGCTg                                         <  1:660334/229‑1 (MQ=255)
                                                                                                                                                                        catcGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTGGGCTgccgc                                     <  1:134152/228‑1 (MQ=255)
                                                                                                                                                                        cagtgCTTTCCTGATTATGGGCCAGGTTGGCGAAATGGTGGTTGAGGGGCCGGTACCCGGTTCATTGGCCAACTTTGCCTATAAATACTGGGGACCGTTGGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTGGGCTgccgc                                     <  1:356875/224‑1 (MQ=255)
                                                                                                                                                                                 ccTGTTTTTGCCCCAGTTTGGCGAAATGTTGTTTGGGGGGCCGTTTTCCGGTTCATTTGCCCACTTTGCCTTTAAATACTGGGGACCGTTTGCGGGTTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTGGGCTGCCGCCTTCTttatt                           <  2:341529/229‑1 (MQ=255)
                                                                                                                                                                                        aTGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGGGGATTTGGGCTGCCCCCTTCTTTTTTatcatc                     >  2:594238/1‑228 (MQ=255)
                                                                                                                                                                                                   tGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTGGGCTGCCGCCTTCTTTATTATCATCAACGCCGTTAAc         >  1:376549/1‑229 (MQ=255)
                                                                                                                                                                                                      cGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTGGGCTGCCGCCTTCTTTATTATCATCAACGCCGTTAACCTg      <  1:594239/229‑1 (MQ=255)
                                                                                                                                                                                                            ggtggtTGAGGGGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTGGGCTGCCGCCTTCTTTATTATCATCAACGCCGTTAACCTGGTGa  <  1:484664/227‑1 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                       
AGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGTTGAGGAGCCGGTATCCGGTTCATTTGCCCACTTTGCCTATAAATACTGGGGACCGTTTGCGGGCTTCCTCTCTGGCTGGAACTACTGGGTAATGTTCGTGCTGGTGGGAATGGCAGAGCTGACCGCTGCGGGCATCTATATGCAGTACTGGTTCCCGGATGTTCCAACGTGGATTTGGGCTGCCGCCTTCTTTATTATCATCAACGCCGTTAACCTGGTGA  >  NC_000913/601985‑602415

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