Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,616 T→C A98A (GCT→GCC)  deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,6160TC96.4% 166.7 / ‑2.3 56A98A (GCT→GCC) deoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (25/29);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

CTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGA‑TCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCC  >  NC_000913/4617371‑4617858
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
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‑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)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
CTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGA‑TCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTT‑GACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGG‑CAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGC‑TGATCCGTAAAGCGTC‑TGAAATCTCCATCAAA‑GCGGGTGCGGAC‑TTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCC  >  NC_000913/4617371‑4617858

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