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

ACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATATTGTTATC  >  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)
     gTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCTCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTATTGGTTACGCGCGCTTc                                                                                                                                                                                                                                                               >  1:142408/1‑278 (MQ=255)
     gTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTTGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGGACTGGTTAtgcccg                                                                                                                                                                                                                                                                    >  2:215656/1‑274 (MQ=255)
     gTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTTGAGCTTCTGTACTGGTTAttcgcg                                                                                                                                                                                                                                                                    >  2:492668/1‑274 (MQ=255)
                    ggCGTGGATGTCTATAGTAAGCCGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGCAATGTCGCTGCGCCAGGTGCGGTTTTCCGATGATCAGCGCCCTCGGGCGGTTGGGCGCCCGCTCGCTTTTGTTTTCTACCGTGGTCTCAGCGTCAGTGCAGCCTCTGTTCTTGTTCTGCTCCGTTCCTATCCCAATccccca                                                                                                                                                                                                                                               >  1:49556/1‑275 (MQ=255)
                    ggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTCGAAGGATGAACGCGGTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCCATGATCCGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTTCCGTGGTCTGAACGACAGTGAAGCTTCCGTTCTGGTTATGCGCGCTTCCGATCACAATcc                                                                                                                                                                                                                                                   >  2:456860/1‑275 (MQ=255)
                      tttggagtgtattaGTAAGTGGTTATGTCCAATGGGCGATCAGATAGCTGACCACACAGGCCCGGTAATTATGGCGGGAGATGTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCGTTGAGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGc                                                                                                                                                                                                                                                 <  4:374384/263‑1 (MQ=255)
                                                    tATTGGCGATCAGATAGTGCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCt                                                                                                                                                                                                                    <  2:499823/274‑1 (MQ=255)
                                                                 gtAGCTCCGCCCAGCGGCCCGGTTATTATGGCGGGTGGTTTCAATGCCTGGGGCCGTAGAAGGATGAACCGGTTATATCGNTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTc                                                                                                                                                                                                   <  3:142408/277‑1 (MQ=255)
                                                                    gcgtaccacACGGGCCCGGTCATTATGGCGGGGGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTc                                                                                                                                                                                                   <  4:521489/271‑1 (MQ=255)
                                                                    gCTCACCACTGCGGCCCGGGCATTATTGCGGGGGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTc                                                                                                                                                                                                   <  4:446853/275‑1 (MQ=255)
                                                                          cacaGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTc                                                                                                                                                                                              <  2:250606/274‑1 (MQ=255)
                                                                                                       tttCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTg                                                                                                                                                                <  2:45185/275‑1 (MQ=255)
                                                                                                                  ggAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCCCGATTTTGTTTTCTACCTTGGTCTGAACGTCAGTGAAGCTTCCGTACTGGTTATGCGCGCTTCCGAGCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTa                                                                                                                                                                                                       >  2:28724/1‑225 (MQ=255)
                                                                                                                         tAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTg                                                                                                                                                                                                                                        >  4:563748/1‑185 (MQ=255)
                                                                                                                             aGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAAtt                                                                                                                                                                                                                                    <  2:563748/185‑1 (MQ=255)
                                                                                                                             aGGATGAACGCGTTATAGCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCggtgg                                                                                                                                           <  2:171994/274‑1 (MQ=255)
                                                                                                                                          tataTCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTCCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATc                                                                                                                           <  1:84360/277‑1 (MQ=255)
                                                                                                                                              tCGTTTTGCGCGGGAAATGTCGCTGCGCCAGTTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGccc                                                                                                                      <  1:81305/278‑1 (MQ=255)
                                                                                                                                               cGCTTTGTGCGGTGAATGTCTCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCcac                                                                                                                    <  1:80960/279‑1 (MQ=255)
                                                                                                                                                 cTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTAc                                                                                                                                                                                                                                              >  3:174289/1‑155 (MQ=255)
                                                                                                                                                 cTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTAc                                                                                                                                                                                                                                              <  1:174289/155‑1 (MQ=255)
                                                                                                                                                  tttGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACAt                                                                                                                  <  1:175539/278‑1 (MQ=255)
                                                                                                                                                    tggcggggAAATGTTGTTGGGCCAGGGCGGTGTTACCGGTGATAAGGGCCGTCCGGCGGTGGGTAGCCCGCTGGATTTTTTTTCCTGCCGTGGCATGAACGTCAGGGAAGCTTCTGGACTGGTTAGGCGCGCTTCCGGTCACAATCCGCTACTCGGTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGGATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCcaca                                                                                                                   <  4:512159/270‑1 (MQ=255)
                                                                                                                                                      cgcgGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAg                                                                                                                                                                                 >  4:309486/1‑211 (MQ=255)
                                                                                                                                                          ggAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCa                                                                                                                                                                             <  2:309486/211‑1 (MQ=255)
                                                                                                                                                                     tggcccAGGTTCGTTTTCACGTTGTGCGGCGCCGTCGGGCATTTGGTCGCCCCCTCGATTTGGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCGTCTGGACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGAAAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGc                                                                                                   <  4:456860/270‑1 (MQ=255)
                                                                                                                                                                       cgcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGcaca                                                                                                >  2:58328/1‑275 (MQ=255)
                                                                                                                                                                                   tttACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCATTCCCGGCAAGCCTGATAAATCAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGGGGCGTATGACTTATCTGCCCACATGCTTGATGTCGTGGCACAACCTTCCGAAACCCgg                                                                                >  1:1713/1‑279 (MQ=255)
                                                                                                                                                                                                       tCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGcgcg                                                                                                                                               >  4:465391/1‑196 (MQ=255)
                                                                                                                                                                                                           gCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCggtgg                                                                                                                                           <  2:465391/196‑1 (MQ=255)
                                                                                                                                                                                                           gCGTTGGGTCTCCCGCTCGATTTTTTTTTCTACCGTGGTCTGAACGTTAGTGAAGCTTCTGTACTGGTTATGGGCGCTTCCGAGCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcacca                                                             <  4:294379/274‑1 (MQ=255)
                                                                                                                                                                                                            cGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATAAcaccac                                                            >  2:177484/1‑274 (MQ=255)
                                                                                                                                                                                                               ttGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGGGGCg                                                                                                                                      >  1:88803/1‑197 (MQ=255)
                                                                                                                                                                                                               ttGGGCGCCCGCTCTGTTTTTTTTTCTCCCTGGGTCTGACCGTCAGGAAGGCTTCGGTACTGGTTACGGCGGTGTTCGATCACATTCCGTTACTCGTTGAATTCAGACCAGGCCAGCCTGTTAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGGTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGc                                                        <  4:58328/275‑1 (MQ=255)
                                                                                                                                                                                                                                              tGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGAc                                                                                                                                >  2:168590/1‑172 (MQ=255)
                                                                                                                                                                                                                                                gTCTGAACGTCAGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCATATGACTTATCTGCCCCCATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCCCCACCCGCCAGGGATATTCCCAAAGTCTGCCATTTGCCGaaa                      >  4:147910/1‑273 (MQ=255)
                                                                                                                                                                                                                                                  cTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTAt                                                                                                                            <  4:168590/172‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaCGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGc                  >  2:546127/1‑275 (MQ=255)
                                                                                                                                                                                                                                                           aGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataa                                                                                                                                                                                                               >  4:610399/1‑80 (MQ=255)
                                                                                                                                                                                                                                                           aGTGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATAtt        >  1:49962/1‑279 (MQ=255)
                                                                                                                                                                                                                                                             tGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCa                                                                                                                                                            >  3:25215/1‑129 (MQ=255)
                                                                                                                                                                                                                                                             tGAAGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCa                                                                                                                                                            <  1:25215/129‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaGCTTCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataaataa                                                                                                                                                                                                           <  2:610399/80‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tCTGTACTGGTTATGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATATTGTTATc  <  2:147910/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATcacca                                                             >  1:80496/1‑217 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                       
ACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGATGTCGTGGCACAAGCTGCCGAAGCCCGGCAACTGAAAAATATCACCACCCGCCAGGGATATGCCGAAAGTCTGCCATTTGCCGATAACGCATTTGATATTGTTATC  >  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: