Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,331,059 2 bp→AT coding (12‑13/2598 nt) topA → DNA topoisomerase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,331,0590TA98.3% 188.7 / ‑6.5 60A4A (GCT→GCA) topADNA topoisomerase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (33/25);  minor base G (1/0);  total (35/25)
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
*NC_0009131,331,0600CT98.3% 189.7 / ‑6.8 60L5F (CTT→TTT) topADNA topoisomerase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/24);  minor base A (0/1);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGCTTTCGTCAATCTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGT  >  NC_000913/1330817‑1331301
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                   
cGCTTTCGTCAATCTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTTTCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCAt                                                                                                                                                                                                                                              <  2:293393/251‑1 (MQ=255)
 gCTTTCGTCAAACTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACATGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGCAGGTTTATCAACGTGCGACGCATTCCTGGTAGAATCAAATTAGGAATGGGGAATATGGGTAAAGCATTTGTCATc                                                                                                                                                                                                                                             >  2:499031/1‑251 (MQ=255)
     tCGTCAATCTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGtt                                                                                                                                                                                                                                          <  1:674093/250‑1 (MQ=255)
      cGTCAATCTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTg                                                                                                                                                                                                                                         >  1:182434/1‑250 (MQ=255)
          aaTCTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTTAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTcc                                                                                                                                                                                                                                    <  2:55642/251‑1 (MQ=255)
                       aTTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGACAGTGGAAGGTTTATCAACGTGCAACGCATTCCTGGAAGAATCAAATTAGGTAAGGGGCATATGGGTAAAGCCTTTGTCATCGTTGAGTCCCCGGCAAAAGcc                                                                                                                                                                                                                        >  2:606770/1‑250 (MQ=255)
                            cTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑Cg                                                                                                                                                                                                                 >  1:3263/1‑251 (MQ=255)
                              tGACATTCGACCAAAATTCCGTCGTAATATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACATGAATTTTCGGGAACAGAGTCACGACAAGGGGTTGATATCCGAAGAGAGCGAGTCCATATCGGTAAATCGATGCCAGTGGAAGGTTTATAAACATGCGACGCATTCCTGGATGCATCACATTAGGTAAGGTGAATGTGGGTAAAGCATTTATCATCGTTTAGTCCCAGGCAAAAGCCAACA‑CGAt                                                                                                                                                                                                               >  2:325930/1‑251 (MQ=255)
                                 cATTCGACCAAAATGCCGTCGTGCTATTGCGCCTGTAGTCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGAGATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATcaa                                                                                                                                                                                                            <  2:666010/251‑1 (MQ=255)
                                        ccAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGGGTCCCCGGCAAAAGCCAAAA‑AGATCAACAAGTAt                                                                                                                                                                                                     >  2:579009/1‑251 (MQ=255)
                                          aaaaTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCt                                                                                                                                                                                                   <  2:3263/251‑1 (MQ=255)
                                                    cgtGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGAc                                                                                                                                                                                         <  1:606770/251‑1 (MQ=255)
                                                    cgtGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGAc                                                                                                                                                                                         <  2:138307/251‑1 (MQ=255)
                                                                       gCCAAGACCTGTTAACTCAGTAACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTTTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCg                                                                                                                                                                      <  1:382991/251‑1 (MQ=255)
                                                                                  ttAACTCAGTCAACTGAATTTTCGTGAACAGAGTCACGACAATGGGTTGATAGCCGCAGCGAGCGAGTCCATATCGGTAACTCGATGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGACGCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGTGTCCCCGGCAACAGCAAAAA‑CGATCAACACGTATCTGGGTAGTGGCTACGCTGTGAAATCCAGCGTCGGTCACATc                                                                                                                                                           >  2:402786/1‑251 (MQ=255)
                                                                                              ccTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCg                                                                                                                                               >  1:630574/1‑251 (MQ=255)
                                                                                                tGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAg                                                                                                                                                                                                                                       <  1:678790/162‑1 (MQ=255)
                                                                                                   aTTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAg                                                                                                                                          <  1:85461/251‑1 (MQ=255)
                                                                                                          tGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAACTTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTc                                                                                                                                    >  1:399058/1‑250 (MQ=255)
                                                                                                              cAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTg                                                                                                                               <  1:310553/251‑1 (MQ=255)
                                                                                                               agagTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGc                                                                                                                              >  1:646131/1‑251 (MQ=255)
                                                                                                                       gACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGGTTATCCACGTGCGACGCATTCCTGGCAGCATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCCAAAGCCAAAA‑CGATCAACAAGTATCTGTGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACCTCCGCGATTGTCAGCCCCGTGGCTCAGCTGCCCAAAagt                                                                                                                      >  2:289447/1‑250 (MQ=255)
                                                                                                                        aCAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTACCGACCAGTGGCTCAGCTGCCAAAAag                                                                                                                       >  1:438408/1‑249 (MQ=255)
                                                                                                                          aaGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGCCGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTg                                                                                                                   >  1:415885/1‑251 (MQ=255)
                                                                                                                           aGGGGATGATATCCGCAGAGAGCGCGTCCATATAGGTAACTCGTTGCCAGTGGAAGGATTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGGAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCACCGGCAACAGCCAAAA‑CGATCAAAAAGTATCTGGGTAGTGACTACGTAGTGAAATCCAGCGTCGGTCACATCCGCGATTTGACGACCAGCGGCTCCGCTGCCAAAAAGAGTg                                                                                                                   >  2:240421/1‑250 (MQ=255)
                                                                                                                            ggggTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑cc                                                                                                                <  1:533860/251‑1 (MQ=255)
                                                                                                                                 tGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGAct                                                                                                            <  2:296809/250‑1 (MQ=255)
                                                                                                                                      tCCGCAGAGAGCGAGTCCATATCGGTAACTCGCTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTACTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGAAAAAGCCAAAA‑AGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCACCTGCCAAAAAGAGTG‑CCAACTCTAccc                                                                                                      >  2:73687/1‑250 (MQ=255)
                                                                                                                                       ccGCAGAGAGCGAGTCCATATCGGGAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAACTTCGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCCCCTCCGCGATTTGCCGACCAGTGGCTCAGCTGCCCAAAAG‑GTGGCCGCCTCTACCTc                                                                                                     >  2:61523/1‑251 (MQ=255)
                                                                                                                                         gaagagagCGAGTACATATAGGAAACTCGTTGCAAATGGAAGGATTATAAACGGGAGAAGCAGTCATGGAAGAATCAAATTAGGTAATGTGAAAATGGGTAAAGCATTTGTCATCGGTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGAAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGAGTGTCCGACCAGTGGCTCAGATGCCAAAAAGAGTG‑ACGACTCTACCTcca                                                                                                   <  2:607461/249‑1 (MQ=255)
                                                                                                                                           agagagCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTccacc                                                                                                 >  1:455202/1‑251 (MQ=255)
                                                                                                                                             agagCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCGCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCaa                                                                                               <  1:543268/251‑1 (MQ=255)
                                                                                                                                              gagCGAGTCCATATCGGAAACTCGTTGCAAGTGGAAGGATTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTAACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAAGG‑CCGACTCTGCCTccacca                                                                                                >  1:271983/1‑249 (MQ=255)
                                                                                                                                                gCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGAc                                                                                            <  1:41755/251‑1 (MQ=255)
                                                                                                                                                  gAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACgg                                                                                          <  1:579009/251‑1 (MQ=255)
                                                                                                                                                  gAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACgg                                                                                          >  1:15791/1‑251 (MQ=255)
                                                                                                                                                  gAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGGGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCAGGCAAACGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTAACATCCGCGATTTGCCGACCAGTGGCTCAACTGCCAAAAAGAGTG‑CCGACTCTACCTCCCCCAAGACgg                                                                                          >  2:557670/1‑251 (MQ=255)
                                                                                                                                                     tCCATATCGATACCTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCCTTCCTGTAAGAATCACATTCGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTTAGTCCCCGGCAAAAGCCAAAA‑CGATCACCAAGTATCTGGGTAGTGACTACGTGGTGCAATCCAGCGTCGAGCACATCCGCCATTTGCCGACCAGTGGCTCAGCTACCACACAGAGTG‑CCACCACTACCTCCACCAACACgcct                                                                                        >  2:444877/1‑247 (MQ=255)
                                                                                                                                                      ccATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTa                                                                                       <  1:54617/250‑1 (MQ=255)
                                                                                                                                                          atCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCCAGACGGCTaaaa                                                                                    >  2:199776/1‑249 (MQ=255)
                                                                                                                                                                 aCTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTaaa                                                                            >  1:238887/1‑250 (MQ=255)
                                                                                                                                                                      ttGCCAGTGGAAGGTTTATCAACGTGCTACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGAt                                                                      >  1:326039/1‑251 (MQ=255)
                                                                                                                                                                              ggAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTgg                                                              >  1:471654/1‑251 (MQ=255)
                                                                                                                                                                              ggAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCAATTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTg                                                               <  1:444877/250‑1 (MQ=255)
                                                                                                                                                                               gAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGAc                                                                                                                                                                                         >  1:514542/1‑128 (MQ=255)
                                                                                                                                                                               gAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGAc                                                                                                                                                                                         <  2:514542/128‑1 (MQ=255)
                                                                                                                                                                                         cAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTACGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCCAAAGCCAAAA‑CGATCAACAAGTATCTGGGGAGTGACTACGTGGTGACATCCAGCGTCGGTCACACCCGCGCTTTGCCGACCAGTAGCTCATCTGCCAAAACGCTGG‑CCGACTCTACCGCCACCAAGACGGCTAAAAAGCCTAAAAAGGAAGAACGCAGCGCTCTCGTc                                                    >  2:539574/1‑250 (MQ=255)
                                                                                                                                                                                          aaCGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCa                                                   >  1:634692/1‑250 (MQ=255)
                                                                                                                                                                                            cGTGAGACGCATTACTGGAAGAATCAACTTAGGTAAGGGGAATATGGGTAAAGCATTTGTCAGCGTTGAGTCCCCGGCAAAAGACAAAA‑CGATCAACAAGGATCTGGGTAGTGACTACGGGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAAcc                                                <  2:471654/251‑1 (MQ=255)
                                                                                                                                                                                                  aCGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATgg                                          <  1:234436/251‑1 (MQ=255)
                                                                                                                                                                                                   cGCATTCCTTGATGAATCATATAAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAATGCCAAAA‑CGATCAACAAGTACCTGGGTAGTGACTACGAGGTGAAATCAAGCGTCGGTCACATCCGCGCTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑ACGACTCTACCTCCACCAAGACGGACAAAAAGCCTAAAAAGGAGGAAcactgtcccctcgtaaaccggacggg                                         >  2:576488/1‑226 (MQ=255)
                                                                                                                                                                                                    gCATTTCTGGAAGAATCAAATTAAGGAAGGGGGAAATGGGTAAAGCATTTGTCAACCTTTGGTCCCCCGCAAAAGCCAAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATggg                                         <  1:234437/251‑1 (MQ=255)
                                                                                                                                                                                                         ccTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGATAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTc                                                    >  1:373707/1‑234 (MQ=255)
                                                                                                                                                                                                                       aTTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAAtt                     >  1:442975/1‑251 (MQ=255)
                                                                                                                                                                                                                          tggTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGCCCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGAGGCCAAATTggg                  >  2:88950/2‑251 (MQ=255)
                                                                                                                                                                                                                           gggaaGGGGAATATGTGTAAAGCATTTATCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGTTTGACTCGTGGCACAATTGGGa                 <  2:455202/248‑1 (MQ=255)
                                                                                                                                                                                                                                  tGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCCCAATTGGGAGGCGCAc          >  1:321552/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCAc          <  2:112604/250‑1 (MQ=255)
                                                                                                                                                                                                                                          ggTAAAGCGTTTGTCATCGTTGAGTCCCAGGCAAATGCCTAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCATCGTCGGGCACATCCGCGCTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑ACGACTCTACCTCCACCAAGACGGCTAAATAGCCTGAAAAGGATGAACGTGGCGCTCTCGTCACCCGTCTGGGGGATGACCCGTGGCACAATTGGGGGGCGCACTATGa     >  2:57837/1‑248 (MQ=255)
                                                                                                                                                                                                                                          ggTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAgg  >  1:470001/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                   
CGCTTTCGTCAATCTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAA‑CGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTG‑CCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGT  >  NC_000913/1330817‑1331301

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