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→GCAtopADNA 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

CGCTTTCGTCAATCTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGT  >  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)
                            cTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACg                                                                                                                                                                                                                 >  1:3263/1‑251 (MQ=255)
                              tGACATTCGACCAAAATTCCGTCGTAATATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACATGAATTTTCGGGAACAGAGTCACGACAAGGGGTTGATATCCGAAGAGAGCGAGTCCATATCGGTAAATCGATGCCAGTGGAAGGTTTATAAACATGCGACGCATTCCTGGATGCATCACATTAGGTAAGGTGAATGTGGGTAAAGCATTTATCATCGTTTAGTCCCAGGCAAAAGCCAACACGAt                                                                                                                                                                                                               >  2:325930/1‑251 (MQ=255)
                                 cATTCGACCAAAATGCCGTCGTGCTATTGCGCCTGTAGTCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGAGATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATcaa                                                                                                                                                                                                            <  2:666010/251‑1 (MQ=255)
                                        ccAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGGGTCCCCGGCAAAAGCCAAAAAGATCAACAAGTAt                                                                                                                                                                                                     >  2:579009/1‑251 (MQ=255)
                                          aaaaTTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCt                                                                                                                                                                                                   <  2:3263/251‑1 (MQ=255)
                                                    cgtGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGAc                                                                                                                                                                                         <  1:606770/251‑1 (MQ=255)
                                                    cgtGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGAc                                                                                                                                                                                         <  2:138307/251‑1 (MQ=255)
                                                                       gCCAAGACCTGTTAACTCAGTAACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTTTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCg                                                                                                                                                                      <  1:382991/251‑1 (MQ=255)
                                                                                  ttAACTCAGTCAACTGAATTTTCGTGAACAGAGTCACGACAATGGGTTGATAGCCGCAGCGAGCGAGTCCATATCGGTAACTCGATGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGACGCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGTGTCCCCGGCAACAGCAAAAACGATCAACACGTATCTGGGTAGTGGCTACGCTGTGAAATCCAGCGTCGGTCACATc                                                                                                                                                           >  2:402786/1‑251 (MQ=255)
                                                                                              ccTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCg                                                                                                                                               >  1:630574/1‑251 (MQ=255)
                                                                                                tGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAg                                                                                                                                                                                                                                       <  1:678790/162‑1 (MQ=255)
                                                                                                   aTTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAg                                                                                                                                          <  1:85461/251‑1 (MQ=255)
                                                                                                          tGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAACTTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTc                                                                                                                                    >  1:399058/1‑250 (MQ=255)
                                                                                                              cAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTg                                                                                                                               <  1:310553/251‑1 (MQ=255)
                                                                                                               agagTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGc                                                                                                                              >  1:646131/1‑251 (MQ=255)
                                                                                                                       gACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGGTTATCCACGTGCGACGCATTCCTGGCAGCATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCCAAAGCCAAAACGATCAACAAGTATCTGTGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACCTCCGCGATTGTCAGCCCCGTGGCTCAGCTGCCCAAAagt                                                                                                                      >  2:289447/1‑250 (MQ=255)
                                                                                                                        aCAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTACCGACCAGTGGCTCAGCTGCCAAAAag                                                                                                                       >  1:438408/1‑249 (MQ=255)
                                                                                                                          aaGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGCCGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTg                                                                                                                   >  1:415885/1‑251 (MQ=255)
                                                                                                                           aGGGGATGATATCCGCAGAGAGCGCGTCCATATAGGTAACTCGTTGCCAGTGGAAGGATTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGGAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCACCGGCAACAGCCAAAACGATCAAAAAGTATCTGGGTAGTGACTACGTAGTGAAATCCAGCGTCGGTCACATCCGCGATTTGACGACCAGCGGCTCCGCTGCCAAAAAGAGTg                                                                                                                   >  2:240421/1‑250 (MQ=255)
                                                                                                                            ggggTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGcc                                                                                                                <  1:533860/251‑1 (MQ=255)
                                                                                                                                 tGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGAct                                                                                                            <  2:296809/250‑1 (MQ=255)
                                                                                                                                      tCCGCAGAGAGCGAGTCCATATCGGTAACTCGCTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTACTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGAAAAAGCCAAAAAGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCACCTGCCAAAAAGAGTGCCAACTCTAccc                                                                                                      >  2:73687/1‑250 (MQ=255)
                                                                                                                                       ccGCAGAGAGCGAGTCCATATCGGGAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAACTTCGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCCCCTCCGCGATTTGCCGACCAGTGGCTCAGCTGCCCAAAAGGTGGCCGCCTCTACCTc                                                                                                     >  2:61523/1‑251 (MQ=255)
                                                                                                                                         gaagagagCGAGTACATATAGGAAACTCGTTGCAAATGGAAGGATTATAAACGGGAGAAGCAGTCATGGAAGAATCAAATTAGGTAATGTGAAAATGGGTAAAGCATTTGTCATCGGTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGAAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGAGTGTCCGACCAGTGGCTCAGATGCCAAAAAGAGTGACGACTCTACCTcca                                                                                                   <  2:607461/249‑1 (MQ=255)
                                                                                                                                           agagagCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacc                                                                                                 >  1:455202/1‑251 (MQ=255)
                                                                                                                                             agagCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCGCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCaa                                                                                               <  1:543268/251‑1 (MQ=255)
                                                                                                                                              gagCGAGTCCATATCGGAAACTCGTTGCAAGTGGAAGGATTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTAACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAAGGCCGACTCTGCCTccacca                                                                                                >  1:271983/1‑249 (MQ=255)
                                                                                                                                                gCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGAc                                                                                            <  1:41755/251‑1 (MQ=255)
                                                                                                                                                  gAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg                                                                                          <  1:579009/251‑1 (MQ=255)
                                                                                                                                                  gAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg                                                                                          >  1:15791/1‑251 (MQ=255)
                                                                                                                                                  gAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGGGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCAGGCAAACGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTAACATCCGCGATTTGCCGACCAGTGGCTCAACTGCCAAAAAGAGTGCCGACTCTACCTCCCCCAAGACgg                                                                                          >  2:557670/1‑251 (MQ=255)
                                                                                                                                                     tCCATATCGATACCTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCCTTCCTGTAAGAATCACATTCGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTTAGTCCCCGGCAAAAGCCAAAACGATCACCAAGTATCTGGGTAGTGACTACGTGGTGCAATCCAGCGTCGAGCACATCCGCCATTTGCCGACCAGTGGCTCAGCTACCACACAGAGTGCCACCACTACCTCCACCAACACgcct                                                                                        >  2:444877/1‑247 (MQ=255)
                                                                                                                                                      ccATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTa                                                                                       <  1:54617/250‑1 (MQ=255)
                                                                                                                                                          atCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCCAGACGGCTaaaa                                                                                    >  2:199776/1‑249 (MQ=255)
                                                                                                                                                                 aCTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTaaa                                                                            >  1:238887/1‑250 (MQ=255)
                                                                                                                                                                      ttGCCAGTGGAAGGTTTATCAACGTGCTACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGAt                                                                      >  1:326039/1‑251 (MQ=255)
                                                                                                                                                                              ggAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTgg                                                              >  1:471654/1‑251 (MQ=255)
                                                                                                                                                                              ggAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCAATTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTg                                                               <  1:444877/250‑1 (MQ=255)
                                                                                                                                                                               gAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGAc                                                                                                                                                                                         >  1:514542/1‑128 (MQ=255)
                                                                                                                                                                               gAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGAc                                                                                                                                                                                         <  2:514542/128‑1 (MQ=255)
                                                                                                                                                                                         cAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTACGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCCAAAGCCAAAACGATCAACAAGTATCTGGGGAGTGACTACGTGGTGACATCCAGCGTCGGTCACACCCGCGCTTTGCCGACCAGTAGCTCATCTGCCAAAACGCTGGCCGACTCTACCGCCACCAAGACGGCTAAAAAGCCTAAAAAGGAAGAACGCAGCGCTCTCGTc                                                    >  2:539574/1‑250 (MQ=255)
                                                                                                                                                                                          aaCGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCa                                                   >  1:634692/1‑250 (MQ=255)
                                                                                                                                                                                            cGTGAGACGCATTACTGGAAGAATCAACTTAGGTAAGGGGAATATGGGTAAAGCATTTGTCAGCGTTGAGTCCCCGGCAAAAGACAAAACGATCAACAAGGATCTGGGTAGTGACTACGGGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAAcc                                                <  2:471654/251‑1 (MQ=255)
                                                                                                                                                                                                  aCGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATgg                                          <  1:234436/251‑1 (MQ=255)
                                                                                                                                                                                                   cGCATTCCTTGATGAATCATATAAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAATGCCAAAACGATCAACAAGTACCTGGGTAGTGACTACGAGGTGAAATCAAGCGTCGGTCACATCCGCGCTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGACGACTCTACCTCCACCAAGACGGACAAAAAGCCTAAAAAGGAGGAAcactgtcccctcgtaaaccggacggg                                         >  2:576488/1‑226 (MQ=255)
                                                                                                                                                                                                    gCATTTCTGGAAGAATCAAATTAAGGAAGGGGGAAATGGGTAAAGCATTTGTCAACCTTTGGTCCCCCGCAAAAGCCAAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATggg                                         <  1:234437/251‑1 (MQ=255)
                                                                                                                                                                                                         ccTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGATAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTc                                                    >  1:373707/1‑234 (MQ=255)
                                                                                                                                                                                                                       aTTAGGTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAAtt                     >  1:442975/1‑251 (MQ=255)
                                                                                                                                                                                                                          tggTAAGGTGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGCCCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGAGGCCAAATTggg                  >  2:88950/2‑251 (MQ=255)
                                                                                                                                                                                                                           gggaaGGGGAATATGTGTAAAGCATTTATCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGTTTGACTCGTGGCACAATTGGGa                 <  2:455202/248‑1 (MQ=255)
                                                                                                                                                                                                                                  tGAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCCCAATTGGGAGGCGCAc          >  1:321552/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggAATATGGGTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCAc          <  2:112604/250‑1 (MQ=255)
                                                                                                                                                                                                                                          ggTAAAGCGTTTGTCATCGTTGAGTCCCAGGCAAATGCCTAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCATCGTCGGGCACATCCGCGCTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGACGACTCTACCTCCACCAAGACGGCTAAATAGCCTGAAAAGGATGAACGTGGCGCTCTCGTCACCCGTCTGGGGGATGACCCGTGGCACAATTGGGGGGCGCACTATGa     >  2:57837/1‑248 (MQ=255)
                                                                                                                                                                                                                                          ggTAAAGCATTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAgg  >  1:470001/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                   
CGCTTTCGTCAATCTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGT  >  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: