Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,400,230 T→C 16.6% R3G (AGA→GGA)  ogt ← O‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,400,2300TC16.6% 117.8 / 2.8 65R3G (AGA→GGA) ogtO‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (25/27);  new base C (7/5);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGG  >  NC_000913/1400010‑1400447
                                                                                                                                                                                                                              |                                                                                                                                                                                                                          
atGCAAAGATTACCGGGAAAATATTACCGAAGGTTGTCCTTTAAACCGCCTGGATTGGTTGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGGACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCAt                                                                                                                                                                                                                      <  2:222069/227‑1 (MQ=255)
   cTAAGATTACCGGCAAAATATTCACGAAACTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCCCTGTACTCTTCCCATTCAACCGCCCGCCGGCGAAATTGCTCATCGCCAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGAAACCCCCGCCTCCGt                                                                                                                                                                                                                  >  1:408646/1‑228 (MQ=255)
      aGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTTAAGTAAATCCAGCATCCGtttt                                                                                                                                                                                                              >  2:73953/1‑226 (MQ=255)
        aTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCCGGCGAAATTGCTCATCGCCAATCACCCACAGTGGCCCCAGTGGCGTGGCAATTTTTTCTTCAATAAAACCCACCATCCGctctct                                                                                                                                                                                                            >  1:137324/1‑228 (MQ=255)
        aTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGCAATTGCTCATCGCAACTCACTCACAGTGGACCCAGTGGCGTGGCAACTTTGTCTTTAAGAAATCTCACCATCCCtttcc                                                                                                                                                                                                             >  2:404377/1‑224 (MQ=255)
        aTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTTTTCAAGTAAtcttcgaatccctttcct                                                                                                                                                                                                             >  2:393636/1‑213 (MQ=255)
        aTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGGGGACCCAGTGGCCGGGCCATTTTTTTTTAAAGTAACTCCCAACACCCGTtctc                                                                                                                                                                                                             >  1:629089/1‑228 (MQ=255)
          ttcgggCAAAATTTTCCCGAATTTGTTCGTTTAACCCGCCTGGATTGGGGGCAGAAATGCGCTCATAGCCTTCTTTGCCATAATGTATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTctctct                                                                                                                                                                                                          <  1:393636/224‑1 (MQ=255)
           aCCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCCCAGTGGACCCAGTGGCGGGGCAATTTTTTCTTCCAGTAATCTCACCATcccgttcctctt                                                                                                                                                                                                          >  2:632979/1‑218 (MQ=255)
            ccGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATGTTTTTCTTCAGTAATTCCAGCACCCCTTGCCTCttc                                                                                                                                                                                                        >  2:251931/1‑227 (MQ=255)
                 aaaaTATTCACGAAGTTTTTCGTTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACa                                                                                                                                                                                                   <  2:92537/228‑1 (MQ=255)
                     tatTCCCGAAGGTTAGCGCTTAAACCGCCTTGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTTCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACg                                                                                                                                                                                               <  2:435284/228‑1 (MQ=255)
                          aCGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAAATTCAGCAACCGTCTTTCCTTAAGCAAAACGTgcg                                                                                                                                                                                           >  2:157026/1‑225 (MQ=255)
                                   tccGCTAAACCCGCCTGGATTGGTGGCAGAAATGCCCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCg                                                                                                                                                                                 <  2:137324/226‑1 (MQ=255)
                                        ttAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGACATTTTTTTTTTCAGGAATCTCCGCCTCCGGTCTCTCTTTAGACACACGTGGGTAAATAccccgtcc                                                                                                                                                                             >  2:136301/1‑222 (MQ=255)
                                                ccTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAATTCACCCACAGTGGACCCAGTGGCGTGGAAATTTTTTTTTCAAGTAATCTCAGCATCCCTCCACCCTTAAAACAACCGTGGGCAAATACCCCGAccatccgtgag                                                                                                                                                                    >  1:73948/1‑219 (MQ=255)
                                                  tGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCCCTTAAGACCAACGTGGGTAAATACCCGCACATATCgggaca                                                                                                                                                                   >  1:419121/1‑223 (MQ=255)
                                                     aTTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCCGTGGCGTGGCAATTTTTTCCTTAAGTAAACTCACCACCCCTTCTCCCTTCAGCAAACCCTGGGTAAATAACCGGACCCACCGCCAAAccc                                                                                                                                                                >  1:9433/1‑226 (MQ=255)
                                                     aTTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGa                                                                                                                                                               <  2:398716/228‑1 (MQ=255)
                                                        ggtggCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGtt                                                                                                                                                                                                                >  2:383899/1‑176 (MQ=255)
                                                         gtggCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTtc                                                                                                                                                                                                               <  1:383899/176‑1 (MQ=255)
                                                          tggCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCCCAGTGGCCCCAGTGGCGTGGCAATTTTTTCTTCAAGGAATCCCCGCATCCGTTCTCTCTTTAGACAAACGTGGGGTAATCCCCGAGCCAAACGGCAGAccccacaac                                                                                                                                                          >  1:609830/1‑221 (MQ=255)
                                                             ccgAAATGCCCTCATAGCCCTCTTTGCGATACAGGCTGCCCAGCAGCTGCCCCATCCGTTCCCTGTACTCTCCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCTTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATAc                                                                                                                                                       <  1:157026/226‑1 (MQ=255)
                                                                  aTGCGCTCATAGCCTTTGTTGCGATAATGGATGTCCAGCAGCTGCCCCATGCGTTCGCTGGACTCTTCCCATACAACCGCCCGGAGGCGAAATTGCTCATCGCAAATCGTCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACg                                                                                                                                                   <  2:394208/227‑1 (MQ=255)
                                                                          aTAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACCTTGGTAAATACCCGGACCTATCGGCAGACCGCCACGATACCCACCCCGCAACa                                                                                                                                          >  2:60525/1‑228 (MQ=255)
                                                                          aTAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAAACTTAGCATCCGTTCTCTCTTAAGACAAACATGGGTGAATACCCGGACCTATCTGCAGACCGACAAGATACCACGCCAGCAACa                                                                                                                                          >  2:489102/1‑228 (MQ=255)
                                                                            tcccTTTTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATc                                                                                                                                        <  1:632979/226‑1 (MQ=255)
                                                                                  cttTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTTTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGGGGGTAAAAACCCGGACCCATCGGCCGACCGACAAGATACCACGCCAGCAACATCagcagc                                                                                                                                  >  2:552114/1‑228 (MQ=255)
                                                                                  cttTGCGATAATGGATGCCCAGCAGCTGAACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGGGGACATTTTTTCTTCAAGGACTCTTATGATCCGTATTGCTCTAATGAAAACGTGGTGGAAAACCCGGACCATACGGCAGAACCTCAAGACAACCCGCCGCCAAAAAAagcagc                                                                                                                                  >  1:500812/1‑228 (MQ=255)
                                                                                       cGATAATGGATGTCCAGCAGCTGCACCATGAGTTCGCGGTACTCTTCCAATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATac                                                                                                                             <  2:408646/228‑1 (MQ=255)
                                                                                           aagggATGCCCACCCGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCGACCCCCCGCAGGCGAAATTGCTCATCGCAAATCCCCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATCCCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATaccacc                                                                                                                         <  1:552114/225‑1 (MQ=255)
                                                                                              ggATGCCCAGCCACTGCACCATGCGTTCGCTTTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCCCCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCaaa                                                                                                                      <  2:141914/228‑1 (MQ=255)
                                                                                              ggATGCCCAGCACCTTCCCCATGCGTTCGATGTCCTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCCCCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCaaa                                                                                                                      <  2:609830/228‑1 (MQ=255)
                                                                                               gATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCCGCATCCGTTTTCTCTTAAAACAAACGTGGGGAAATACCCGGACCTATCGGCATGCCGACCAAATACAACGCCAGCAACATCAGCAGCACTACCACCaaaa                                                                                                                     >  2:354633/1‑228 (MQ=255)
                                                                                                          gctgccctatcgcttCGCTTTCCCCTCCCCATTAAACCGCCCGCAGGCTAAATTGCTCATCGCAAATCCCCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTAAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGa                                                                                                          <  1:242540/215‑1 (MQ=255)
                                                                                                          gcTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGCCCAGATACCACGCCAGCAACATCAGCAGCCATaccactgtctcttctacaaa                                                                                                           >  2:349470/1‑212 (MQ=255)
                                                                                                                  atctgtccGCTGCACTCTTCCCAATCAACCGCCCGCCGGCGAAATTGGTCATCGCAAATCACCCCCAGGTGACCCTGTGGCGTGGCAATTTTGTCGTGAAGTAACCCCAGCATCCGTTCACTCTTAAAAATAACGGGGGTAAATACCCGGACCTAACGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCaa                                                                                                  <  1:621086/221‑1 (MQ=255)
                                                                                                                   tGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAg                                                                                                 >  2:242540/1‑228 (MQ=255)
                                                                                                                   tGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGCGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGACAGCAACATCAGCCACCATAACACCAAAAAGATAAGTGGACAgggaaca                                                                                                 >  2:302854/1‑224 (MQ=255)
                                                                                                                   tGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGCTCTCTCTTAAGACAAACGTGGGTAAATACCCGAGCCTATCGCACAACCGAAAAAATACCACGCCCACAAAATCCACACACACACCACCCAAAAAATaaaaggaaaaggccaa                                                                                                 >  1:447457/1‑214 (MQ=255)
                                                                                                                       ttgccTGTACTCCTCCCTATCACCCGCCCCCAGGTGAAATTTCTAAGCTAAAATCACCCACAGCGTACCCAGTGGCGTGTCAATTTTTTCCTCAAGTAACCGCAGAATCCGTTCGCTCATAAGACAAACGGGGGTAAATACCCGGCCCTATCGGCAGACCGACAAGATGCCACTCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGa                                                                                                <  1:489102/221‑1 (MQ=255)
                                                                                                                        tCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCCCAGCCTCCGTTCTCTCTTAAGACTAACGTGGGTAAAAACACGGACCTATCGTCACACCGTCAAGTTACAACGCCAGCAACATCAGCAGCCATCCCCCCCAAAAGATAAGAGGATAGGGCAAGACtca                                                                                            >  2:129853/1‑226 (MQ=255)
                                                                                                                           ctttccctctCCCATTCAACCCCCCCCAGGCGAAATTTCTTCACACAAATCACCCCCAGTTGACCCAGTTGCGTGGCAATTTTTGCTTCAAGTACTCCCAGCATCCGTTCGCCCGTAAGACAACCGTGGGTAAATACCCGGACCCATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATTAGTGGATAGGGCAAGACTAACGa                                                                                         <  1:628493/219‑1 (MQ=255)
                                                                                                                                ccctcccccTTTAACCCCCCGGAGGCGAAATTGCTCTTCGCAAATCCCCCACAGTGGCCCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCCCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATagt                                                                                    <  1:60525/223‑1 (MQ=255)
                                                                                                                                   ttCCCATTCAACCGCCCGCAGGCGAACTTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCCATTTTTTCTTCAAGTAATCTCAGCATCCGTtctccc                                                                                                                                                                                                           >  1:658006/1‑104 (MQ=255)
                                                                                                                                    tCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTctctct                                                                                                                                                                                                          <  2:658006/106‑1 (MQ=255)
                                                                                                                                      ccATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCAGTGCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCCGCCGACCGACAAGATACCACGCCAGCAACATCAGCAGCCACTCCAACAAAAAGATAAGTGGAAAGGGCAAGACTAACGAAAAGTAAATAc                                                                              >  2:688586/1‑228 (MQ=255)
                                                                                                                                          tCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATTTCAGCATCCGTTCTCTCTTAAGAc                                                                                                                                                                                                    >  2:340193/1‑106 (MQ=255)
                                                                                                                                           cAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCAGCCGTTCTCTCTTAAGACa                                                                                                                                                                                                   <  1:340193/106‑1 (MQ=255)
                                                                                                                                             aCCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTgg                                                                                                           >  2:87112/1‑192 (MQ=255)
                                                                                                                                              ccgcccgcAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCGTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGa                                                                                                          <  1:87112/192‑1 (MQ=255)
                                                                                                                                                gcccgcAGGCTAAATTGCTCATCGAAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCCCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTcgc                                                                    <  2:186276/228‑1 (MQ=255)
                                                                                                                                                  ccgcAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGt                                                                  <  2:9433/228‑1 (MQ=255)
                                                                                                                                                  ccgcAGGCGAAATTGCTCATCGCAAATCACCCACAGTGCACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTGAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGt                                                                  <  2:419121/228‑1 (MQ=255)
                                                                                                                                                   cgaaGGGGAAATTGCATATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGGGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTc                                                                 <  2:665140/225‑1 (MQ=255)
                                                                                                                                                                                           cccgtggcgtggcAATTTTTTCTAAAAGTAATCTCAGCCTCCGTTCTCGCTAAAAAAAAACGTGGGGAAAGACCCGGCCCTACCGGCAAACCAACAAGATGCCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACa                         <  1:136301/225‑1 (MQ=255)
                                                                                                                                                                                             aGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCCCGCCAGCAACATCAGCAGCCCTACCACCAAAAACATAAGTGGATAGGGCAAGACTAACGCAAAGTAAATACCCCCGTTCCCAGCTGGTTTGTAACGTTGCAGGAAACCAAAAAACAGAGGaaaaaa                       >  1:143680/1‑228 (MQ=255)
                                                                                                                                                                                              gtggcgtggcAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAAGGCAAGACTAACGAACAGTAAGTACCCAGTTTCGCGGCTGGTTTGTAGGGTTTCGAGAAAACAAAAAAAAGAGGGACaaaa                      >  1:527879/1‑228 (MQ=255)
                                                                                                                                                                                              gtggcgtggcAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAAAATCAACAGCCATACCACCAAAACAAAAAGTGGAAAGGGCAAGAATAACGAATAAAAAATGCCCAAGTTCCCGTCCGGTGTGTAGCGGTGCAGGAATCCAAAAAAAAAGGGAaaaaaa                      >  1:358876/1‑228 (MQ=255)
                                                                                                                                                                                                 gcgtggcAATTTTTTCTTCAAGTAATCCCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGt                                                                  <  2:305975/181‑1 (MQ=255)
                                                                                                                                                                                                      gcAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATAc              <  2:629089/228‑1 (MQ=255)
                                                                                                                                                                                                               gtcttcAAGTAACCTCCGCCTCCCTCCCCCCGAAAAACAAACGTGGATAAAGCCCCGGACCTACCGGCAAACCGAAAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCt       <  1:548249/226‑1 (MQ=255)
                                                                                                                                                                                                                ttttaaaGTAAGCCCACCATCCGTTCTCTCTTAAGACAACCGGGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTaa     <  1:395859/223‑1 (MQ=255)
                                                                                                                                                                                                                 cttcAAGTAATCTCAGCATCCGTTCTCCCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTaaa    <  2:527879/228‑1 (MQ=255)
                                                                                                                                                                                                                  ttcAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGTACCTATCGGCAGACCGACAAGATCCCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTaaa    <  2:462813/227‑1 (MQ=255)
                                                                                                                                                                                                                   tcAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGATTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACATTGGAACAAAAGGAGAAACCGGCGACAAAgg  >  2:295604/1‑228 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                          
ATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGG  >  NC_000913/1400010‑1400447

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