Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 718,669 C→T 14.0% E598K (GAA→AAA)  speF ← ornithine decarboxylase isozyme, inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913718,6690CT14.0% 214.8 / 6.6 76E598K (GAA→AAA) speFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base C (9/55);  new base T (9/2);  total (18/58)
Fisher's exact test for biased strand distribution p-value = 1.64e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCACCCCAGATTTCACCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAAGATCGCATTTTTCCGGCACTACGCCATTTTCGCGCAGGAAGTTAGCAAGAATC  >  NC_000913/718436‑718907
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               
accaccCCAGATTTCACCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGac                                                                                                                                                                                                                                            <  2:189336/239‑1 (MQ=255)
    ccccAGATTTCACCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTAAAAAACATTCATTTCCTggacca                                                                                                                                                                                                                                         >  1:718812/1‑233 (MQ=255)
    ccccAGATTTCACCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAATCATGCATTTCCTTTcaca                                                                                                                                                                                                                                         >  1:243295/1‑238 (MQ=255)
          aTTTCACCCGGAACAACACACGCCCCACCCGGAGGATAAGGAAGCTCACCTTCGGCAGCGATACGGCCTGCGGCATCCTGCCGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACg                                                                                                                                                                                                                                 <  2:52428/239‑1 (MQ=255)
           tttCACCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACg                                                                                                                                                                                                                                 <  1:992026/238‑1 (MQ=255)
           tttCACCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACg                                                                                                                                                                                                                                 <  1:412067/238‑1 (MQ=255)
           tttCACCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACg                                                                                                                                                                                                                                 <  1:404436/238‑1 (MQ=255)
           tttCACCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGc                                                                                                                                                                                                                                <  2:6103/239‑1 (MQ=255)
           tttCACCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGc                                                                                                                                                                                                                                <  2:1026482/239‑1 (MQ=255)
           tttCACCCCGAACAACACACAGCACACCCGGAGGATAAGGAAGCGAACCGTCGGCAGCGATACGGCCTTCTGCATCCGGCAGGCGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCGTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACg                                                                                                                                                                                                                                 <  1:1037977/238‑1 (MQ=255)
           ttccACCCGGACCAACACCCAGAACCCCCCGAGGATAAGGCAGCGCACCTTCGGAAGCCATACGGCCTTGTGCATCCGGCAGACGAACCAGTCCCCCTTCACCCCTTAAATATTCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCTTTACGGAACATCTCTTTTTTCTTTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACACCTGACg                                                                                                                                                                                                                                 <  2:787388/235‑1 (MQ=255)
           tgtCACCCGGAGCAACACACAGCACACCCGGAGGATAAGGAGGCGCAACTTCGGCAGCGATACGGCCGTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGAGTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACg                                                                                                                                                                                                                                 <  2:950437/236‑1 (MQ=255)
           tgtCACCCGGAACATCACACAGCACGCCCGGAGGATAAGGAAGCGCCCCCTCGGCATCGATACGGCCTTCTGCATCCGGCATACTAACCTGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGc                                                                                                                                                                                                                                <  2:725590/237‑1 (MQ=255)
               aCCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGTATACAAAACACGCATTTTCTTGCACAAccttcgccggg                                                                                                                                                                                                                             >  1:352910/1‑229 (MQ=255)
               aCCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATACCTTTTCTGTCACAAATGTCTCAggg                                                                                                                                                                                                                            >  1:25216/1‑239 (MQ=255)
               aCCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAAATCATTTTCTCACCCAAccgtaacacgg                                                                                                                                                                                                                             >  1:420180/1‑229 (MQ=255)
                    gAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCAATAAAAACAAGCATTTCCTCTGTCAACAGACGTAGGGtgaaa                                                                                                                                                                                                                       >  2:411996/1‑235 (MQ=255)
                          cacacaGCACACCCGGAGGATAAGGAAGCGCCCCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCg                                                                                                                                                                                                                 <  2:1043125/239‑1 (MQ=255)
                          cacacaGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCg                                                                                                                                                                                                                 <  1:297561/239‑1 (MQ=255)
                              caGCACACCCGGAAGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAgc                                                                                                                                                                                                             <  1:633203/239‑1 (MQ=255)
                                gCACACCCTGAGGATAAGGAAGCGCCCCTGCGGCATCGATTCGTCCTTCGGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAgcgc                                                                                                                                                                                                           <  2:1189640/239‑1 (MQ=255)
                                  acacCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAgcgc                                                                                                                                                                                                           <  2:71931/237‑1 (MQ=255)
                                  acacCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAgcg                                                                                                                                                                                                            <  1:854637/236‑1 (MQ=255)
                                  acacCCGGAGGATAAGGAAGCGCACCTGCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAgcg                                                                                                                                                                                                            <  1:624581/236‑1 (MQ=255)
                                    acCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCTTTTTCTTACACAACCGACGCgggggggaaacggccgaccgcgctt                                                                                                                                                                                                       >  2:149003/1‑218 (MQ=255)
                                        ggatgaTAAGGAAGCGCCACGTCGGCAGCCGTACGGCCGGTTGCATCCGTCAGCCGAGCCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGGTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCACGTTTTGCAGTTGTGTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGTCACAATTGACGCAGGGTATAACCGGCGGAGTGCTCTTCATg                                                                                                                                                                                                   <  2:462595/239‑1 (MQ=255)
                                        ggaggaTAAGGAAGCGCATCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACGTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACCCAGGGTATAACCGGCGTAGCGCTCTTCATg                                                                                                                                                                                                   <  1:876006/239‑1 (MQ=255)
                                        ggaggaTAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATg                                                                                                                                                                                                   <  2:48591/239‑1 (MQ=255)
                                        ggaggaTAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATg                                                                                                                                                                                                   <  1:1191326/239‑1 (MQ=255)
                                         gaggaTAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAATAATTCATTTTCTGACACAAAAGAGGCAGGGGATAACCGGCGTAACGCGCTTcct                                                                                                                                                                                                    >  1:137979/1‑235 (MQ=255)
                                         gaggaTAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCTTTTCATGGCACAACTGACGCTGGGGGTAACCGGCGGAGCGCTCTTCAtc                                                                                                                                                                                                   >  2:396475/1‑237 (MQ=255)
                                                      cgcACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATGCATTTTCTGCGACAACAGCAGCAGGGGGTAACCGCCGGAGCGCTCTTCacgtcttttgttcaggg                                                                                                                                                                                     >  2:970700/1‑223 (MQ=255)
                                                       gcACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGa                                                                                                                                                                                    <  1:363376/239‑1 (MQ=255)
                                                       gcACCGTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATg                                                                                                                                                                                      <  1:586723/237‑1 (MQ=255)
                                                        cACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATATGCGTAGCTGGCTTCTTGCGGATTCATGCTGACGCGGGGGAAGTGCTCCTTACGGCACATCTCTTTTTGCAGTTGTTTCACGTTGCGGCGGGCGTAGAAATCCTGCATTTTTTCACACAAACCACTCCGGGTGTGACAGCCGTAGCGccctctatatccctcgattatagca                                                                                                                                                                                   >  2:698280/1‑218 (MQ=255)
                                                            ttCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAACTGCTTTTTCTGTCGCAACACACGCCGGGGGTGAACGGCGTAGCGCTGTTCTTGCTGTTTTTAGTTGGAgggga                                                                                                                                                                               >  2:312179/1‑237 (MQ=255)
                                                             tCGGCAGCGATACGGCCTTCTGCATCGGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAg                                                                                                                                                                              <  1:633229/239‑1 (MQ=255)
                                                             tCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAg                                                                                                                                                                              <  2:509223/239‑1 (MQ=255)
                                                             tCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAg                                                                                                                                                                              <  1:54950/239‑1 (MQ=255)
                                                              cGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGc                                                                                                                                                                             <  1:312179/239‑1 (MQ=255)
                                                               ggCAGCGATACGGCCTTCTTCAGCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCa                                                                                                                                                                            <  1:115021/239‑1 (MQ=255)
                                                                  aGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACt                                                                                                                                                                          <  2:1027538/238‑1 (MQ=255)
                                                                    cGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCg                                                                                                                                                                                                                 <  1:720175/197‑1 (MQ=255)
                                                                     gATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAATCAAGCATTTCCTTGCCCAAACGAACCAGGGGATAACCCGCCGAGCGCTCTTCATTCTTTTTTTAGATTGAAGGTAGCACTTCt                                                                                                                                                                       >  1:804585/1‑238 (MQ=255)
                                                                     gATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGAGATAAAAAACATACATTTTCTTTAACAAATCACGGCGGGGTAAAACCGCGTTACGCCTCTCATTCTGTTTTTATATgggaaggagaccttg                                                                                                                                                                        >  1:210789/1‑224 (MQ=255)
                                                                                 tGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGt                                                                                                                                                          <  1:521415/239‑1 (MQ=255)
                                                                                              cGAATCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTACTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTGGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAg                                                                                                                                                <  1:999210/236‑1 (MQ=255)
                                                                                              cGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCa                                                                                                                                              <  1:269007/238‑1 (MQ=255)
                                                                                                       tCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGCTTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGa                                                                                                                                     <  2:1144809/238‑1 (MQ=255)
                                                                                                       tCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGa                                                                                                                                     <  1:798728/238‑1 (MQ=255)
                                                                                                                 cgcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACcagc                                                                                                                          <  2:1198503/239‑1 (MQ=255)
                                                                                                                          aGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCGTGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca                                                                                                                 <  2:324776/239‑1 (MQ=255)
                                                                                                                                 gTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCt                                                                                                          <  1:764351/239‑1 (MQ=255)
                                                                                                                                       ttcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGt                                                                                                     <  1:1131588/238‑1 (MQ=255)
                                                                                                                                       ttcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGGTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGt                                                                                                     <  2:701316/238‑1 (MQ=255)
                                                                                                                                       gtcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGt                                                                                                     <  1:1032484/237‑1 (MQ=255)
                                                                                                                                                aTTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCa                                                                                                                                              <  1:1178000/188‑1 (MQ=255)
                                                                                                                                                            gcgcGGGAAGTGCTCGTGACGGAACATCTCGTGTTGCAGGTGTTGCAAAGTGTGGCGGGCAGACAAATTCTGCATGTGCTGACACAACTAACGCCGGGTATACCCGGCGTGGCGCTCTTCATGCTGTTTGTAGATGGAAGGGAGCATTTCTGCTAAGGGCGCGTCGGACACAAGCAGGTTTGCGAAGCGTACGGGGAGGGCAACAAGTTGCTGAAGTTGGGCCATAGCTTCTGCCGGAg                                                                               <  2:138105/239‑1 (MQ=255)
                                                                                                                                                                  gAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTcagc                                                                          <  1:836398/238‑1 (MQ=255)
                                                                                                                                                                       gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGaa                                                                     <  1:540891/238‑1 (MQ=255)
                                                                                                                                                                          ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAgg                                                                 <  1:363146/239‑1 (MQ=255)
                                                                                                                                                                           cTTACGGAACCTATCTTTTTGCAGTTGTTTCGCGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGa                                                                <  2:352910/239‑1 (MQ=255)
                                                                                                                                                                           cTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGa                                                                <  1:360570/239‑1 (MQ=255)
                                                                                                                                                                           cTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGa                                                                <  2:456483/239‑1 (MQ=255)
                                                                                                                                                                           cTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGa                                                                <  2:704960/239‑1 (MQ=255)
                                                                                                                                                                            ttACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATGCAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGAt                                                               <  1:1051155/239‑1 (MQ=255)
                                                                                                                                                                            ttACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAgatgat                                                               <  1:426844/239‑1 (MQ=255)
                                                                                                                                                                            ttACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGAt                                                               <  1:396475/239‑1 (MQ=255)
                                                                                                                                                                                      agacgttttataGTGGTTGCACGTTGTGGCGGGCATATAAAGCATGCACTTTCTGACACAACGGACGCAGGGTATAACCGGCGTAGCGCTCCTCATGCTGTTTGTAGATGGAAGGTAGCTCTTCTGATAATGGCGCGTGGGACTCAAGATGTTTTGCGGAGCGTACCAGCAGGGCAACGAGTTGCTGAAGTTTGGCCATATTTTTTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAGGa                                                     <  2:692207/228‑1 (MQ=255)
                                                                                                                                                                                                                        aTACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAAGATCGCATTTTTCCGGCACTACGCCATTTTCGCGCa                   <  2:301384/239‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCATTTCCTGACTCCACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGGTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGGTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAAGATCGCATTTTTCCGGCACTACGCCATTTTCGCGCAGGAAGTTAGc         <  2:1005805/239‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAAGATCGCATTTTTCCGGCACTACGCCATTTTCGCGCAGGAAGTTAGc         >  2:16591/1‑239 (MQ=255)
                                                                                                                                                                                                                                  tGCATTGTCTGATACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAAGATCGCATTTTTCCGGCACTACGCCATTTTCGCGCAGGAAGTTAGc         <  1:1094722/239‑1 (MQ=255)
                                                                                                                                                                                                                                   gCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAAGATCGCATTTTTCCGGCACTACGCCATTTTCGCGCAGGAAGTTAGCa        >  1:614072/1‑239 (MQ=255)
                                                                                                                                                                                                                                    cATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAAGATCGCATTTTTCCGGCACTACGCCATTTTCGCGCAGGAAGTTAGCaa       >  1:59798/1‑239 (MQ=255)
                                                                                                                                                                                                                                        tCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAAGATCGCATTTTTCCGGCACTACGCCATTTTCGCGCAGGAAGTTAGCAAGaa    >  1:1072759/1‑238 (MQ=255)
                                                                                                                                                                                                                                         ccTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAAGATCGCATTTTTCCGGCACTACGCCATTTTCGCGCAGGAAGTTAGCAAGAAtc  >  1:571559/1‑239 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               
ACCACCCCAGATTTCACCCGGAACAACACACAGCACACCCGGAGGATAAGGAAGCGCACCTTCGGCAGCGATACGGCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGAAGAGGATGGAGTTAAGATCGCATTTTTCCGGCACTACGCCATTTTCGCGCAGGAAGTTAGCAAGAATC  >  NC_000913/718436‑718907

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