Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,565 C→T D81D (GAC→GAT)  deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,5650CT98.0% 141.7 / ‑4.3 51D81D (GAC→GAT) deoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (29/21);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

TGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGA‑TCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGC  >  NC_000913/4617319‑4617810
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
tgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACgatgac                                                                                                                                                                                                                                                     >  1:341916/1‑250 (MQ=255)
   aatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCACCTTCCCACACGGTAACGAtgacatcg                                                                                                                                                                                                                                                 >  2:255707/1‑251 (MQ=255)
    atgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGGATTGCTACGGAAACCAACTTCCCACACGGTAACGAtgccatcg                                                                                                                                                                                                                                                 >  2:395715/1‑250 (MQ=255)
               gAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTgg                                                                                                                                                                                                                                     >  1:144480/1‑251 (MQ=255)
                  agcaagcaGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTgcc                                                                                                                                                                                                                                    >  1:400883/1‑247 (MQ=255)
                   gcaagcaGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGc                                                                                                                                                                                                                                    >  2:142028/1‑248 (MQ=255)
                         agcCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTACCCCGATTGCTCGCAAAACACTGAAAGAGCAGTGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCg                                                                                                                                                                                                                           >  1:43902/1‑251 (MQ=255)
                                    aCTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGAAGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATcc                                                                                                                                                                                                                >  1:11926/1‑250 (MQ=255)
                                         aaTTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAC‑TTGCTCTGAGTCATCAGGCCAAAACCCCGGTAGGCAATACCGCCGCTATCTGTATCTATCATCGCTTTATCCCGATTGCTAGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTCCGCTAACAAACTTCCAACACGGTAACGATGACATCGAAATAACGCTGGCAGAAACCCGTGCGGCCATCGcatcc                                                                                                                                                                                                           >  2:174079/1‑247 (MQ=255)
                                             gATGGACCCGCCCACCCTAAATGACGACGACACCGACGAAAAAGTAA‑TTGCTTTGTGTCACCAGACCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGt                                                                                                                                                                                                        <  1:299227/250‑1 (MQ=255)
                                                ggACCTGACCACCCTGAATGACGACGACCCCGACGAGAAAGCAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATATGTATCTATCCTCGCTTTATCCCGATTGCACGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCAGCATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAACTCGAGCTGCCAGGAACACTTGCGGCAATCGCCTACGGTGccg                                                                                                                                                                                                    >  1:672464/1‑249 (MQ=255)
                                                    cTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgatg                                                                                                                                                                                                 >  1:593478/1‑250 (MQ=255)
                                                     catctccccTGAATGACGACGACACCGACGA‑AAAATAAGTTGCTCTTGGTCATCAGCCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGaa                                                                                                                                                                                               <  1:39305/245‑1 (MQ=255)
                                                          accCTGAATGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGaa                                                                                                                                                                                               >  2:39305/1‑246 (MQ=255)
                                                                aaTGACGACGACACCGACGAGAAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGtggg                                                                                                                                                                                    >  1:219543/1‑248 (MQ=255)
                                                                                  tagaAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGTAAATCCGTATTGCTACGGTAACCAACTTCCAACACGGTAACGATGACCTCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGAAGTGGAGTTCCCGTAACGCGCGCt                                                                                                                                                                  >  2:509213/2‑251 (MQ=255)
                                                                                    gaAAGTAA‑TTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTg                                                                                                                                                                 >  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‑GACCTGGGGAAAGCCTGTAAAGAGGCTTTCGCGGc                                                                                                  >  1:483978/1‑249 (MQ=255)
                                                                                                                                                            gCTTTTTCCCGATTGCTCGCAAACCACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTAAGGTCACCAAATTCCCACACGGTCACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGAAATCGCCTACGGAGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACCAGCCGGCTGGTTTT‑GACCTGGTGCAAGCCAGTACAGAGGCTTgcgggggcgcggatgt                                                                                         >  2:612311/1‑237 (MQ=255)
                                                                                                                                                                 aTCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGAGCCGGTTGGGTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTActgct                                                                                   >  1:427582/1‑251 (MQ=255)
                                                                                                                                                                             cGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGatca                                                                        <  1:36444/250‑1 (MQ=255)
                                                                                                                                                                               cAAAACACTGCAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAAATTCCCACACGGTAACGATGACATCGTAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGCTGGCGGGAAACGAGCAGGTAGGGTTT‑GACCTGGTGAAAGCCTGTGAAGAGGATTGGGCGGCAGGGGATGGACTGCTGCAGATGATCATCg                                                                     >  1:370434/1‑251 (MQ=255)
                                                                                                                                                                               cAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGGTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTg                                                                                                        >  1:196722/1‑216 (MQ=255)
                                                                                                                                                                                  aacactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGaaa                                                                  <  2:548546/247‑1 (MQ=255)
                                                                                                                                                                                  aacactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGAGAGCATCGaa                                                                   >  1:622815/5‑250 (MQ=255)
                                                                                                                                                                                        gAAAGAGCAGTGCACCCCGGAAATCCGTAATGCTACGGTAACCAACTCCCCACACGGTAACGATTACATCGAAATCTGGCTGGCGGAACCCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAAGGTACTGCTTAAAGTGATCATCGAAACCGGCg                                                            <  2:597761/251‑1 (MQ=255)
                                                                                                                                                                                        gAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGGTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCgg                                                              <  1:612311/249‑1 (MQ=255)
                                                                                                                                                                                          aaGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGGATGTACTGCTGAAAGGGATCATCGAAAACGGCGa                                                           >  1:416419/1‑250 (MQ=255)
                                                                                                                                                                                          aaGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTGGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGaa                                                          <  2:483978/251‑1 (MQ=255)
                                                                                                                                                                                           agcgCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAAGAGCAGGGTGGTTTT‑TACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGGATGGACTGCTGAAAGTGATCATCAAAACCGGCGAAc                                                         >  2:294121/4‑251 (MQ=255)
                                                                                                                                                                                                            aaaTCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTg                                        <  2:244064/251‑1 (MQ=255)
                                                                                                                                                                                                                           cGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGGGAAAGCCTGTAAAGAGGCTTGCGCGGCAGGGGATGTAATGCTGAAAGTGATCATCGAAACCGGCGCACCGAAAAAGGAAGCGCTGagcccgaaagaggat                         >  2:312757/1‑237 (MQ=255)
                                                                                                                                                                                                                            ggTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTAACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGCGCACGTTGTGTTT‑GACCTGGt                                                                                                                             >  1:232235/1‑150 (MQ=255)
                                                                                                                                                                                                                                   aaCTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCCAAACCGGCGAACTTGAAGACGAAGCGCGGATCCGTAAAGCGTCTGaaaattc                 >  1:296277/1‑247 (MQ=255)
                                                                                                                                                                                                                                            cacGGTAATGATGACATCGAAATCGCGCTGGCAGAAACCAGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGc        <  1:500037/251‑1 (MQ=255)
                                                                                                                                                                                                                                             aaggTAACGATGACATCGAAATCGCGTTGGCAGAATCCCGTGCGGAAATCGCCTACGGTGCCGATGAAGGTGACGTGGTGTGCCCGTACCGCGCGCTGAGGGAGGGTAACGAGCAGGTAGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGAGCAACGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGGCTGAAATCTCCAGCAAAgag       <  2:43902/249‑3 (MQ=255)
                                                                                                                                                                                                                                                  aaCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGc  <  1:116409/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
TGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGA‑TCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGC  >  NC_000913/4617319‑4617810

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