Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,932,096 C→G 5.9% intergenic (+37/+19) purT → / ← eda phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,932,0960CG 5.9% 186.8 / 2.1 58intergenic (+37/+19)purT/edaphosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base G (2/2);  total (29/29)
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

TCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGCCGGGAAGAATTTGAACTCTTTCAAACCGTAGTCCATACCCAGCAT  >  NC_000913/1931803‑1932393
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                         
tCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCagg                                                                                                                                                                                                                                                                                                    >  2:140402/1‑301 (MQ=255)
    ccTTCCGTTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtt                                                                                                                                                                                                                                                                                                <  2:104559/300‑1 (MQ=255)
       tCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTg                                                                                                                                                                                                                                                                                             <  2:422846/301‑1 (MQ=255)
         cGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCCCGCCGCCGGACAGGTACAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATACCCGATCCTcaaaccggccattttac                                                                                                                                                                                                                                                                                           >  1:259198/1‑285 (MQ=255)
                      ccGTTAGTTTGTTCTTGCAGTTTTGGGTGTTATTTTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAg                                                                                                                                                                                                                                                                              <  1:613954/301‑1 (MQ=255)
                       cGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTCCAGGGTTAAACCCGGCAAAAAAACGCTACAAACATGCCCGAACCTCGATCGGGCATTTTGACTTTAACAGCTTAgc                                                                                                                                                                                                                                                                             >  2:503591/1‑301 (MQ=255)
                            gTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGCTTTGCCGATTCGTTTATTTGGTAAGCCGGCACTTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCCGAGAGTGTTGTTGCCGCCGTTGAACGCGCTACGCCCGCCGCCCGCCCGGTAAAAGTACAGGGCTAAACCCGGCAAAAAAACGCTCCAAAAATGCACGATCCTCGCTCGGCCATTCTACCTTTTACCGCCTAGAGCCtt                                                                                                                                                                                                                                                                        >  1:515318/1‑301 (MQ=255)
                                   ccTGCAGTTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTGTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTAGTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGGGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGc                                                                                                                                                                                                                                                                 <  2:179601/301‑1 (MQ=255)
                                             cggctgTTATTCTGCTACAACTTACCAGTCAGTATGTCATGTTTGTTGATTTGCAGATTGCCGTAGTGGCATATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGTGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGcc                                                                                                                                                                                                                                                       <  1:268759/296‑1 (MQ=255)
                                                tcGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGc                                                                                                                                                                                                                                                    <  2:457991/300‑1 (MQ=255)
                                                ccGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGc                                                                                                                                                                                                                                                    >  2:170168/1‑301 (MQ=255)
                                                 cttttttCTGCCACAACTTACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTGAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACGGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                                                                                                                                                                                                                                   <  1:169971/296‑1 (MQ=255)
                                                 cGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                                                                                                                                                                                                                                   >  2:343613/1‑301 (MQ=255)
                                                          gCCACCACTGACCCGTCCGAATGTCCCGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGCTTCGTTTATTTGGTCAGCCGGAAATTGCTGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTCCAGCTCGTGTTGTTGACCCCCTTGAACGCGCGCAGCACTCCGCCGGACAGGTAACCGTACCGGGTCACACCCGCCACAACCACGCTACCAAAATCCCCCATCATCGCCCGGGCATTTTCACTTTTACACCTTACCGCCTTCtcccgcatccctccccaaccccgtcattc                                                                                                                                                                                                                                          >  1:128994/1‑273 (MQ=255)
                                                                                tttaCGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGGCgattacgaccgcattact                                                                                                                                                                                                                    <  2:307113/298‑18 (MQ=255)
                                                                                   aCGTTTGTTAATGTGCAGAATGCCGTAGGCGTAGTTTTGCAGATTCGTTTATTTGTTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgcgc                                                                                                                                                                                                                 <  1:343150/301‑1 (MQ=255)
                                                                                           tAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCg                                                                                                                                                                                                         <  2:333933/301‑1 (MQ=255)
                                                                                           aaaggttagtagtgtggTAGGGGGATATTTTCTGTTTGGTTTATTTGGTAGGTCGGAAATTGTGGGCGGCCGTCGGCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCg                                                                                                                                                                                                         <  1:250036/285‑1 (MQ=255)
                                                                                                      aTGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagccag                                                                                                                                                                                              <  2:316192/301‑1 (MQ=255)
                                                                                                       tgtctttggCGCAGATTTGTAGATTTGTTTATTTGTTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAgg                                                                                                                                                                                             <  1:86386/294‑1 (MQ=255)
                                                                                                       ctctgTAGGCGCAGATTTGCATATTCGTTTATTTGGTAAGCCGGAAATTGTTGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAgg                                                                                                                                                                                             <  1:157367/297‑1 (MQ=255)
                                                                                                         ccGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGaa                                                                                                                                                                                           >  2:157518/1‑301 (MQ=255)
                                                                                                              gtgggtgttttgtgtgttCGTTTAGTGTGTAAGCCGGAATTTGGTGGTAGTCGTCGGCTGGGGGTGGCACTGGCTACTGAAGAGAGGGTTGTTGACGCAATGGACCGCGCGAAGCAGGCGGCCGGAAAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAGAAATGCCCGATCCTCGATCGGGCATTTTGACTTGTACAGCTTAGCGCCTGCTACAGCGTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAaccacc                                                                                                                                                                                      <  1:403954/285‑1 (MQ=255)
                                                                                                                                   ttATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAg                                                                                                                                                                 <  2:259466/301‑1 (MQ=255)
                                                                                                                                   ttATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAACCAGGAACCACCGACGCACACCCCCCTTTTCac                                                                                                                                                                 >  1:398139/1‑300 (MQ=255)
                                                                                                                                    tATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAgc                                                                                                                                                                >  2:402143/1‑301 (MQ=255)
                                                                                                                                    tATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAgc                                                                                                                                                                >  2:10264/1‑301 (MQ=255)
                                                                                                                                                cGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCa                                                                                                                                                    <  2:398653/301‑1 (MQ=255)
                                                                                                                                                    attcttggtcAGCTGTGGTCTGGGGTTGGCACTGGCTACTGCAGAGGGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGt                                                                                                                                                <  1:636895/292‑1 (MQ=255)
                                                                                                                                                     agtGAAGGCGGCCGTGGGGTGGGGGGGGCGTTGGTGACTGCGGGGAGTGTTGGTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTa                                                                                                                                               <  1:10264/299‑1 (MQ=255)
                                                                                                                                                         aTGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGATa                                                                                                                                           >  2:637772/1‑301 (MQ=255)
                                                                                                                                                             cAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCg                                                                                                                                       <  2:474263/301‑1 (MQ=255)
                                                                                                                                                                         gggggTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCACCAGGTAGTAACGGTAGTTAGCCGGAGAAATaacaa                                                                                                                           >  2:82084/1‑300 (MQ=255)
                                                                                                                                                                             gtagCACTGCCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGACGCACGCCGCCGGACAGGTACAAGTACAGGGGTAAACCCGGCAAAAAAACGCAACAAAACTGCTCGATCCACCATCGGGCATTTTGACTATTACAGCTTAGCGCCTTCTACCGCTTCCCGCGCCAGCTTAGTAATGCGGTAGTACTCGCCCGCTTCCAGCGCATCTGCCGGAACAAGCCAGGAACAACCGATGCACAGCACGCTTTTCAGCGCCAGGAAGTCAAGGTAGACAGCAGGAGAAATaccacaccc                                                                                                                       >  1:171153/4‑297 (MQ=255)
                                                                                                                                                                               ggCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACCCTTTTCCGCGCCAGGTAGTCACGGTAATTAGCCGGAGAACAACCCCCCGTCg                                                                                                                     >  1:577568/1‑301 (MQ=255)
                                                                                                                                                                                gCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCATAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGCAACCACCGATGCACAGCCCGCTTTTCAGCGCCAGGTAGTCACGGTACTTACCCGGCAAAATACCAccccacag                                                                                                                    >  1:538269/1‑296 (MQ=255)
                                                                                                                                                                                  acggggtACTGGAGAGAGGGTTGTTGACGCCATTGAACTCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGc                                                                                                                  <  1:249593/295‑1 (MQ=255)
                                                                                                                                                                                               gagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTCACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATTCGGTCGTAATCGCCCGCTTCCAGCCCCTCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAATCACGGTAATTAGCCGGAAAAAACACAACCCTAGGCCAGACACGGAccac                                                                                                     >  1:490386/1‑299 (MQ=255)
                                                                                                                                                                                                                   ccttttccgatctcACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTCCAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCg                                                                                 <  1:181796/288‑1 (MQ=255)
                                                                                                                                                                                                                     aaCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTCGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCCGGGAGAACGGACCCGCAAACGcc                                                                               >  1:281872/1‑301 (MQ=255)
                                                                                                                                                                                                                         cgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCCAACGCCCGcc                                                                           >  1:208064/1‑300 (MQ=255)
                                                                                                                                                                                                                           ggagttcagacgtgtgctcttccgatctACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAgg                                                                         <  2:629272/274‑1 (MQ=255)
                                                                                                                                                                                                                                cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCgagatcggaagagc                                                                    >  2:181993/1‑288 (MQ=255)
                                                                                                                                                                                                                                        cgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCAcgcagc                                                            >  2:612118/1‑301 (MQ=255)
                                                                                                                                                                                                                                          gACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACgccgccc                                                          >  2:549311/1‑300 (MQ=255)
                                                                                                                                                                                                                                            cAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCAGCGATCGCCTGCAGGGCTTTCACGCCGCCGtt                                                        >  2:403144/1‑301 (MQ=255)
                                                                                                                                                                                                                                                tAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGNNNNCGCCCGCTTCCNGCGCATCTGCCNNNNNCNNNNAGGAACNNNNNNNNNACAGCACNNNNTTCAGCGCCAGGTNNNNCNGGTAGNNNNNNNGAGAAATANNNNNCGTCGGGCAGAAACGGACCTGGGAGAACGAACACACGAACGCCCGCCGAGCTTTCCCACCGCCGtccacc                                                    >  1:81733/1‑296 (MQ=255)
                                                                                                                                                                                                                                                      tACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGcc                                              <  2:282175/301‑1 (MQ=255)
                                                                                                                                                                                                                                                       aCAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTCGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAggagatcggaaaagcggcgaggagggaaag                                             >  1:628395/1‑273 (MQ=255)
                                                                                                                                                                                                                                                         aGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACTGACCTGGGAGCACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTCCAGCCggg                                           >  1:570619/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         aGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCGTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTGTGCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGCCggg                                           <  2:171350/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGCACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCCGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGACACCGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCCTTCCCGCCGCCGTCAGCTTCACCCGGgaaggagt                                    >  1:437761/1‑297 (MQ=255)
                                                                                                                                                                                                                                                                  cccGGCAAAAAAACGCTACAAAAAGGCCCGATCCTCGATCGGGCATTTTGACGTTGACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCGTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGCAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGCCGGGAAGAATTTg                                  <  1:140277/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gggCAAAAAAAGGGTACAAAATTGCCCGTTCTTGGGTCGGGAATTGTGTCTGTTACAGCTTAGCGCCTTCTACAGCTGCACGAGCCAGCGTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGCCGGGAAGAATTTGaa                                <  1:401629/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGCCGGGAAGAATTTGAACTCTTTCAAAc                     <  2:578284/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        aaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGCCGGGAAGAATTTGAACTCTTTCAAACCGTAGTCCa            >  2:187892/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           gCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGCCGGGAAGAATTTGAACTCTTTCAAACCGTAGTCCATaa         >  2:174284/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                  tcTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGCCGGGAAGAATTTGAACTCTTTCAAACCGTAGTCCATACCCAGCAt  <  2:315583/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                  ccTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGCCGGGAAGAATTTGAACTCTTTCAAACCGTAGTCCATACCCAGCac  >  2:133007/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                         
TCGGACTTCCGGTTGGCGGGATCCGTCAGTATGGTCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAGTCACGGTAGTTAGCCGGAGAAATACCACCCGTCGGGCAGAAACGGACCTGGGAGAACGGACCCGCGATCGCCTGCAGGGCTTTCACGCCGCCGTTAGCTTCAGCCGGGAAGAATTTGAACTCTTTCAAACCGTAGTCCATACCCAGCAT  >  NC_000913/1931803‑1932393

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: