Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,400,026 C→A 13.4% A71S (GCC→TCC)  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,0260CA13.4% 187.7 / 2.0 70A71S (GCC→TCC) ogtO‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (36/23);  new base A (2/9);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATGGTGCCGTTTCGGCCAATAACCCGATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGT  >  NC_000913/1399789‑1400259
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
tattGGTGCCGTTTCGGCCAATAACCCGATGGCAAGGGACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGc                                                                                                                                                                                                                                           <  1:502968/235‑1 (MQ=255)
        ccGTTTCGGCCAATAACCCGATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCCTTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGCCTTCGCGCTGAAATGGCGTCCCCCCCGTTGCAGTGGGAAGCGTACACAAAATGCCAAGATTACCGCCAAAAAATTc                                                                                                                                                                                                                                  >  2:79861/1‑239 (MQ=255)
                 cccaTACCCCGGTGGCAAGGTTCGACGATGCGTATGGGATTGAGTCCGTTGCCCGCCGCAACGCCCCGCGCCGGCCCATGACGGACCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                                                                                                                                                                                         <  2:239623/236‑1 (MQ=255)
                 cccaTACCCCGGTGGCAAGGTACGCCGATGCTGATGGGATTCGATCCGGTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                                                                                                                                                                                         <  2:517548/236‑1 (MQ=255)
                   aaTAACCCGATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGt                                                                                                                                                                                                                       <  1:716234/239‑1 (MQ=255)
                     tAACCCGATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGCCTTCGCGCTGACATGGCGTCCCCCCCGTCGCAGTGGGAAGCGTATCAATAATGCTAAGGTTCCCGGCCAAATATTCACCAACATTGTCg                                                                                                                                                                                                                     >  1:144785/1‑239 (MQ=255)
                           gATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGGGCGTAGTGTTTTCCAGACTTCCCGCTGAAATGGCGTCCCCCCCGTAGCAGGGGGAAGGGTATCAATAATGCAAAGATTACCGGAAAAATATTTACGAAGCTTGTCCCTTaaa                                                                                                                                                                                                               >  2:11452/1‑239 (MQ=255)
                                   ggTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTgg                                                                                                                                                                                                       >  1:897295/1‑239 (MQ=255)
                                   ggTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGGATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTgg                                                                                                                                                                                                       >  1:478367/1‑239 (MQ=255)
                                             tGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCCGGGGCTAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAACATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCa                                                                                                                                                                                             >  1:377628/1‑239 (MQ=255)
                                               cTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCCAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGa                                                                                                                                                                                           >  1:379631/1‑239 (MQ=255)
                                                       aTTCGGTCCGTTTGCCCGACGAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                                                                                                                                                                                    <  1:574491/238‑1 (MQ=255)
                                                       aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGGGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCCGACAAATATTCCCGAAGCTTGTCGCTTAAACCGCCCGGATGGGTGGCAGAAATgggc                                                                                                                                                                                    >  2:796722/1‑238 (MQ=255)
                                                             ctcttttgccacccacgcGGCGCGCGCCGGGCCGGGCCGGCCAATTTGCTCAGCCAGTTGGCTGAAATGCTTTCCCTGCCCGCAGGGGATAGTGCGTAGTGTTTTACAGACGTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATCCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGc                                                                                                                                                                             <  2:377628/222‑1 (MQ=255)
                                                                 tttGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCACTAATGCTAAGCTTACCGGCAAAATATTCACGGAGCTTTTCGCTTAAACCCCCTGGACTGGTGGCCGAAATGCGCTCATAGCctct                                                                                                                                                                         >  2:239875/1‑237 (MQ=255)
                                                                      cGCACCAACGGCACGCGCCGCGCCGGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAGGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTg                                                                                                                                                                     <  1:124184/238‑1 (MQ=255)
                                                                      cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTg                                                                                                                                                                     <  1:11452/238‑1 (MQ=255)
                                                                               ggccggcgtcgcgccATGACGGCAAATGTGATCGGCCAGTGGCCGGAATGGCATTCCCTGCCGGCAGGGGATAGTGAGGAGGGGTGTGAAGACTTCGGGAGGAAAGGGCGTCCCCCCAGTAGGAGGGGGAAGCGTATCAATAATGCAAAGATTACCGGAAAAATATTCCCGAAGCTGGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATgg                                                                                                                                                            <  2:144785/233‑1 (MQ=255)
                                                                                    gcgccgcgccAGGAGGGCCCAATTGCACAGCCAGTTGGCCGTAATGCATTCCCGGCCCGCAGGGTATAGGGCGTAGGGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGGGGGAAGCGTATCAATAATGCTAAGATTACCGGAAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc                                                                                                                                                      <  1:239875/239‑1 (MQ=255)
                                                                                    gcgccgagccAGGCCGGCCCAATTGCTCAGCCAGTTGGCCGTAAGTCCTTACCTGCCCGGAGGGGATAGTGCGTAGTTTTTTCCAGACTTCGCGCTGAAAGGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGGCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc                                                                                                                                                      <  1:645167/239‑1 (MQ=255)
                                                                                    cgcccgcgccAGGCCGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGAAGGGGATAGTGCGTAGTGTTTCCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGGGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc                                                                                                                                                      <  1:576446/236‑1 (MQ=255)
                                                                                    cgcccccgccgggcctgccataTTTCTCAGCCTGTTGGCCGTAATGCATTACCTGCCCCGGGGGGATAGGGCTTAGTGTTTTCCAGCCTTCGCGCTGAAATGGCGTCCCCCCCGTAGCATTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc                                                                                                                                                      <  2:133715/218‑1 (MQ=255)
                                                                                     cgccgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAAAGCGCCCCTAGCCTTCTTTGCGAAAAAGGGTTTccc                                                                                                                                                     >  2:836007/1‑238 (MQ=255)
                                                                                         gcgccGGGCCGGCCCAATTGCCCGGCCGTTTGGCCGTAAAGCATTCCCGGCCCGGAGGGGATAGTGCGGAGTGTTTTCCAGACTTCGCGCTGAATGGGCGTCCCCCCCGTAGCAGTGGGAAGCGTTTCAATAATGCTAAGATGACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCcagcag                                                                                                                                                 <  2:635126/239‑1 (MQ=255)
                                                                                           gccggggcgggcgatttGCTGAGCCGTGTGGCCGACATTCATTACCGGCCCGCAGGGGATGGGGCGTGGTGTTGTCCAGACGTCGCGCGGAAATGGCCTCCCCCCAGTGGAAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGAAAAATATTCAAGAAGCTTGCCGCTTAAACCGCCTGGATTGGTGGCAGAAAGGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCagcagc                                                                                                                                                <  2:635637/223‑1 (MQ=255)
                                                                                                    gCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGTGCTCCTAGCCTTCTTTTCGATAATGGATGTCCAGCACCTGCCCCCTGc                                                                                                                                      >  2:642113/1‑239 (MQ=255)
                                                                                                     cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATTGTGCGTAGGGTTTTCCAGACTTCGCGCTGTAATGGCGTCACCCCCGTAGCCGTGGGAAGCGTATCAATAATGCTCAGATTACCTGCAAAATATTCACGCAGCTTGTCCCTTAAACCGCCTGGATTGGTGGCAGAAAATCGCGCAATCCCTTCTTTGCGAAAATGGATggcaaacacgtgccccatacc                                                                                                                                     >  2:856764/1‑219 (MQ=255)
                                                                                                     cccAATTGCCCAGCCAGGTGGCCGTAATGCATTCCCTGCCCGCAGGGGATAGTGCGTAGTTTTTTCCAGTCTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGc                                                                                                                                      <  1:561066/238‑1 (MQ=255)
                                                                                                         aTTGGTCAGCCAGTTGGCCGTAATGCATTGCCGGCCCGCAGGGGTTAGTGCGTGGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCg                                                                                                                                 <  1:79861/239‑1 (MQ=255)
                                                                                                               cAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTCAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGt                                                                                                                             >  1:279183/1‑237 (MQ=255)
                                                                                                               cAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAAAATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATCCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCCGCACCATGCGTTCCCTgga                                                                                                                            >  2:604463/1‑236 (MQ=255)
                                                                                                                  ccAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTActc                                                                                                                         >  1:91656/1‑238 (MQ=255)
                                                                                                                       cggctctAATGCATTCCCTGCCCGCAGGGCATAGGGCGTAGTGTTTTCCAGACTTCGCGTTGAAATGGCGTCCCCCCCGTAGCGGTGGGAAGCGTATCAATAATGCTAAGTTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCa                                                                                                                   <  1:253734/233‑1 (MQ=255)
                                                                                                                           cctAATTCCTTTCCCGCCCGGAGGGGGTAGTGGGTGTTTTGTTCCAGACTTGCGGCTTGAATGGGCGCCCCCCCCGTGGAGGGGGGAGCGTTACCAAAAAGTCTAAGATTCCGGGAAAAATATTCAGGAAGTTTGCCGTTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCAtt                                                                                                                 <  2:897295/236‑1 (MQ=255)
                                                                                                                                    tttcccTCCCCGGGGGGGATTGGGGGTTGTGTTTCCCAGTCTTCGCCGTGAAATGGCGTCCCCCCCGGAGCAGGGGGAAGCTTATCAATAATGCTAAGATCCCCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAAccgcccg                                                                                                       <  2:408567/235‑1 (MQ=255)
                                                                                                                                      tACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCa                                                                                                     >  1:62955/1‑239 (MQ=255)
                                                                                                                                           gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCCTGCGTTCGCTGTACTCTTCCCATTCAACCGACCCAAGGCGa                                                                                                >  2:747195/1‑239 (MQ=255)
                                                                                                                                            cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGaa                                                                                               >  1:277708/1‑239 (MQ=255)
                                                                                                                                                cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAAtt                                                                                            >  1:590322/1‑238 (MQ=255)
                                                                                                                                                        gtgggtgtgtgtTTTCCAGACTTCGCGGTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTAGCAATAATGCTAAGAGTCCCGGAAAAATATTCGCGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGc                                                                                   <  2:377873/231‑1 (MQ=255)
                                                                                                                                                             ttAGTTTTTTCCGGTCTTCGCGTTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTAACAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAAt                                                                               <  2:379631/237‑1 (MQ=255)
                                                                                                                                                               aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCAACGCAAATCa                                                                             >  2:12473/1‑238 (MQ=255)
                                                                                                                                                                 tgtttgccagtattcggcgtggAAATGGCTTCCCCCCCGTGTCGGGGGGAAGGGTATCAATAATGCTAATATTACCGGCAATATTTTAGCGAAGCGTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCAcc                                                                           <  2:279183/218‑1 (MQ=255)
                                                                                                                                                                  tttttCCAGCCTTCGGCCTGAAATGGCTCCCCCCCCGTAGCAGTGGGAAGCGTAACAATAATGCTAAGATTACCGGAAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCAccc                                                                          <  1:836007/237‑1 (MQ=255)
                                                                                                                                                                    tttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCcagcag                                                                                                                                                 >  2:686751/1‑165 (MQ=255)
                                                                                                                                                                    tttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCcagcag                                                                                                                                                 <  1:686751/165‑1 (MQ=255)
                                                                                                                                                                      tCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATTCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGt                                                                     >  2:650341/1‑239 (MQ=255)
                                                                                                                                                                      tCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGt                                                                     >  1:499322/1‑239 (MQ=255)
                                                                                                                                                                           aCTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCCCTGTACTCTTCCCATTAAACCGCCCGCAGGCGAAATTGCCCACGCaaaaacaccaccagggcac                                                                >  2:833304/1‑222 (MQ=255)
                                                                                                                                                                                cgcgCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCAACAGTGGACCCAgt                                                            >  1:777822/1‑239 (MQ=255)
                                                                                                                                                                                                 ccccccGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGCAATTGCTCATCGCAAATACCCCACAGTGGACCCAGTGGCGGGGCAAatttttt                                           >  1:559283/1‑238 (MQ=255)
                                                                                                                                                                                                         aGCAGAGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGACTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGACCTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCCCAGTGGACCCAGTGGCGTGGCAAATTTttcttacaaaa                                   >  2:152441/1‑233 (MQ=255)
                                                                                                                                                                                                           cAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAACCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTAAACCGCCCCCAGGCGAAATTGCTCATCGCCAATCACCCACACTGGAACCAGTAGCGGGGTGATTTTTTTTTTAAGTaa                                  >  1:761539/1‑238 (MQ=255)
                                                                                                                                                                                                              tGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGGTTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCACATCACCCACAGTGGACCAAGGGGCGTGGCAATTTTTTTTTCaaaaaaaccc                              >  2:374145/1‑231 (MQ=255)
                                                                                                                                                                                                              tGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCCCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCCCATCGCAAATCACCCACAGTGGACCCCGTGGGGTGGCAAATTTTTTTTCaaaaacactc                              >  2:5669/1‑231 (MQ=255)
                                                                                                                                                                                                               gggAAGCTCTTCAATAATGTTAAGATTCCGGGAAAAATATTCCCGAAGCTTGGCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAAtctc                              <  2:78616/238‑1 (MQ=255)
                                                                                                                                                                                                               gggAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAAtctc                              >  1:771672/1‑238 (MQ=255)
                                                                                                                                                                                                                ggAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCa                             >  1:78616/1‑238 (MQ=255)
                                                                                                                                                                                                                ggAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGGGGCAAATTTTTCTTTAAGTAATCTCa                             >  1:563765/1‑238 (MQ=255)
                                                                                                                                                                                                                 gAAGCGTAAAAAAAATGCTAAGGTTACCGGCAAAATATTCCCGAAGCTTGCCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAg                            <  1:374145/238‑1 (MQ=255)
                                                                                                                                                                                                                         aTCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGCCGGGGCAATTTTTTTTTCAAGAAACCCCAACAACCGt                     >  2:466579/1‑238 (MQ=255)
                                                                                                                                                                                                                          tcaatttttcTAAGTTTCCCGGAAAAATTTTCACGAACTTTGTCTCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTtc                   <  2:91656/230‑1 (MQ=255)
                                                                                                                                                                                                                          tCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTtc                   <  1:604463/239‑1 (MQ=255)
                                                                                                                                                                                                                          tCAATAATGCAAAGATTACCGGAAAAATATTCACGAACCTTGGCCCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTtc                   <  1:642113/239‑1 (MQ=255)
                                                                                                                                                                                                                          tCAATAAGGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTtc                   <  1:796722/239‑1 (MQ=255)
                                                                                                                                                                                                                                 tGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGGTCATCGCAAATCAACCCCACGGGACCCAGTGGCGGCGACATTTTTTCTTCAAATAATCccaccacacagtctccctta            >  2:840993/1‑225 (MQ=255)
                                                                                                                                                                                                                                     aaGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGa         <  1:833304/238‑1 (MQ=255)
                                                                                                                                                                                                                                      aGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAACATCCGTTCTCTCTTAAGACa       >  1:588830/1‑239 (MQ=255)
                                                                                                                                                                                                                                        aTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCCCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGGAATCCCAAAATCCGTTCTCCCTTAAGAaaaa     >  2:819025/1‑239 (MQ=255)
                                                                                                                                                                                                                                          ttccGGCAAAATTTTCACGAGTCTTTTCTTTTAAACCGCCTGGATTGGGGGCAGAAATGCGCTCATAGCCGTCTTTTCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACg   <  2:478367/237‑1 (MQ=255)
                                                                                                                                                                                                                                           aCCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCCATTTTTTTTTCCAGTAATCTCAACATCCGTTCTCCCTTAAGACAAACGt  >  1:818161/1‑239 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
TCATGGTGCCGTTTCGGCCAATAACCCGATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGT  >  NC_000913/1399789‑1400259

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