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

TCATGGTGCCGTTTCGGCCAATAACCCGATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGT  >  NC_000913/1399789‑1400259
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
tattGGTGCCGTTTCGGCCAATAACCCGATGGCAAGGGACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGc                                                                                                                                                                                                                                           <  1:502968/235‑1 (MQ=255)
        ccGTTTCGGCCAATAACCCGATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCCTTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGCCTT‑CGCGCTGAAATGGCGTCCCCCCCGTTGCAGTGGGAAGCG‑TACACAAAATGCCAAGATTACCGCCAAAAAATTc                                                                                                                                                                                                                                  >  2:79861/1‑239 (MQ=255)
                 cccaTACCCCGGTGGCAAGGTTCGACGATGCGTATGGGATTGAGTCCGTTGCCCGCCGCAACGCCCCGCGCCGGCCCATGACGGACCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                                                                                                                                                                                         <  2:239623/236‑1 (MQ=255)
                 cccaTACCCCGGTGGCAAGGTACGCCGATGCTGATGGGATTCGATCCGGTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                                                                                                                                                                                         <  2:517548/236‑1 (MQ=255)
                   aaTAACCCGATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGt                                                                                                                                                                                                                       <  1:716234/239‑1 (MQ=255)
                     tAACCCGATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGCCTT‑CGCGCTGACATGGCGTCCCCCCCGTCGCAGTGGGAAGCG‑TATCAATAATGCTAAGGTTCCCGGCCAAATATTCACCAACATTGTCg                                                                                                                                                                                                                     >  1:144785/1‑239 (MQ=255)
                           gATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGGGCGTAGTGTTTTCCAGACTT‑CCCGCTGAAATGGCGTCCCCCCCGTAGCAGGGGGAAGGG‑TATCAATAATGCAAAGATTACCGGAAAAATATTTACGAAGCTTGTCCCTTaaa                                                                                                                                                                                                               >  2:11452/1‑239 (MQ=255)
                                   ggTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTgg                                                                                                                                                                                                       >  1:897295/1‑239 (MQ=255)
                                   ggTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑GATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTgg                                                                                                                                                                                                       >  1:478367/1‑239 (MQ=255)
                                             tGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCCGGGGCTAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAACATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCa                                                                                                                                                                                             >  1:377628/1‑239 (MQ=255)
                                               cTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCCAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGa                                                                                                                                                                                           >  1:379631/1‑239 (MQ=255)
                                                       aTTCGGTCCGTTTGCCCGACGAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                                                                                                                                                                                    <  1:574491/238‑1 (MQ=255)
                                                       aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGGGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCCGACAAATATTCCCGAAGCTTGTCGCTTAAACCGCCCGGATGGGTGGCAGAAATgggc                                                                                                                                                                                    >  2:796722/1‑238 (MQ=255)
                                                             ctcttttgccacccacgcGGCGCGCGCCGGGCCGGGCCGGCCAATTTGCTCAGCCAGTTGGCTGAAAT‑GCTTTCCCTGCCCGCAGGGGATAGTGCGTAGTGTTTTACAGACGT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATCCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGc                                                                                                                                                                             <  2:377628/222‑1 (MQ=255)
                                                                 tttGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCACTAATGCTAAGCTTACCGGCAAAATATTCACGGAGCTTTTCGCTTAAACCCCCTGGACTGGTGGCCGAAATGCGCTCATAGCctct                                                                                                                                                                         >  2:239875/1‑237 (MQ=255)
                                                                      cGCACCAACGGCACGCGCCGCGCCGGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAG‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTg                                                                                                                                                                     <  1:124184/238‑1 (MQ=255)
                                                                      cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTg                                                                                                                                                                     <  1:11452/238‑1 (MQ=255)
                                                                               ggccggcgtcgcgccATGACGGCAAATGTGATCGGCCAG‑TGGCCGGAATGGCATTCCCTGCCGGCAGGGGATAGTGAGGAGGGGTGTGAAGACTT‑CGGGAGGAAAGGGCGTCCCCCCAGTAGGAGGGGGAAGCG‑TATCAATAATGCAAAGATTACCGGAAAAATATTCCCGAAGCTGGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATgg                                                                                                                                                            <  2:144785/233‑1 (MQ=255)
                                                                                    gcgccgcgccAGGAGGGCCCAATTGCACAGCCAGTTGGCCGTAAT‑GCATTCCCGGCCCGCAGGGTATAGGGCGTAGGGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGGGGGAAGCG‑TATCAATAATGCTAAGATTACCGGAAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc                                                                                                                                                      <  1:239875/239‑1 (MQ=255)
                                                                                    gcgccgagccAGGCCGGCCCAATTGCTCAGCCAGTTGGCCGTAAG‑TCCTTACCTGCCCGGAGGGGATAGTGCGTAGTTTTTTCCAGACTT‑CGCGCTGAAAGGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGGCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc                                                                                                                                                      <  1:645167/239‑1 (MQ=255)
                                                                                    cgcccgcgccAGGCCGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGAAGGGGATAGTGCGTAGTGTTTCCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGGGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc                                                                                                                                                      <  1:576446/236‑1 (MQ=255)
                                                                                    cgcccccgccgggcctgccataTTTCTCAGCCTGTTGGCCGTAAT‑GCATTACCTGCCCCGGGGGGATAGGGCTTAGTGTTTTCCAGCCTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCATTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc                                                                                                                                                      <  2:133715/218‑1 (MQ=255)
                                                                                     cgccgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAAAGCGCCCCTAGCCTTCTTTGCGAAAAAGGGTTTccc                                                                                                                                                     >  2:836007/1‑238 (MQ=255)
                                                                                         gcgccGGGCCGGCCCAATTGCCCGGCCGTTTGGCCGTAAA‑GCATTCCCGGCCCGGAGGGGATAGTGCGGAGTGTTTTCCAGACTT‑CGCGCTGAATGGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TTTCAATAATGCTAAGATGACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCcagcag                                                                                                                                                 <  2:635126/239‑1 (MQ=255)
                                                                                           gccggggcgggcgatttGCTGAGCCGTGTGGCCGACAT‑TCATTACCGGCCCGCAGGGGATGGGGCGTGGTGTTGTCCAGACGT‑CGCGCGGAAATGGCCTCCCCCCAGTGGAAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGAAAAATATTCAAGAAGCTTGCCGCTTAAACCGCCTGGATTGGTGGCAGAAAGGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCagcagc                                                                                                                                                <  2:635637/223‑1 (MQ=255)
                                                                                                    gCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGTGCTCCTAGCCTTCTTTTCGATAATGGATGTCCAGCACCTGCCCCCTGc                                                                                                                                      >  2:642113/1‑239 (MQ=255)
                                                                                                     cccAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATTGTGCGTAGGGTTTTCCAGACTT‑CGCGCTGTAATGGCGTCACCCCCGTAGCCGTGGGAAGCG‑TATCAATAATGCTCAGATTACCTGCAAAATATTCACGCAGCTTGTCCCTTAAACCGCCTGGATTGGTGGCAGAAAATCGCGCAATCCCTTCTTTGCGAAAATGGATggcaaacacgtgccccatacc                                                                                                                                     >  2:856764/1‑219 (MQ=255)
                                                                                                     cccAATTGCCCAGCCAGGTGGCCGTAAT‑GCATTCCCTGCCCGCAGGGGATAGTGCGTAGTTTTTTCCAGTCTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGc                                                                                                                                      <  1:561066/238‑1 (MQ=255)
                                                                                                         aTTGGTCAGCCAGTTGGCCGTAAT‑GCATTGCCGGCCCGCAGGGGTTAGTGCGTGGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCg                                                                                                                                 <  1:79861/239‑1 (MQ=255)
                                                                                                               cAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTCAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGt                                                                                                                             >  1:279183/1‑237 (MQ=255)
                                                                                                               cAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAAAATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATCCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCCGCACCATGCGTTCCCTgga                                                                                                                            >  2:604463/1‑236 (MQ=255)
                                                                                                                  ccAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTActc                                                                                                                         >  1:91656/1‑238 (MQ=255)
                                                                                                                       cggctctAAT‑GCATTCCCTGCCCGCAGGGCATAGGGCGTAGTGTTTTCCAGACTT‑CGCGTTGAAATGGCGTCCCCCCCGTAGCGGTGGGAAGCG‑TATCAATAATGCTAAGTTTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCa                                                                                                                   <  1:253734/233‑1 (MQ=255)
                                                                                                                           cctAAT‑TCCTTTCCCGCCCGGAGGGGGTAGTGGGTGTTTTGTTCCAGACTTGCGGCTTGAATGGGCGCCCCCCCCGT‑GGAGGGGGGAGCGTTACCAAAAAGTCTAAGATTCCGGGAAAAATATTCAGGAAGTTTGCCGTTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCAtt                                                                                                                 <  2:897295/236‑1 (MQ=255)
                                                                                                                                    tttcccTCCCCGGGGGGGATTGGGGGTTGTGTTTCCCAGTCTT‑CGCCGTGAAATGGCGTCCCCCCCGGAGCAGGGGGAAGCT‑TATCAATAATGCTAAGATCCCCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAAccgcccg                                                                                                       <  2:408567/235‑1 (MQ=255)
                                                                                                                                      tACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCa                                                                                                     >  1:62955/1‑239 (MQ=255)
                                                                                                                                           gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCCTGCGTTCGCTGTACTCTTCCCATTCAACCGACCCAAGGCGa                                                                                                >  2:747195/1‑239 (MQ=255)
                                                                                                                                            cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGaa                                                                                               >  1:277708/1‑239 (MQ=255)
                                                                                                                                                cAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAAtt                                                                                            >  1:590322/1‑238 (MQ=255)
                                                                                                                                                        gtgggtgtgtgtTTTCCAGACTT‑CGCGGTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCG‑TAGCAATAATGCTAAGAGTCCCGGAAAAATATTCGCGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGc                                                                                   <  2:377873/231‑1 (MQ=255)
                                                                                                                                                             ttAGTTTTTTCCGGTCTT‑CGCGTTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TAACAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAAt                                                                               <  2:379631/237‑1 (MQ=255)
                                                                                                                                                               aGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCAACGCAAATCa                                                                             >  2:12473/1‑238 (MQ=255)
                                                                                                                                                                 tgtttgccagtattcggcgtggAAATGGCTTCCCCCCCGTGTCGGGGGGAAGGG‑TATCAATAATGCTAATATTACCGGCAATATTTTAGCGAAGCGTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCAcc                                                                           <  2:279183/218‑1 (MQ=255)
                                                                                                                                                                  tttttCCAGCCTT‑CGGCCTGAAATGGCTCCCCCCCCGTAGCAGTGGGAAGCG‑TAACAATAATGCTAAGATTACCGGAAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCAccc                                                                          <  1:836007/237‑1 (MQ=255)
                                                                                                                                                                    tttCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCcagcag                                                                                                                                                 >  2:686751/1‑165 (MQ=255)
                                                                                                                                                                    tttCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCcagcag                                                                                                                                                 <  1:686751/165‑1 (MQ=255)
                                                                                                                                                                      tCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATTCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGt                                                                     >  2:650341/1‑239 (MQ=255)
                                                                                                                                                                      tCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGt                                                                     >  1:499322/1‑239 (MQ=255)
                                                                                                                                                                           aCTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCCCTGTACTCTTCCCATTAAACCGCCCGCAGGCGAAATTGCCCA‑CGCaaaaacaccaccagggcac                                                                >  2:833304/1‑222 (MQ=255)
                                                                                                                                                                                cgcgCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCAACAGTGGACCCAgt                                                            >  1:777822/1‑239 (MQ=255)
                                                                                                                                                                                                 ccccccGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGCAATTGCTCATCGCAAATACCCCACAGTGGACCCAGTGGCGGGGCAAatttttt                                           >  1:559283/1‑238 (MQ=255)
                                                                                                                                                                                                         aGCAGAGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGACTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGACCTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCCCAGTGGACCCAGTGGCGTGGCAAATTTttcttacaaaa                                   >  2:152441/1‑233 (MQ=255)
                                                                                                                                                                                                           cAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAACCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTAAACCGCCCCCAGGCGAAATTGCTCATCGCCAATCACCCACACTGGAACCAGTAGCGGGGTGATTTTTTTTTTAAGTaa                                  >  1:761539/1‑238 (MQ=255)
                                                                                                                                                                                                              tGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGGTTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCACATCACCCACAGTGGACCAAGGGGCGTGGCAATTTTTTTTTCaaaaaaaccc                              >  2:374145/1‑231 (MQ=255)
                                                                                                                                                                                                              tGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCCCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCCCATCGCAAATCACCCACAGTGGACCCCGTGGGGTGGCAAATTTTTTTTCaaaaacactc                              >  2:5669/1‑231 (MQ=255)
                                                                                                                                                                                                               gggAAGCT‑CTTCAATAATGTTAAGATTCCGGGAAAAATATTCCCGAAGCTTGGCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAAtctc                              <  2:78616/238‑1 (MQ=255)
                                                                                                                                                                                                               gggAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAAtctc                              >  1:771672/1‑238 (MQ=255)
                                                                                                                                                                                                                ggAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCa                             >  1:78616/1‑238 (MQ=255)
                                                                                                                                                                                                                ggAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGGGGCAAATTTTTCTTTAAGTAATCTCa                             >  1:563765/1‑238 (MQ=255)
                                                                                                                                                                                                                 gAAGCG‑TAAAAAAAATGCTAAGGTTACCGGCAAAATATTCCCGAAGCTTGCCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAg                            <  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)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
TCATGGTGCCGTTTCGGCCAATAACCCGATGGCAAGGTACGACGATGCTGATGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAAT‑GCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTT‑CGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCG‑TATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGCAGCTGCACCATGCGTTCGCTGTACTCTTCCCATTCAACCGCCCGCAGGCGAAATTGCTCATCGCAAATCACCCACAGTGGACCCAGTGGCGTGGCAATTTTTTCTTCAAGTAATCTCAGCATCCGTTCTCTCTTAAGACAAACGT  >  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: ‑