Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,399,985 (C)7→8 coding (252/516 nt) 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,399,9781.C98.0% 177.9 / ‑3.9 51coding (259/516 nt)ogtO‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (19/31);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTG  >  NC_000913/1399842‑1400096
                                                                                                                                         |                                                                                                                        
ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                       <  1:488002/139‑1 (MQ=255)
ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                       <  1:1447616/139‑1 (MQ=255)
 gATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑ccc                                                                                                                      <  1:1626128/139‑1 (MQ=255)
  aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCGCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGC‑‑cccc                                                                                                                     <  2:1121610/139‑1 (MQ=255)
    tCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                       <  2:31854/135‑1 (MQ=255)
    tCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                       >  1:31854/1‑135 (MQ=255)
       aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc‑ccccg                                                                                                               <  2:1102364/139‑3 (MQ=255)
       aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc‑ccccg                                                                                                               <  2:1441482/139‑3 (MQ=255)
       aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc‑ccccg                                                                                                               <  1:1551861/139‑3 (MQ=255)
          cGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc‑ccccgtag                                                                                                            >  2:781947/1‑134 (MQ=255)
           gTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGc                                                                                                            <  1:519992/139‑1 (MQ=255)
                 cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTggg                                                                                                      <  1:383099/139‑1 (MQ=255)
                  gCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGa                                                                                                     <  1:1068036/139‑1 (MQ=255)
                  gCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGa                                                                                                     <  1:1447910/139‑1 (MQ=255)
                    aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAg                                                                                                   <  1:771443/139‑1 (MQ=255)
                    aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAg                                                                                                   <  1:452867/139‑1 (MQ=255)
                       aaCGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGt                                                                                                <  2:1632037/139‑1 (MQ=255)
                         cGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTAt                                                                                              <  2:722823/139‑1 (MQ=255)
                         cGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                       >  1:1684920/1‑114 (MQ=255)
                         cGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                       <  2:1684920/114‑1 (MQ=255)
                            cACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCaa                                                                                           <  1:781947/139‑1 (MQ=255)
                            cACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCaa                                                                                           <  2:133903/139‑1 (MQ=255)
                            cACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCaa                                                                                           <  1:932195/139‑1 (MQ=255)
                             aCGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCaat                                                                                          <  2:717810/139‑1 (MQ=255)
                              cgcgccgcgcGAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCATAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc‑ccccg                                                                                                               >  1:1259318/1‑114 (MQ=255)
                              cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc‑ccccg                                                                                                               <  2:1259318/116‑3 (MQ=255)
                                   cgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCt                                                                                    <  1:1102476/139‑1 (MQ=255)
                                    gcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTa                                                                                   >  2:175906/1‑139 (MQ=255)
                                      gccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGt                                                                                                >  2:447643/1‑124 (MQ=255)
                                      gccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGt                                                                                                <  1:447643/124‑1 (MQ=255)
                                               gCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGc                                                                        <  2:432449/139‑1 (MQ=255)
                                                cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                       >  1:888035/1‑139 (MQ=255)
                                                cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc                                                                                                                     >  2:1826031/1‑93 (MQ=255)
                                                cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc                                                                                                                     <  1:1826031/93‑1 (MQ=255)
                                                        cTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTc                                                               >  1:1054933/1‑139 (MQ=255)
                                                          cAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc                                                                                                                     <  1:1165022/83‑1 (MQ=255)
                                                          cAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc                                                                                                                     >  2:1165022/1‑83 (MQ=255)
                                                                          aTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc                                                                                                                     <  1:1041206/67‑1 (MQ=255)
                                                                          aTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc                                                                                                                     >  2:1041206/1‑67 (MQ=255)
                                                                           tGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaaa                                            >  2:1729555/1‑139 (MQ=255)
                                                                                tACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                       >  1:1076448/1‑59 (MQ=255)
                                                                                tACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                       <  2:1076448/59‑1 (MQ=255)
                                                                                 aCCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCt                                      >  1:1157674/1‑139 (MQ=255)
                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCaa                                                                                           >  1:756067/1‑82 (MQ=255)
                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCaa                                                                                           <  2:756067/82‑1 (MQ=255)
                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCaa                                                                                           >  2:63675/1‑82 (MQ=255)
                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCaa                                                                                           <  1:63675/82‑1 (MQ=255)
                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGAt                                  >  2:320132/1‑139 (MQ=255)
                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑c                                                                                                                        >  1:736064/1‑53 (MQ=255)
                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑c                                                                                                                        <  2:736064/53‑1 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCa                                                                                                           <  2:978002/64‑1 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCa                                                                                                           >  1:978002/1‑64 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATg                                                                                      <  2:263568/85‑1 (MQ=255)
                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATg                                                                                      >  1:263568/1‑85 (MQ=255)
                                                                                        cGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGa                                                                                                     <  1:220718/69‑1 (MQ=255)
                                                                                        cGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGa                                                                                                     >  2:220718/1‑69 (MQ=255)
                                                                                            ggggATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc                                                                                                                     >  1:1848001/1‑49 (MQ=255)
                                                                                            ggggATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc                                                                                                                     <  2:1848001/49‑1 (MQ=255)
                                                                                                   gTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                    <  1:642330/139‑1 (MQ=255)
                                                                                                       gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                       <  2:638450/36‑1 (MQ=255)
                                                                                                       gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                       >  1:638450/1‑36 (MQ=255)
                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAAAAATGCGCTc               <  2:581195/139‑1 (MQ=255)
                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTggg                                                                                                      <  1:1064566/51‑1 (MQ=255)
                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTggg                                                                                                      >  2:1064566/1‑51 (MQ=255)
                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGt                                                                                                >  2:1337504/1‑57 (MQ=255)
                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGt                                                                                                <  1:1337504/57‑1 (MQ=255)
                                                                                                         aGTGTTTCCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTCGCAGTGGGAAGCGTATCAATAATGCTAACATTACCGGCAACATATGCACGAAGCTTGTCGCTTAAAGCGCCTGCTTTGGTGGCAGACATGCACTCa              >  1:1211659/1‑139 (MQ=255)
                                                                                                             ttttCCAGACTTCGCGCTGAAGTGGCGG‑‑CCCCTCCCGTAACAGTTTGACGCGTATCAATAATGCTAAGATTACCGGCTAAATATTCACGAAACTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAGTGCGCTCATAGc          >  1:1381522/1‑139 (MQ=255)
                                                                                                              tttCCAGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGcc         >  1:950856/1‑139 (MQ=255)
                                                                                                                   aGACTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCcttctt    <  2:950856/139‑1 (MQ=255)
                                                                                                                     aCTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                      >  2:1223618/1‑87 (MQ=255)
                                                                                                                     aCTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                      <  1:1223618/87‑1 (MQ=255)
                                                                                                                     aCTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTg  <  1:1798960/139‑1 (MQ=255)
                                                                                                                     aCTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTg  <  1:434970/139‑1 (MQ=255)
                                                                                                                     aCTTCGCGCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTg  <  1:175906/139‑1 (MQ=255)
                                                                                                                         cgcgCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGa                                                                                                     <  1:758980/36‑1 (MQ=255)
                                                                                                                         cgcgCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGa                                                                                                     >  2:758980/1‑36 (MQ=255)
                                                                                                                         cgcgCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCaat                                                                                          <  2:409369/47‑1 (MQ=255)
                                                                                                                         cgcgCTGAAATGGCGTC‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCaat                                                                                          >  1:409369/1‑47 (MQ=255)
                                                                                                                                         |                                                                                                                        
GGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑CCCC‑CCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTG  >  NC_000913/1399842‑1400096

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑