Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,631 T→G V103V (GTT→GTG)  deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,6310TG92.7% 137.7 / 1.9 55V103V (GTT→GTG) deoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (23/28);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

GAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATG  >  NC_000913/4617401‑4617877
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                            
tagaAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGTAAATCCGTATTGCTACGGTAACCAACTTCCAACACGGTAACGATGACCTCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGAAGTGGAGTTCCCGTAACGCGCGCt                                                                                                                                                                                                                                          >  2:509213/2‑251 (MQ=255)
  gaAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTg                                                                                                                                                                                                                                         >  1:244064/1‑250 (MQ=255)
     agtaattgctctgagtCATCAGGCCAAAACCCCGGTCGGCAATACCACCGCTATCTGTCTCTATCCTCGCTTTATCCCGATAGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGCAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGGTGACATCGAAATCGCGCGGGCAGCAACCCGTGCGGCAATCTCCTACGGTGACGATGAAGTATACGTGGTGGTCCCGAACCGCGCGCTGacgg                                                                                                                                                                                                                                     >  2:219610/15‑248 (MQ=255)
           tgctcggtgtCATCAGCCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGt                                                                                                                                                                                                                                <  1:281345/244‑1 (MQ=255)
            gctcTGTGTCACCAGGCCAAAACCCCGGTCGTCAATTCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTa                                                                                                                                                                                                                               <  1:255707/247‑1 (MQ=255)
               cTGTGACATCAGCCGAAACCCCCGGTAGACAATACCGCCGCTAGCTGTATCTATCCTCGCCTTCTCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCGACATACCACACGGTTACGATGCAATCGAAATCGCGCTGACAGAAACCCGTACGGCAATCGCATTCGGGGCCGATGACGTTGATGTGGTGTtaccgtaccgaccgctaccggcggggtacg                                                                                                                                                                                                                            >  2:124547/1‑230 (MQ=255)
                   gtCATCAGGCCAAAACCCCGGTCGGCAAATCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAg                                                                                                                                                                                                                       <  1:630271/251‑1 (MQ=255)
                   gtCAACAGGCAAAACCCCCGTTCGGCATAACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAg                                                                                                                                                                                                                       <  1:224296/251‑1 (MQ=255)
                      atcaGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGa                                                                                                                                                                                                                                        >  1:380394/1‑231 (MQ=255)
                        cggtccacacccccGGTCGGAAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggttgg                                                                                                                                                                                                                  <  2:146221/239‑1 (MQ=255)
                                acccccgcggggaaacaCCGCCGCGAGCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GAcc                                                                                                                                                                                                         <  2:341916/235‑1 (MQ=255)
                                          ggCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGc                                                                                                                                                                                               >  1:103645/1‑251 (MQ=255)
                                           gaaaTACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGcc                                                                                                                                                                                              <  1:211407/249‑1 (MQ=255)
                                                 ccgccgCTAACTGTATCTAACCTCGCCTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCGCCCCGTAAATCCGTACTGCTACGTGAACCAACTTCCCACGCGGTAACCATGACATCGAAATCGCGCTGGCAGAAACCCGCGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGAGGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTaaa                                                                                                                                                                                        <  2:275221/251‑1 (MQ=255)
                                                  cgccgcTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGtca                                                                                                                                                                                         >  1:557682/1‑247 (MQ=255)
                                                    ccgcTATCGGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAagag                                                                                                                                                                                     <  2:219543/251‑1 (MQ=255)
                                                             gTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGTAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGAGCCGATGCAGTTGACGTGGTGTTCACGTCCAGCGCGCTGATGTCGGGAAACGCGCAGGGGGGGTTT‑TACCTGGTGAAAGACTGTCAAGAGGCTTGcgcg                                                                                                                                                                            >  2:483050/1‑251 (MQ=255)
                                                             gTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGcgcg                                                                                                                                                                            <  1:395715/251‑1 (MQ=255)
                                                               atctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGGGCAGGTTGGTTTTTGCCCGGGT‑AAAGCCTGAAAAGAgggtttgcgggg                                                                                                                                                                          >  2:232742/1‑241 (MQ=255)
                                                                 ctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGGGAAAGCCTGTAAAGAGGCTTTCGCGG‑c                                                                                                                                                                         >  1:483978/1‑249 (MQ=255)
                                                                         gCTTTTTCCCGATTGCTCGCAAACCACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTAAGGTCACCAAATTCCCACACGGTCACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGAAATCGCCTACGGAGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACCAGCCGGCTGGTTTT‑GACCTGGTGCAAGCCAGTACAGAGGCTTgcgggggcgcggatgt                                                                                                                                                                 >  2:612311/1‑237 (MQ=255)
                                                                              aTCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGAGCCGGTTGGGTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTActgct                                                                                                                                                          >  1:427582/1‑251 (MQ=255)
                                                                                          cGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGatca                                                                                                                                               <  1:36444/250‑1 (MQ=255)
                                                                                            cAAAACACTGCAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAAATTCCCACACGGTAACGATGACATCGTAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGCTGGCGGGAAACGAGCAGGTAGGGTTT‑GACCTGGTGAAAGCCTGTGAAGAGGATTGGGCGG‑CAGGGGATGGACTGCTGCAGATGATCATCg                                                                                                                                            >  1:370434/1‑251 (MQ=255)
                                                                                            cAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGGTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTg                                                                                                                                                                                >  1:196722/1‑216 (MQ=255)
                                                                                               aacactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGaaa                                                                                                                                         <  2:548546/247‑1 (MQ=255)
                                                                                               aacactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGAGAGCATCGaa                                                                                                                                          >  1:622815/5‑250 (MQ=255)
                                                                                                     gAAAGAGCAGTGCACCCCGGAAATCCGTAATGCTACGGTAACCAACTCCCCACACGGTAACGATTACATCGAAATCTGGCTGGCGGAACCCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAAGGTACTGCTTAAAGTGATCATCGAAACCGGCg                                                                                                                                   <  2:597761/251‑1 (MQ=255)
                                                                                                     gAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGGTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCgg                                                                                                                                     <  1:612311/249‑1 (MQ=255)
                                                                                                       aaGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGGATGTACTGCTGAAAGGGATCATCGAAAACGGCGa                                                                                                                                  >  1:416419/1‑250 (MQ=255)
                                                                                                       aaGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTGGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGaa                                                                                                                                 <  2:483978/251‑1 (MQ=255)
                                                                                                        agcgCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAAGAGCAGGGTGGTTTT‑TACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGGATGGACTGCTGAAAGTGATCATCAAAACCGGCGAAc                                                                                                                                >  2:294121/4‑251 (MQ=255)
                                                                                                                         aaaTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑Tg                                                                                                              <  2:244064/251‑1 (MQ=255)
                                                                                                                                        cGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGGGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGGGGATGTAATGCTGAAAGTGATCATCGAAACCGGCGCACCGAAAAAGGAAGCGC‑TGagcccgaaagaggat                                                                                               >  2:312757/1‑237 (MQ=255)
                                                                                                                                         ggTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTAACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGCGCACGTTGTGTTT‑GACCTGGt                                                                                                                                                                                                     >  1:232235/1‑150 (MQ=255)
                                                                                                                                                aaCTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCCAAACCGGCGAACTTGAAGACGAAGCGC‑GGATCCGTAAAGCGTC‑TGaaaattc                                                                                      >  1:296277/1‑247 (MQ=255)
                                                                                                                                                         cacGGTAATGATGACATCGAAATCGCGCTGGCAGAAACCAGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑Gc                                                                            <  1:500037/251‑1 (MQ=255)
                                                                                                                                                          aaggTAACGATGACATCGAAATCGCGTTGGCAGAATCCCGTGCGGAAATCGCCTACGGTGCCGATGAAGGTGACGTGGTGTGCCCGTACCGCGCGCTGAGGGAGGGTAACGAGCAGGTAGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGAGCAACGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGGC‑TGAAATCTCCAGCAAA‑gag                                                                           <  2:43902/249‑3 (MQ=255)
                                                                                                                                                               atcgaggactgcgaaatcgCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGAAGTGGTTTTCCCGTACCGCGCGCTGAGGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGACTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGCCGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTg                                                                       <  2:92871/235‑1 (MQ=255)
                                                                                                                                                               aaCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGc                                                                      <  1:116409/251‑1 (MQ=255)
                                                                                                                                                                  gatgacatcgaaatcGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGCTGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAAAGACTTGAGCGGACGGGGAAT‑TACTTCTGAAAGGGATGATCGAAACCGGCGAACTGACAGGCGCAGAGC‑CGAGCCGGAAAACGGC‑TGAAATCTCCATCCAA‑GACGGTGCggg                                                                   >  2:578196/4‑250 (MQ=255)
                                                                                                                                                                   atgacatcgaaatcGCGCTGGCAGAAGCCCGGGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAc                                                                  <  2:367427/249‑1 (MQ=255)
                                                                                                                                                                    tgacatcgaaatcGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTAATGCTGAAAGGGGTCATCGAAACCGGCGAACTGAAAGACGCAGCGC‑TGGTCCGTAAAGCGTC‑TGAAATCTCCATCCAA‑ACGGGGGCgggct                                                                 >  1:597892/2‑248 (MQ=255)
                                                                                                                                                                                            aaCCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCCTCAAAACCTCTACCGGTaaa                                         >  1:248797/1‑250 (MQ=255)
                                                                                                                                                                                               ccGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCTCTACCGGTAAAGTgg                                     <  1:623554/251‑1 (MQ=255)
                                                                                                                                                                                                  tGCGGCAATCGCCTACGGGGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGATAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATTTACTTCTGAAAGAGATCATCCACACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAACGTC‑TTAAATCTCCATCCAA‑GCGTGTGCGGAC‑TTCCTCAAAACCTCTAcccgtcaagatgctt                                  >  1:193063/1‑238 (MQ=255)
                                                                                                                                                                                                   gCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCTCTACCGGTAAAGTGGCtgt                                 <  1:635825/251‑1 (MQ=255)
                                                                                                                                                                                                   gCGGCAATCGCCTACAGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGATGAGGGCGGGTAACGAGCAGGGTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCTCTACCGGTAAAGTGGCtgt                                 <  2:315042/251‑1 (MQ=255)
                                                                                                                                                                                                    cGGCAATCGCCTACGGTGCCGATGAAGATGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCCGGG‑CCACGAATTTACTGCTGGAAGGGGTCATTGAAATCGGCGAAATCAAAAACGAAACGCGTG‑TCCAGAAAGCGTCGTGAAA‑CTCCATCAAAAGCGGG‑GCGGACCTCCAACAAAACCTTTAC‑GGTTAAGTgggtggg                                >  2:523746/1‑246 (MQ=255)
                                                                                                                                                                                                       cAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCCCTACCGGTAAAGTGGCTGTGAAc                             <  2:144480/251‑1 (MQ=255)
                                                                                                                                                                                                       cAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTT‑‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAAcg                            <  1:142028/251‑1 (MQ=255)
                                                                                                                                                                                                                aCGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGcc                     <  1:604532/250‑1 (MQ=255)
                                                                                                                                                                                                                 cGGTGCCGATGAAGGTTAAGTGGTGTTCCCGTACCGGGCTCTTATGGCGGGGAACGAGGAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCCGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGTGTC‑TGAAATCTCCATCAAA‑GCGGGGGCGGAC‑GTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGcc                     <  2:416419/249‑1 (MQ=255)
                                                                                                                                                                                                                              gTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGACTGTAAAGAGGCTTGCTCGG‑CAGCTAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCAACAAAACCTCTACTGGTAAAGTGGCTGTGAACTCGACGCCGTAAAGCTCGCGcat      <  2:427582/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aCGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCTCTACCGGTAAAGTGGATGTGAACGCGACGCCGGAAAGCGCGCGCATCatg  <  1:312757/251‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                            
GAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATG  >  NC_000913/4617401‑4617877

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