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.2% 197.7 / ‑3.9 56coding (259/516 nt)ogtO‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (26/29);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATA  >  NC_000913/1399841‑1400101
                                                                                                                                           |                                                                                                                           
gggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTc                                                                                                                            <  1:53941/139‑1 (MQ=255)
 ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc                                                                                                                           <  1:1066367/139‑1 (MQ=255)
 ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc                                                                                                                           <  1:315453/139‑1 (MQ=255)
 ggATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcc                                                                                                                           <  2:562957/139‑1 (MQ=255)
   aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcccc                                                                                                                         <  1:479779/139‑1 (MQ=255)
   aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcccc                                                                                                                         <  1:775662/139‑1 (MQ=255)
   aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTcccc                                                                                                                         <  1:1574112/139‑1 (MQ=255)
        atacGTTTGCCACACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCGGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTccccccccg                                                                                                                    >  1:934502/4‑137 (MQ=255)
        aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTccccccccg                                                                                                                    <  1:420363/139‑3 (MQ=255)
        aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTccccccccg                                                                                                                    <  1:986539/139‑3 (MQ=255)
          ccGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTccccccccgta                                                                                                                  <  1:929865/139‑5 (MQ=255)
             tttGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCa                                                                                                                >  2:520141/1‑139 (MQ=255)
             tttGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCa                                                                                                                >  1:1188706/1‑139 (MQ=255)
                  cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                           <  2:934502/139‑1 (MQ=255)
                  cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                           <  1:283287/139‑1 (MQ=255)
                   gCACCAACGGTACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                          <  1:726847/139‑1 (MQ=255)
                     aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                        <  2:1021577/139‑1 (MQ=255)
                             cACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                <  1:94860/139‑1 (MQ=255)
                                  gccgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATg                                                                                           <  1:520141/139‑1 (MQ=255)
                                                 cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                          <  2:724574/109‑1 (MQ=255)
                                                 cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                          >  1:724574/1‑109 (MQ=255)
                                                 cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaat                                                                                               >  2:437862/1‑120 (MQ=255)
                                                 cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaat                                                                                               <  1:437862/120‑1 (MQ=255)
                                                 cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATg                                                                                           <  1:1080954/124‑1 (MQ=255)
                                                 cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                            >  1:1140440/1‑139 (MQ=255)
                                                 cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                            >  1:213188/1‑139 (MQ=255)
                                                 cccAATTGCACAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATg                                                                                           >  2:1080954/1‑124 (MQ=255)
                                                     aTTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAat                                                                        <  2:1258019/139‑1 (MQ=255)
                                                         cTCAGCCAGTTGGCCGTAATGCATTACCTTCGCGAAGTGGGTAGTGCGTAGTGTTTTCCAGTATTCGCGCTGAATTCACGTTCGCCCTCGTAGAAGTGGGAAGCGCATCAATAATGGTGAGATTTCCGGCAAAATATTc                                                                    <  1:92454/139‑1 (MQ=255)
                                                          tCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCa                                                                   <  2:1188706/139‑1 (MQ=255)
                                                            aGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCa                                                                                                                <  2:136712/92‑1 (MQ=255)
                                                            aGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCa                                                                                                                >  1:136712/1‑92 (MQ=255)
                                                                   tGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                <  2:581358/101‑1 (MQ=255)
                                                                   tGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                >  1:581358/1‑101 (MQ=255)
                                                                   tGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTg                                                          >  1:39090/1‑139 (MQ=255)
                                                                      ccGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCg                                                       >  2:446156/1‑139 (MQ=255)
                                                                           aTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                                  <  1:446156/139‑1 (MQ=255)
                                                                           aTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                                  <  1:238652/139‑1 (MQ=255)
                                                                           aTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                                  <  2:1140440/139‑1 (MQ=255)
                                                                                ttACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGc                                             >  1:50639/1‑139 (MQ=255)
                                                                                   ccTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                           <  2:587935/74‑1 (MQ=255)
                                                                                   ccTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTggg                                                                                                           >  1:587935/1‑74 (MQ=255)
                                                                                     tGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                        <  1:811447/139‑1 (MQ=255)
                                                                                       cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTc                                                                                                                            <  1:84899/52‑1 (MQ=255)
                                                                                       cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTc                                                                                                                            >  2:84899/1‑52 (MQ=255)
                                                                                           cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                        <  1:30790/69‑1 (MQ=255)
                                                                                           cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                        >  2:30790/1‑69 (MQ=255)
                                                                                           cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                            <  2:371115/97‑1 (MQ=255)
                                                                                           cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                            >  1:371115/1‑97 (MQ=255)
                                                                                                gATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTccccccccg                                                                                                                    >  2:1155533/1‑49 (MQ=255)
                                                                                                gATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTccccccccg                                                                                                                    <  1:1155533/51‑3 (MQ=255)
                                                                                                    gTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                         <  2:1426645/139‑1 (MQ=255)
                                                                                                     tgagTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCt                                           <  2:947902/117‑1 (MQ=255)
                                                                                                     tGCGTAGTGTTTTCCAGACTACGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCt                                           >  1:947902/1‑120 (MQ=255)
                                                                                                      gCGTAGTGTTTTCCAGACTTCGCTCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcg                       >  2:1210088/1‑139 (MQ=255)
                                                                                                       cGTAGTGTTTTCCAGACTTCGCTCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                      <  1:1210088/139‑1 (MQ=255)
                                                                                                       cGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                      <  2:1448271/139‑1 (MQ=255)
                                                                                                        gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                         <  2:122343/135‑1 (MQ=255)
                                                                                                        gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                         >  1:122343/1‑135 (MQ=255)
                                                                                                        gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTccccccccg                                                                                                                    >  2:556785/1‑41 (MQ=255)
                                                                                                        gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTccccccccg                                                                                                                    <  1:556785/43‑3 (MQ=255)
                                                                                                          aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                 <  2:1277896/45‑1 (MQ=255)
                                                                                                          aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                 >  1:1277896/1‑45 (MQ=255)
                                                                                                          aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                 <  2:722352/45‑1 (MQ=255)
                                                                                                          aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                 >  1:722352/1‑45 (MQ=255)
                                                                                                             gtTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAg                >  2:783517/1‑139 (MQ=255)
                                                                                                              ttttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaataat                                                                                            >  2:392582/1‑62 (MQ=255)
                                                                                                              ttttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaataat                                                                                            <  1:392582/62‑1 (MQ=255)
                                                                                                               tttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGcc              >  2:4517/1‑139 (MQ=255)
                                                                                                                 tCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCTTAGCAGTGGGAAGCGTATCAATAATACTAAGATTACCTGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCctt            <  1:187077/139‑1 (MQ=255)
                                                                                                                      aCTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTg       >  2:592074/1‑139 (MQ=255)
                                                                                                                          cgcgCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTAt                                                                                                   >  1:1613188/1‑43 (MQ=255)
                                                                                                                          cgcgCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTAt                                                                                                   <  2:1613188/43‑1 (MQ=255)
                                                                                                                           gcgcTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATa  >  2:1166490/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                           
GGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATA  >  NC_000913/1399841‑1400101

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: