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

ATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGG  >  NC_000913/1400010‑1400447
                                                                                                                                                                                                                              |                                                                                                                                                                                                                          
atGCAAAGATTACCGGGAAAATATTACCGAAGGTTGTCCTTTAAACCGCCTGGATTGGTTGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGGACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCAt                                                                                                                                                                                                                      <  1:222069/227‑1 (MQ=255)
   cTAAGATTACCGGCAAAATATTCACGAAACTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCCCTGTACTCTTCCCATTCAACCGCCCGCCGGCGAAATTGCTCATCGCCAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGAAA‑CCCCCGCCTCCGt                                                                                                                                                                                                                  >  2:408646/1‑228 (MQ=255)
      aGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTTAAGTAA‑ATCCAGCATCCGt‑ttt                                                                                                                                                                                                              >  1:73953/1‑226 (MQ=255)
        aTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCCGGCGAAATTGCTCATCGCCAATCACCCACAGTGGCCCCAGTGGCGTGGCAAT‑TTTTTCTTCAATAAA‑ACCCACCATCCGc‑tctct                                                                                                                                                                                                            >  2:137324/1‑228 (MQ=255)
        aTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGCAATTGCTCATCGCAACTCACTCACAGTGGACCCAGTGGCGTGGCAAC‑TTTGTCTTTAAGAAA‑TCTCACCATCCCt‑ttcc                                                                                                                                                                                                             >  1:404377/1‑224 (MQ=255)
        aTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGGGGACCCAGTGGCCGGGCCAT‑TTTTTTTTAAAGTAACTCCCAACACCCGT‑tctc                                                                                                                                                                                                             >  2:629089/1‑228 (MQ=255)
          ttcgggCAAAATTTTCCCGAATTTGTTCGTTTAACCCGCCTGGATTGGGGGCAGAAATGCGCTCATAGCCTTCTTTGCCATAATGTATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑Tctctct                                                                                                                                                                                                          <  2:393636/224‑1 (MQ=255)
           aCCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCCCAGTGGACCCAGTGGCGGGGCAAT‑TTTTTCTTCCAGTAA‑TCTCACCATcccgttcctctt                                                                                                                                                                                                          >  1:632979/1‑218 (MQ=255)
            ccGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATGTTTTTCTTC‑AGTAA‑TTCCAGCACCCCT‑TGCCTCttc                                                                                                                                                                                                        >  1:251931/1‑227 (MQ=255)
                 aaaaTATTCACGAAGTTTTTCGTTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACa                                                                                                                                                                                                   <  1:92537/228‑1 (MQ=255)
                     tatTCCCGAAGGTTAGCGCTTAAACCGCCTTGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTTCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACg                                                                                                                                                                                               <  1:435284/228‑1 (MQ=255)
                          aCGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑ATTCAGCAACCGT‑CTTTCCTTAAGCAAAACGTgcg                                                                                                                                                                                           >  1:157026/1‑225 (MQ=255)
                                   tccGCTAAACCCGCCTGGATTGGTGGCAGAAATGCCCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCg                                                                                                                                                                                 <  1:137324/226‑1 (MQ=255)
                                        ttAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGACAT‑TTTTTTTTTCAGGAA‑TCTCCGCCTCCGG‑TCTCTCTTTAGACACACGTGGGTAAATAccccgtcc                                                                                                                                                                             >  1:136301/1‑222 (MQ=255)
                                                ccTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAATTCACCCACAGTGGACCCAGTGGCGTGGAAAT‑TTTTTTTTCAAGTAA‑TCTCAGCATCCCT‑CCACCCTTAAAACAACCGTGGGCAAATACCCCGAccatccgtgag                                                                                                                                                                    >  2:73948/1‑219 (MQ=255)
                                                  tGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCCCTTAAGACCAACGTGGGTAAATACCCGCACATATCgggaca                                                                                                                                                                   >  2:419121/1‑223 (MQ=255)
                                                     aTTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCCGTGGCGTGGCAAT‑TTTTTCCTTAAGTAA‑ACTCACCACCCCT‑TCTCCCTTCAGCAAACCCTGGGTAAATAACCGGACCCACCGCCAAAccc                                                                                                                                                                >  2:9433/1‑226 (MQ=255)
                                                     aTTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGa                                                                                                                                                               <  1:398716/228‑1 (MQ=255)
                                                        ggtggCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGt‑t                                                                                                                                                                                                                >  1:383899/1‑176 (MQ=255)
                                                         gtggCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑tc                                                                                                                                                                                                               <  2:383899/176‑1 (MQ=255)
                                                          tggCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCCCAGTGGCCCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGGAA‑TCCCCGCATCCGT‑TCTCTCTTTAGACAAACGTGGGGTAATCCCCGAGCCAAACGGCAGAccccacaac                                                                                                                                                          >  2:609830/1‑221 (MQ=255)
                                                             ccgAAATGCCCTCATAGCCCTCTTTGCGATACAGGCTGCCCAGCAGCTGCCCCATCCGTTCCCTGTACTCTCCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCTTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATAc                                                                                                                                                       <  2:157026/226‑1 (MQ=255)
                                                                  aTGCGCTCATAGCCTTTGTTGCGATAATGGATGTCCAGCAGCTGCCCCATGCGTTCGCTGGACTCTTCCCATACAACCGCCCGGAGGCGAAATTGCTCATCGCAAATCGTCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACg                                                                                                                                                   <  1:394208/227‑1 (MQ=255)
                                                                          aTAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACCTTGGTAAATACCCGGACCTATCGGCAGACCGCCACGATACCCACCCCGCAACa                                                                                                                                          >  1:60525/1‑228 (MQ=255)
                                                                          aTAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑ACTTAGCATCCGT‑TCTCTCTTAAGACAAACATGGGTGAATACCCGGACCTATCTGCAGACCGACAAGATACCACGCCAGCAACa                                                                                                                                          >  1:489102/1‑228 (MQ=255)
                                                                            tcccTTTTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATc                                                                                                                                        <  2:632979/226‑1 (MQ=255)
                                                                                  cttTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTTTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGGGGGTAAAAACCCGGACCCATCGGCCGACCGACAAGATACCACGCCAGCAACATCagcagc                                                                                                                                  >  1:552114/1‑228 (MQ=255)
                                                                                  cttTGCGATAATGGATGCCCAGCAGCTGAACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGGGGACAT‑TTTTTCTTCAAGGAC‑TCTTATGATCCGTATTGCTCTAATGA‑AAACGTGGTGGAAAACCCGGACCATACGGCAGAACCTCAAGACAACCCGCCGCCAAAAAAagcagc                                                                                                                                  >  2:500812/1‑228 (MQ=255)
                                                                                       cGATAATGGATGTCCAGCAGCTGCACCATGAGTTCGCGGTACTCTTCCAATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATac                                                                                                                             <  1:408646/228‑1 (MQ=255)
                                                                                           aagggATGCCCACCCGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCGACCCCCCGCAGGCGAAATTGCTCATCGCAAATCCCCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATCCCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATaccacc                                                                                                                         <  2:552114/225‑1 (MQ=255)
                                                                                              ggATGCCCAGCCACTGCACCATGCGTTCGCTTTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCCCCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCaaa                                                                                                                      <  1:141914/228‑1 (MQ=255)
                                                                                              ggATGCCCAGCACCTTCCCCATGCGTTCGATGTCCTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCCCCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCaaa                                                                                                                      <  1:609830/228‑1 (MQ=255)
                                                                                               gATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCCGCATCCGT‑TTTCTCTTAAAACAAACGTGGGGAAATACCCGGACCTATCGGCATGCCGACCAAATACAACGCCAGCAACATCAGCAGCACTACCACCaaaa                                                                                                                     >  1:354633/1‑228 (MQ=255)
                                                                                                                  atctgtccGCTGCACTCTTCCCAATCAACCGCCCGCCGGCGAAATTGGTCATCGCAAATCACCCCCAGGTGACCCTGTGGCGTGGCAAT‑TTTGTCGTGAAGTAA‑CCCCAGCATCCGT‑TCACTCTTAAAAATAACGGGGGTAAATACCCGGACCTAACGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCaa                                                                                                  <  2:621086/221‑1 (MQ=255)
                                                                                                                   tGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAg                                                                                                 >  1:242540/1‑228 (MQ=255)
                                                                                                                   tGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGCGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGACAGCAACATCAGCCACCATAACACCAAAAAGATAAGTGGACAgggaaca                                                                                                 >  1:302854/1‑224 (MQ=255)
                                                                                                                       ttgccTGTACTCCTCCCTATCACCCGCCCCCAGGTGAAATTTCTAAGCTAAAATCACCCACAGCGTACCCAGTGGCGTGTCAAT‑TTTTTCCTCAAGTAA‑CCGCAGAATCCGT‑TCGCTCATAAGACAAACGGGGGTAAATACCCGGCCCTATCGGCAGACCGACAAGATGCCACTCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGa                                                                                                <  2:489102/221‑1 (MQ=255)
                                                                                                                        tCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCCCAGCCTCCGT‑TCTCTCTTAAGACTAACGTGGGTAAAAACACGGACCTATCGTCACACCGTCAAGTTACAACGCCAGCAACATCAGCAGCCATCCCCCCCAAAAGATAAGAGGATAGGGCAAGACtca                                                                                            >  1:129853/1‑226 (MQ=255)
                                                                                                                           ctttccctctCCCATTCAACCCCCCCCAGGCGAAATTTCTTCACACAAATCACCCCCAGTTGACCCAGTTGCGTGGCAAT‑TTTTGCTTCAAGTAC‑TCCCAGCATCCGT‑TCGCCCGTAAGACAACCGTGGGTAAATACCCGGACCCATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATTAGTGGATAGGGCAAGACTAACGa                                                                                         <  2:628493/219‑1 (MQ=255)
                                                                                                                                ccctcccccTTTAACCCCCCGGAGGCGAAATTGCTCTTCGCAAATCCCCCACAGTGGCCCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCCCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATagt                                                                                    <  2:60525/223‑1 (MQ=255)
                                                                                                                                   ttCCCATTCAACCGCCCGCAGGCGAACTTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCCAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑tctccc                                                                                                                                                                                                           >  2:658006/1‑104 (MQ=255)
                                                                                                                                    tCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑Tctctct                                                                                                                                                                                                          <  1:658006/106‑1 (MQ=255)
                                                                                                                                      ccATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCAGT‑GCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCCGCCGACCGACAAGATACCACGCCAGCAACATCAGCAGCCACTCCAACAAAAAGATAAGTGGAAAGGGCAAGACTAACGAAAAGTAAATAc                                                                              >  1:688586/1‑228 (MQ=255)
                                                                                                                                          tCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TTTCAGCATCCGT‑TCTCTCTTAAGAc                                                                                                                                                                                                    >  1:340193/1‑106 (MQ=255)
                                                                                                                                           cAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCAGCCGT‑TCTCTCTTAAGACa                                                                                                                                                                                                   <  2:340193/106‑1 (MQ=255)
                                                                                                                                             aCCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTgg                                                                                                           >  1:87112/1‑192 (MQ=255)
                                                                                                                                              ccgcccgcAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCGTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGa                                                                                                          <  2:87112/192‑1 (MQ=255)
                                                                                                                                                gcccgcAGGCTAAATTGCTCATCGAAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCCCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTcgc                                                                    <  1:186276/228‑1 (MQ=255)
                                                                                                                                                  ccgcAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGt                                                                  <  1:9433/228‑1 (MQ=255)
                                                                                                                                                  ccgcAGGCGAAATTGCTCATCGCAAATCACCCACAGTGCACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTGAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGt                                                                  <  1:419121/228‑1 (MQ=255)
                                                                                                                                                   cgaaGGGGAAATTGCATATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGGGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTc                                                                 <  1:665140/225‑1 (MQ=255)
                                                                                                                                                                                           cccgtggcgtggcAAT‑TTTTTCTAAAAGTAA‑TCTCAGCCTCCGT‑TCTCGCTAAAAAAAAACGTGGGGAAAGACCCGGCCCTACCGGCAAACCAACAAGATGCCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACa                         <  2:136301/225‑1 (MQ=255)
                                                                                                                                                                                             aGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCCCGCCAGCAACATCAGCAGCCCTACCACCAAAAACATAAGTGGATAGGGCAAGACTAACGCAAAGTAAATACCCCCGTTCCCAGCTGGTTTGTAACGTTGCAGGAAACCAAAAAACAGAGGaaaaaa                       >  2:143680/1‑228 (MQ=255)
                                                                                                                                                                                              gtggcgtggcAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAAGGCAAGACTAACGAACAGTAAGTACCCAGTTTCGCGGCTGGTTTGTAGGGTTTCGAGAAAACAAAAAAAAGAGGGACaaaa                      >  2:527879/1‑228 (MQ=255)
                                                                                                                                                                                              gtggcgtggcAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAAAATCAACAGCCATACCACCAAAACAAAAAGTGGAAAGGGCAAGAATAACGAATAAAAAATGCCCAAGTTCCCGTCCGGTGTGTAGCGGTGCAGGAATCCAAAAAAAAAGGGAaaaaaa                      >  2:358876/1‑228 (MQ=255)
                                                                                                                                                                                                 gcgtggcAAT‑TTTTTCTTCAAGTAA‑TCCCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGt                                                                  <  1:305975/181‑1 (MQ=255)
                                                                                                                                                                                                      gcAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATAc              <  1:629089/228‑1 (MQ=255)
                                                                                                                                                                                                               gtcttcAAGTAA‑CCTCCGCCTCCCT‑CCCCCCGAAAAACAAACGTGGATAAAGCCCCGGACCTACCGGCAAACCGAAAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCt       <  2:548249/226‑1 (MQ=255)
                                                                                                                                                                                                                ttttaaaGTAA‑GCCCACCATCCGT‑TCTCTCTTAAGACAACCGGGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTaa     <  2:395859/223‑1 (MQ=255)
                                                                                                                                                                                                                 cttcAAGTAA‑TCTCAGCATCCGT‑TCTCCCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTaaa    <  1:527879/228‑1 (MQ=255)
                                                                                                                                                                                                                  ttcAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGTACCTATCGGCAGACCGACAAGATCCCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTaaa    <  1:462813/227‑1 (MQ=255)
                                                                                                                                                                                                                   tcAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGATTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACATTGGAACAAAAGGAGAAACCGGCGACAAAgg  >  1:295604/1‑228 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                          
ATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAAT‑TTTTTCTTCAAGTAA‑TCTCAGCATCCGT‑TCTCTCTTAAGACAAACGTGGGTAAATACCCGGACCTATCGGCAGACCGACAAGATACCACGCCAGCAACATCAGCAGCCATACCACCAAAAAGATAAGTGGATAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGG  >  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: ‑