Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,400,230 T→C 17.5% 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,2300TC17.5% 110.8 / 3.2 62R3G (AGA→GGA) ogtO‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (23/26);  new base C (7/5);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

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