Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 231,861 C→T 71.8% T247M (ACG→ATG)  yafD → endo/exonuclease/phosphatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,8610CT71.8% 64.8 / 27.2 43T247M (ACG→ATG) yafDendo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base T (13/17);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

ACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATATTGTTATC  >  NC_000913/231588‑232123
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                       
acGCCGTCAACTTCAGTCTGGGCGTTGATGTATATAGTAAGCCGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCTGCCCGGTCATTATGGCGGGAGATTTCAATGCCGGGAGCCGTAGAATGATGAACGCGTTATAGCGCTTTGCGCTGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcgcgcg                                                                                                                                                                                                                                                                    <  1:54056/279‑1 (MQ=255)
 cGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGGAGAAGGATGAACGCGGTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCCGGTGCGGTTTACCGCTGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCGCGATTTTGTTTTCTACCGTGGTCTGACCGTCAGTGAAGATTCTGTTCTGGTTAtgc                                                                                                                                                                                                                                                                       >  2:247108/1‑275 (MQ=255)
 cGCCGGCAAGGTCAGGGTGGGAGTGGATGTCTATAGTAAGCAGGTAGTGCCTATGGGCGATCAGATAGCGCACCAGAGCGGCCCGGTCATTATGGCGGGAGATGTCAATGCCTGGAGCCGTAGAAGGATGAACGCGGTATATCGCTGTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTAtgcgcg                                                                                                                                                                                                                                                                    <  1:181759/278‑1 (MQ=255)
    cGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTAtgcgcg                                                                                                                                                                                                                                                                    >  2:521489/1‑275 (MQ=255)
     gTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCTCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTATTGGTTACGCGCGCT‑Tc                                                                                                                                                                                                                                                               >  1:142408/1‑278 (MQ=255)
     gTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTTGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGGACTGGTTAtgcccg                                                                                                                                                                                                                                                                    >  2:215656/1‑274 (MQ=255)
     gTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTTGAGCTTCTGTACTGGTTAttcgcg                                                                                                                                                                                                                                                                    >  2:492668/1‑274 (MQ=255)
                    ggCGTGGATGTCTATAGTAAGCCGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGCAATGTCGCTGCGCCAGGTGCGGTTTTCCGATGATCAGCGCCCTCGGGCGGTTGGGCGCCCGCTCGCTTTTGTTTTCTACCGTGGTCTCAGCGTCAGTGCAGCCTCTGTTCTTGTTCTGCTCCGT‑TCCTATCCCAATccccca                                                                                                                                                                                                                                               >  1:49556/1‑275 (MQ=255)
                    ggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTCGAAGGATGAACGCGGTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCCATGATCCGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTTCCGTGGTCTGAACGACAGTGAAGCTTCCGTTCTGGTTATGCGCGCT‑TCCGATCACAATcc                                                                                                                                                                                                                                                   >  2:456860/1‑275 (MQ=255)
                      tttggagtgtattaGTAAGTGGTTATGTCCAATGGGCGATCAGATAGCTGACCACACAGGCCCGGTAATTATGGCGGGAGATGTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCGTTGAGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGc                                                                                                                                                                                                                                                 <  4:374384/263‑1 (MQ=255)
                                                    tATTGGCGATCAGATAGTGCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCt                                                                                                                                                                                                                    <  2:499823/274‑1 (MQ=255)
                                                                 gtAGCTCCGCCCAGCGGCCCGGTTATTATGGCGGGTGGTTTCAATGCCTGGGGCCGTAGAAGGATGAACCGGTTATATCGNTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTc                                                                                                                                                                                                   <  3:142408/277‑1 (MQ=255)
                                                                    gcgtaccacACGGGCCCGGTCATTATGGCGGGGGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTc                                                                                                                                                                                                   <  4:521489/271‑1 (MQ=255)
                                                                    gCTCACCACTGCGGCCCGGGCATTATTGCGGGGGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTc                                                                                                                                                                                                   <  4:446853/275‑1 (MQ=255)
                                                                          cacaGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTc                                                                                                                                                                                              <  2:250606/274‑1 (MQ=255)
                                                                                                       tttCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTg                                                                                                                                                                <  2:45185/275‑1 (MQ=255)
                                                                                                                  ggAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCCCGATTTTGTTTTCTACCTTGGTCTGAACGTCAGTGAAGCTTCCGTACTGGTTATGCGCGCT‑TCCGAGCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTa                                                                                                                                                                                                       >  2:28724/1‑225 (MQ=255)
                                                                                                                         tAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTg                                                                                                                                                                                                                                        >  4:563748/1‑185 (MQ=255)
                                                                                                                             aGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAAtt                                                                                                                                                                                                                                    <  2:563748/185‑1 (MQ=255)
                                                                                                                             aGGATGAACGCGTTATAGCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCggtgg                                                                                                                                           <  2:171994/274‑1 (MQ=255)
                                                                                                                                          tataTCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTCCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATc                                                                                                                           <  1:84360/277‑1 (MQ=255)
                                                                                                                                              tCGTTTTGCGCGGGAAATGTCGCTGCGCCAGTTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGccc                                                                                                                      <  1:81305/278‑1 (MQ=255)
                                                                                                                                               cGCTTTGTGCGGTGAATGTCTCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCcac                                                                                                                    <  1:80960/279‑1 (MQ=255)
                                                                                                                                                 cTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTAc                                                                                                                                                                                                                                              >  3:174289/1‑155 (MQ=255)
                                                                                                                                                 cTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTAc                                                                                                                                                                                                                                              <  1:174289/155‑1 (MQ=255)
                                                                                                                                                  tttGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACAt                                                                                                                  <  1:175539/278‑1 (MQ=255)
                                                                                                                                                    tggcggggAAATGTTGTTGGGCCAGGGCGGTGTTACCGGTGATAAGGGCCGTCCGGCGGTGGGTAGCCCGCTGGATTTTTTTTCCTGCCGTGGCATGAACGTCAGGGAAGCTTCTGGACTGGTTAGGCGCGCT‑TCCGGTCACAATCCGCTACTCGGTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGGATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCcaca                                                                                                                   <  4:512159/270‑1 (MQ=255)
                                                                                                                                                      cgcgGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAg                                                                                                                                                                                 >  4:309486/1‑211 (MQ=255)
                                                                                                                                                          ggAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCa                                                                                                                                                                             <  2:309486/211‑1 (MQ=255)
                                                                                                                                                                     tggcccAGGTTCGTTTTCACGTTGTGCGGCGCCGTCGGGCATTTGGTCGCCCCCTCGATTTGGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCGTCTGGACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGAAAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGc                                                                                                   <  4:456860/270‑1 (MQ=255)
                                                                                                                                                                       cgcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGcaca                                                                                                >  2:58328/1‑275 (MQ=255)
                                                                                                                                                                                   tttACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCATTCCCGGCAAGCCTGATAAATCAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGGGGCGTATGACTTATCTGCCCACATGCTTGATGTCGTGGCACAACCTTCCGAAACCCgg                                                                                >  1:1713/1‑279 (MQ=255)
                                                                                                                                                                                                       tCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGcgcg                                                                                                                                               >  4:465391/1‑196 (MQ=255)
                                                                                                                                                                                                           gCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCggtgg                                                                                                                                           <  2:465391/196‑1 (MQ=255)
                                                                                                                                                                                                           gCGTTGGGTCTCCCGCTCGATTTTTTTTTCTACCGTGGTCTGAACGTTAGTGAAGCTTCTGTACTGGTTATGGGCGCT‑TCCGAGCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcacca                                                             <  4:294379/274‑1 (MQ=255)
                                                                                                                                                                                                            cGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATAAcaccac                                                            >  2:177484/1‑274 (MQ=255)
                                                                                                                                                                                                               ttGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGGGGCg                                                                                                                                      >  1:88803/1‑197 (MQ=255)
                                                                                                                                                                                                               ttGGGCGCCCGCTCTGTTTTTTTTTCTCCCTGGGTCTGACCGTCAGGAAGGCTTCGGTACTGGTTACG‑GCGGTGTTCGATCACATTCCGTTACTCGTTGAATTCAGACCAGGCCAGCCTGTTAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGGTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGc                                                        <  4:58328/275‑1 (MQ=255)
                                                                                                                                                                                                                                              tGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGAc                                                                                                                                >  2:168590/1‑172 (MQ=255)
                                                                                                                                                                                                                                                gTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGAT‑AATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCATATGACTTATCTGCCCCCATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCCCCACCCGCCAGGGATATTCCCAAAGTCTGCCATTTGCCGaaa                      >  4:147910/1‑273 (MQ=255)
                                                                                                                                                                                                                                                  cTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTAt                                                                                                                            <  4:168590/172‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaCGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGc                  >  2:546127/1‑275 (MQ=255)
                                                                                                                                                                                                                                                           aGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataa                                                                                                                                                                                                               >  4:610399/1‑80 (MQ=255)
                                                                                                                                                                                                                                                           aGTGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATAtt        >  1:49962/1‑279 (MQ=255)
                                                                                                                                                                                                                                                             tGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCa                                                                                                                                                            >  3:25215/1‑129 (MQ=255)
                                                                                                                                                                                                                                                             tGAAGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCa                                                                                                                                                            <  1:25215/129‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaGCTTCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataaataa                                                                                                                                                                                                           <  2:610399/80‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tCTGTACTGGTTATGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGAT‑AATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATATTGTTATc  <  2:147910/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcacca                                                             >  1:80496/1‑217 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                       
ACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCT‑TCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATATTGTTATC  >  NC_000913/231588‑232123

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑