Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,489,862 A→G 21.7% S107P (TCC→CCC)  argD ← bifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,489,8620AG21.7% 129.4 / 3.6 61S107P (TCC→CCC) argDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (13/32);  new base G (13/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCATCCATCACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCC‑TTGCCTTGCTGATCCCAGATTCGGCTGCCCT  >  NC_000913/3489625‑3490086
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      
gATCATCCATCACCGCTTTCCATGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAg                                                                                                                                                                                                                                       <  1:515716/239‑1 (MQ=255)
    atccatcACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCCTTTGCTT‑TC‑CGGCCGGAGtcat                                                                                                                                                                                                                                   >  1:278595/1‑236 (MQ=255)
      ccatcACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGGTTTGTTAGCTT‑CC‑GTGCCGGAGTT‑CCTa                                                                                                                                                                                                                                >  1:498884/1‑239 (MQ=255)
      ccatcACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGGTTCGTTAGCTT‑CC‑GTGCCGGGGGT‑CATa                                                                                                                                                                                                                                >  1:45986/1‑239 (MQ=255)
       catcACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGTCCACCCACCGAACCGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCATGCCAGTTTAAAGGCGGTTTCGTTAGGTT‑CC‑CGGCCGGaaggtagaa                                                                                                                                                                                                                                >  2:622752/1‑231 (MQ=255)
       catcACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGATT‑CC‑CGGCCGGAGTG‑CCTaa                                                                                                                                                                                                                               >  1:360073/1‑239 (MQ=255)
        tcccccccTTTCCCTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATaaa                                                                                                                                                                                                                              <  2:328033/232‑1 (MQ=255)
           aCCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTTGTTTTGAACGGGCTATGACGACCACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTGCGTAGATT‑CC‑ATCCCCGGGTTCCATAAac                                                                                                                                                                                                                             >  1:916056/1‑238 (MQ=255)
               cTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GGGCCGGGGTT‑CTTAAAca‑c‑cac                                                                                                                                                                                                                      >  1:132972/1‑238 (MQ=255)
                       gCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGACCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCAACCAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAACTT‑CCGGTGCCGGGGTC‑AATAAACC‑C‑CCcggttacgg                                                                                                                                                                                                              >  1:880530/1‑232 (MQ=255)
                        cATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTCAGTTT‑CC‑GTGCCCGAGTT‑CATAAa‑a‑c‑cacaccgttctgac                                                                                                                                                                                                           >  1:118516/1‑228 (MQ=255)
                             agaTCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTTGTTAGCTT‑TC‑CTGCCGGGGTT‑CATAAACA‑C‑CCCCCGTTCGGCAAAcct                                                                                                                                                                                                       >  1:638912/1‑237 (MQ=255)
                              gaTTGTTAAAGGTGACGTGGATGATGTGTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTGGGCTGCCCACCCCACGAAACGGTAAACCAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTGTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACg                                                                                                                                                                                                        <  1:251518/238‑1 (MQ=255)
                                  gTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTTGGTAACTT‑CC‑GTGCCGGGGTT‑CCTAAAAA‑C‑CCACCGTTCGGCAAACGTTgtc                                                                                                                                                                                                   >  1:283433/1‑236 (MQ=255)
                                  gTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTTGCTTCCC‑GTCCCGG‑GTT‑CCTAAAAA‑A‑CACCCGTTCGGCAAACGGTGcc                                                                                                                                                                                                   >  1:785043/1‑238 (MQ=255)
                                      aaGGGAACGGGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑At                                                                                                                                                                                             <  1:622752/239‑1 (MQ=255)
                                           aaCGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CCTAAAAA‑C‑CACCCGTTCGGCAAACCTTGCCTCA‑Atacgtt                                                                                                                                                                                        >  1:544857/1‑234 (MQ=255)
                                           aaCGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTGTCGCTTTCTT‑TC‑GGGCCCGGGTT‑TATAAACA‑CACAC‑CGTTTGGCAAACGTTTcccccaacaagt                                                                                                                                                                                         >  2:702750/1‑229 (MQ=255)
                                           aaCGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGGTTTGGTAGATT‑TC‑CGGCCGGGGTT‑CATAAAAC‑C‑CACTGGGTCGGCAAACGTTGCCTCC‑ATAAGtt                                                                                                                                                                                        >  1:616463/1‑239 (MQ=255)
                                           aaCGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGGTTCGTTAGCTT‑CC‑GTGCCGGAGGT‑CATAAACA‑C‑CACGCGGTCGGCAAACGTTGTCTCA‑ATCAggt                                                                                                                                                                                        >  1:543868/1‑237 (MQ=255)
                                           aaCGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGCTTCTTTTGCCT‑CC‑GGTCCCGGGGT‑CATAACAA‑C‑AACCCGTGTGGCCAACGGTGACTCa‑aacc                                                                                                                                                                                           >  2:126879/1‑233 (MQ=255)
                                             cGTTGATGATGTCTGCCGGTTGCGGCCCAAAGCCGTCGGAATATTGTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTa                                                                                                                                                                                      <  1:265070/239‑1 (MQ=255)
                                             cGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTG‑CCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCCAACAGGCGTAATGGCGTGCCAGTTTAAAAGCGGGTTTGCTTGCTTCCC‑GTCCCGGGGTT‑‑AAAAAAA‑C‑CACGCGGTCGTCCAACGTTTGCTCACACCAGttt                                                                                                                                                                                       >  1:882237/1‑238 (MQ=255)
                                              gTGGGTGATGTGTGCCGTTTTCGGCCCAAAGCCGGTGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTAc                                                                                                                                                                                     <  2:482279/239‑1 (MQ=255)
                                                     atgTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑TC‑GTGCCGGGGTT‑CATAAACA‑C‑AACGCGTTCGGCAAACGTTGCCTCC‑ATCCGTTTTCGccccaa                                                                                                                                                                              >  1:131724/1‑236 (MQ=255)
                                                            ccGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑TGGCCGGAGTT‑CATAAACA‑C‑CACGCGTTTGGCAAACGTTGCCCCC‑ACCAGTTTACGCCCCACACGCac                                                                                                                                                                        >  1:881255/1‑237 (MQ=255)
                                                                       ccAACGCCGTCGGATTATTTTGGCGTCCCACCCGCCGAAAGGGTAAA‑AAGCGAGCGAACATGAAAAGCGTTATGGAAGGCAATAAGTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCAtt                                                                                                                                                            <  2:715732/239‑1 (MQ=255)
                                                                       ccAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCAt                                                                                                                                                             <  1:48944/238‑1 (MQ=255)
                                                                          aaGCCGTGGGAATATTTTGGGTGCCCACCCCCCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGt                                                                                                                                                         <  2:881255/239‑1 (MQ=255)
                                                                            cccGGCGGAATTTTTTGGCCGCCCCCCCACCGAAACGGTAAA‑CAGGGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGa                                                                                                                                                       <  2:132972/238‑1 (MQ=255)
                                                                              cggaggAATATGTTGGCTGCCCCCCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGa                                                                                                                                                       <  1:126879/233‑1 (MQ=255)
                                                                              cGTCGGAATTTTTTGGCCGCCCACCCAACGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGaa                                                                                                                                                      <  2:283433/238‑1 (MQ=255)
                                                                              cGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGaaa                                                                                                                                                     <  1:143149/239‑1 (MQ=255)
                                                                               ggcGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGaaaa                                                                                                                                                    <  1:279403/237‑1 (MQ=255)
                                                                               gTCGGAATTTTTTGGCTGCCCACCCCCCGAAACGGTAAA‑AAGCGAGCGACCATGAAAAGGTTTAAGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGaa                                                                                                                                                      <  2:785043/237‑1 (MQ=255)
                                                                                       tttttGGCCGCCCCCCCACCGAAACGGTAAA‑CAGCGAGCGACCAGGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTg                                                                                                                                            <  2:544857/238‑1 (MQ=255)
                                                                                          ttGGCCGCCCCCCCACCGAAACGGTAAA‑CAGCGGGCGACCATGAAAGGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGata                                                                                                                                         <  2:278595/239‑1 (MQ=255)
                                                                                                 cccacccaCCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCcaca                                                                                                                                  <  1:823203/239‑1 (MQ=255)
                                                                                                  ccccccaacGAAACGGTAAG‑AAGCGATGGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCcaca                                                                                                                                  <  1:369613/235‑1 (MQ=255)
                                                                                                            aaaCGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGccc                                                                                                                       <  2:716294/239‑1 (MQ=255)
                                                                                                                ggTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGGTCGGCAAACGTTGCCTCA‑ATCAGGTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGGTATGCCCCAGAGTTTCGCCCC‑ggg                                                                                                                  >  1:223783/1‑239 (MQ=255)
                                                                                                                           gagcgaCCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGGTTACGCCCAAGACGCAGCGCCGGGTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCTGGGGTTTTTaagc                                                                                                         >  1:80305/1‑237 (MQ=255)
                                                                                                                             gcgaCCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAAcgcg                                                                                                       <  1:702750/238‑1 (MQ=255)
                                                                                                                                    tGAAAAGCGTTATGGAAGGCAATAATTCTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGaa                                                                                                                                                      <  1:503054/185‑1 (MQ=255)
                                                                                                                                           ggTTAGGGAAGGCAATAATTTTGGTTTTGAACGGGCTATTACGCCCACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGa                                                                                        <  1:284918/238‑1 (MQ=255)
                                                                                                                                             ttttGGAAGGCAATAATTTTGGGTTTGAACGGGTTTTGACGCACCCAGGCGTAAGGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGAt                                                                                       <  2:638912/235‑1 (MQ=255)
                                                                                                                                                         attatTTTGGTTTTGAACGGGTTTGGGCCCACAAGGGCGTAATGGCGTGCCAGTTTAAAGGGGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACC‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGCTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCCaa                                                                          <  2:360073/236‑1 (MQ=255)
                                                                                                                                                          tgatttgtgTTTTGAACGGTCTATGACGCACCCGGGCGTACGGGGCTGCCAGTTTAAGGTCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACC‑C‑CGCGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTCCCGCCCAAGACGCCTCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCAAAc                                                                         <  2:157509/231‑1 (MQ=255)
                                                                                                                                                                                     gCACCCAGGCGTATTGGGGTGACAGTTTAAAGGGTGTTTACTTAGCTT‑CC‑GGGCCGGGGTT‑CATAAACC‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGaaa                                              <  2:45986/239‑1 (MQ=255)
                                                                                                                                                                                      cacacaGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAAt                                             <  1:131272/239‑1 (MQ=255)
                                                                                                                                                                                       acacaGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATa                                                                                                                                                                                                                                <  2:442508/63‑1 (MQ=255)
                                                                                                                                                                                       acacaGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATa                                                                                                                                                                                                                                >  1:442508/1‑63 (MQ=255)
                                                                                                                                                                                       acacaGGCGTAATGGCGTGCCAGTTTAAAGGCGGGTTCGTTAGCTT‑CC‑GTGCCCGGGTT‑CATAAACA‑A‑CACCCGGTCCGCAAACCTTGCCTCA‑ATCAGTTTACGCCCCAGACGCAGCCACGGGTCATTGGGGAAAACGTTACTAAAATGCCACAAGGATTCGCCCT‑GGGTTTTTAAAGCGTTCACACCCCCCGGAACGGAAAGGCCCCACCCCGTACATCCAATGCCAccccccgaaac                                            >  2:518889/1‑231 (MQ=255)
                                                                                                                                                                                           aGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGGGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGGTTACGCCCAAGACCCCGCGCCGGGTCATTGGTGAAAACGTTACTGATATTCCACAGAGGTTCGCCCT‑GGGGTTTTAAAGCGTTCACCCACCCCGGGTGGCAAAGGCCCCACCCCCTAAATACAATGACACCCCCGCAAACGACa                                        >  2:118276/1‑239 (MQ=255)
                                                                                                                                                                                                tAATGGCGTGCCAGTTTAAAGGCGTTTTCCTTAGGTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGGGTTCGGCAAACGTTGCCTCA‑ATCAGGTTACGCCCAAGCCGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACt                                    <  1:652886/238‑1 (MQ=255)
                                                                                                                                                                                                       gTGCCAGTTTAAAGGGGGTTTCGTTTGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CCCGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTcc‑ttgcc                            <  2:498884/238‑1 (MQ=255)
                                                                                                                                                                                                         gCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGGTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATTCCACCCCCCGAATCGACAAAATcc‑ttgcctt                          >  1:715000/1‑238 (MQ=255)
                                                                                                                                                                                                                        gTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGGTTACGCCCAAGACACAGCCCCGGTTCATTGGGGGAAACGTTACTGTTATGCCACAGAGTTTCGCCCT‑GGGGTTTTAACTCGATCACCCGCGCCGAGTGGGAATGGGCCCACGGCGTAACTGCAAAGCCACCCCCGCAATCGACAAACCCCTTTCCCTTCCTGACCCCacgta            >  2:833067/1‑234 (MQ=255)
                                                                                                                                                                                                                                tgtct‑CC‑GTGCCGGGTTT‑CA‑AAACACC‑CAACCGGTCGGGAAAGGTGTCCTCA‑ATCAGTTTACGCGCAAGTAGCAGGTCCGTTTCAATGGTGAAAACGTTACTGATATGCCACAGAGGTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCC‑TTGCCTTGCTGATCCCAGATTCGGCTGCCCt  <  2:916056/235‑1 (MQ=255)
                                                                                                                                                                                                                                aGCTT‑CC‑GTGCCGGAGTT‑CATAAACC‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCC‑TTGCCTTGCTGATCCCAGATTCGGCTGCCCt  <  1:138094/239‑1 (MQ=255)
                                                                                                                                                                                                                                   tt‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCAc                                                                                                  <  2:16028/141‑1 (MQ=255)
                                                                                                                                                                                                                                   tt‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCAc                                                                                                  >  1:16028/1‑141 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      
GATCATCCATCACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGTTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCC‑TTGCCTTGCTGATCCCAGATTCGGCTGCCCT  >  NC_000913/3489625‑3490086

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