Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,489,860 A→G 21.9% 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,8600AG21.9% 137.1 / 4.5 68S107P (TCC→CCC) argDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (14/36);  new base G (15/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 2.55e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCATCCATCACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCC‑TTGCCTTGCTGATCCCAGATTCGGCTGCCCT  >  NC_000913/3489623‑3490084
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
gATCATCCATCACCGCTTTCCATGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAg                                                                                                                                                                                                                                        <  1:515716/239‑1 (MQ=255)
    atccatcACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCCTTTGCTT‑TC‑CGGCCGGAGtcat                                                                                                                                                                                                                                    >  1:278595/1‑236 (MQ=255)
      ccatcACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑GTTTGTTAGCTT‑CC‑GTGCCGGAGTT‑CCTa                                                                                                                                                                                                                                 >  1:498884/1‑239 (MQ=255)
      ccatcACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑GTTCGTTAGCTT‑CC‑GTGCCGGGGGT‑CATa                                                                                                                                                                                                                                 >  1:45986/1‑239 (MQ=255)
       catcACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGTCCACCCACCGAACCGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCATGCCAGTTTAAAGGCGG‑TTTCGTTAGGTT‑CC‑CGGCCGGaaggtagaa                                                                                                                                                                                                                                 >  2:622752/1‑231 (MQ=255)
       catcACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGATT‑CC‑CGGCCGGAGTG‑CCTaa                                                                                                                                                                                                                                >  1:360073/1‑239 (MQ=255)
        tcccccccTTTCCCTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATaaa                                                                                                                                                                                                                               <  2:328033/232‑1 (MQ=255)
           aCCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTTGTTTTGAACGGGCTATGACGACCACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTGCGTAGATT‑CC‑ATCCCCGGGTTCCATAAac                                                                                                                                                                                                                              >  1:916056/1‑238 (MQ=255)
            cctctttcctgcgaTGGAGATCGTTAAAGGGAACGGGGATGAGGCCTGCCGGTTCCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAAca‑c‑c                                                                                                                                                                                                                         <  2:632608/226‑1 (MQ=255)
               cTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GGGCCGGGGTT‑CTTAAAca‑c‑cac                                                                                                                                                                                                                       >  1:132972/1‑238 (MQ=255)
                       gCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGACCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCAACCAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAACTT‑CCGGTGCCGGGGTC‑AATAAACC‑C‑CCcggttacgg                                                                                                                                                                                                               >  1:880530/1‑232 (MQ=255)
                        cATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTCAGTTT‑CC‑GTGCCCGAGTT‑CATAAa‑a‑c‑cacaccgttctgac                                                                                                                                                                                                            >  1:118516/1‑228 (MQ=255)
                             agaTCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTTGTTAGCTT‑TC‑CTGCCGGGGTT‑CATAAACA‑C‑CCCCCGTTCGGCAAAcct                                                                                                                                                                                                        >  1:638912/1‑237 (MQ=255)
                              gaTTGTTAAAGGTGACGTGGATGATGTGTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTGGGCTGCCCACCCCACGAAACGGTAAACCAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTGTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACg                                                                                                                                                                                                         <  1:251518/238‑1 (MQ=255)
                                  gTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGGCTTTCG‑TAGCTT‑CC‑CTGCCGGGGGT‑CACAAAca‑c‑caacccgtcgggaaacagtggc                                                                                                                                                                                                    >  2:279403/1‑218 (MQ=255)
                                  gTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTTGGTAACTT‑CC‑GTGCCGGGGTT‑CCTAAAAA‑C‑CCACCGTTCGGCAAACGTTgtc                                                                                                                                                                                                    >  1:283433/1‑236 (MQ=255)
                                  gTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTTGCTTCCC‑GTCCCGG‑GTT‑CCTAAAAA‑A‑CACCCGTTCGGCAAACGGTGcc                                                                                                                                                                                                    >  1:785043/1‑238 (MQ=255)
                                      aaGGGAACGGGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑At                                                                                                                                                                                              <  1:622752/239‑1 (MQ=255)
                                           aaCGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CCTAAAAA‑C‑CACCCGTTCGGCAAACCTTGCCTCA‑Atacgtt                                                                                                                                                                                         >  1:544857/1‑234 (MQ=255)
                                           aaCGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TGTCGCTTTCTT‑TC‑GGGCCCGGGTT‑TATAAACA‑CACAC‑CGTTTGGCAAACGTTTcccccaacaagt                                                                                                                                                                                          >  2:702750/1‑229 (MQ=255)
                                           aaCGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑GTTTGGTAGATT‑TC‑CGGCCGGGGTT‑CATAAAAC‑C‑CACTGGGTCGGCAAACGTTGCCTCC‑ATAAGtt                                                                                                                                                                                         >  1:616463/1‑239 (MQ=255)
                                           aaCGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑GTTCGTTAGCTT‑CC‑GTGCCGGAGGT‑CATAAACA‑C‑CACGCGGTCGGCAAACGTTGTCTCA‑ATCAggt                                                                                                                                                                                         >  1:543868/1‑237 (MQ=255)
                                           aaCGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑CTTCTTTTGCCT‑CC‑GGTCCCGGGGT‑CATAACAA‑C‑AACCCGTGTGGCCAACGGTGACTCa‑aacc                                                                                                                                                                                            >  2:126879/1‑233 (MQ=255)
                                             cGTTGATGATGTCTGCCGGTTGCGGCCCAAAGCCGTCGGAATATTGTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTa                                                                                                                                                                                       <  1:265070/239‑1 (MQ=255)
                                             cGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTG‑CCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCCAACAGGCGTAATGGCGTGCCAGTTTAAAAGCGG‑GTTTGCTTGCTTCCC‑GTCCCGGGGTT‑‑AAAAAAA‑C‑CACGCGGTCGTCCAACGTTTGCTCACACCAGttt                                                                                                                                                                                        >  1:882237/1‑238 (MQ=255)
                                              gTGGGTGATGTGTGCCGTTTTCGGCCCAAAGCCGGTGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTAc                                                                                                                                                                                      <  2:482279/239‑1 (MQ=255)
                                                     atgTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑TC‑GTGCCGGGGTT‑CATAAACA‑C‑AACGCGTTCGGCAAACGTTGCCTCC‑ATCCGTTTTCGccccaa                                                                                                                                                                               >  1:131724/1‑236 (MQ=255)
                                                            ccGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑TGGCCGGAGTT‑CATAAACA‑C‑CACGCGTTTGGCAAACGTTGCCCCC‑ACCAGTTTACGCCCCACACGCac                                                                                                                                                                         >  1:881255/1‑237 (MQ=255)
                                                                       ccAACGCCGTCGGATTATTTTGGCGTCCCACCCGCCGAAAGGGTAAA‑AAGCGAGCGAACATGAAAAGCGTTATGGAAGGCAATAAGTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCAtt                                                                                                                                                             <  2:715732/239‑1 (MQ=255)
                                                                       ccAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCAt                                                                                                                                                              <  1:48944/238‑1 (MQ=255)
                                                                          aaGCCGTGGGAATATTTTGGGTGCCCACCCCCCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGt                                                                                                                                                          <  2:881255/239‑1 (MQ=255)
                                                                            cccGGCGGAATTTTTTGGCCGCCCCCCCACCGAAACGGTAAA‑CAGGGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGa                                                                                                                                                        <  2:132972/238‑1 (MQ=255)
                                                                              cggaggAATATGTTGGCTGCCCCCCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGa                                                                                                                                                        <  1:126879/233‑1 (MQ=255)
                                                                              cGTCGGAATTTTTTGGCCGCCCACCCAACGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGaa                                                                                                                                                       <  2:283433/238‑1 (MQ=255)
                                                                              cGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGaaa                                                                                                                                                      <  1:143149/239‑1 (MQ=255)
                                                                               ggcGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGaaaa                                                                                                                                                     <  1:279403/237‑1 (MQ=255)
                                                                               gTCGGAATTTTTTGGCTGCCCACCCCCCGAAACGGTAAA‑AAGCGAGCGACCATGAAAAGGTTTAAGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGaa                                                                                                                                                       <  2:785043/237‑1 (MQ=255)
                                                                                       tttttGGCCGCCCCCCCACCGAAACGGTAAA‑CAGCGAGCGACCAGGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTg                                                                                                                                             <  2:544857/238‑1 (MQ=255)
                                                                                          ttGGCCGCCCCCCCACCGAAACGGTAAA‑CAGCGGGCGACCATGAAAGGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGata                                                                                                                                          <  2:278595/239‑1 (MQ=255)
                                                                                                 cccacccaCCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCcaca                                                                                                                                   <  1:823203/239‑1 (MQ=255)
                                                                                                  ccgcccgcggagaCGGGAAA‑CAGCGGGGGCCCATGAAAAGGGTTAGGGAAGGCAATAAGTTTGGTTTGGAAAGGGATATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCTCCGGTTCATTGGTGAAAACGTTACTGATATGCaaca                                                                                                                                   <  2:326347/226‑1 (MQ=255)
                                                                                                  ccccccaacGAAACGGTAAG‑AAGCGATGGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCcaca                                                                                                                                   <  1:369613/235‑1 (MQ=255)
                                                                                                            aaaCGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGccc                                                                                                                        <  2:716294/239‑1 (MQ=255)
                                                                                                                ggTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTTGCTT‑CC‑GTGCCGCAGTT‑CGTAAAAC‑C‑CCCCCGCTCGGCAAACGTTGCCTCA‑CTCAGTTTAGCCCCAAGACCCACAGCCGGTTCCATTGTGCAAAAGTCACTGATAAGcacccacaactacccccggg                                                                                                                    >  2:408246/1‑220 (MQ=255)
                                                                                                                ggTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGGTCGGCAAACGTTGCCTCA‑ATCAGGTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGGTATGCCCCAGAGTTTCGCCCC‑ggg                                                                                                                   >  1:223783/1‑239 (MQ=255)
                                                                                                                           gagcgaCCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGGTTACGCCCAAGACGCAGCGCCGGGTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCTGGGGTTTTTaagc                                                                                                          >  1:80305/1‑237 (MQ=255)
                                                                                                                             gcgaCCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAAcgcg                                                                                                        <  1:702750/238‑1 (MQ=255)
                                                                                                                                    tGAAAAGCGTTATGGAAGGCAATAATTCTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGaa                                                                                                                                                       <  1:503054/185‑1 (MQ=255)
                                                                                                                                           ggTTAGGGAAGGCAATAATTTTGGTTTTGAACGGGCTATTACGCCCACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGa                                                                                         <  1:284918/238‑1 (MQ=255)
                                                                                                                                             ttttGGAAGGCAATAATTTTGGGTTTGAACGGGTTTTGACGCACCCAGGCGTAAGGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑t                                                                                       <  2:638912/235‑1 (MQ=255)
                                                                                                                                                         attatTTTGGTTTTGAACGGGTTTGGGCCCACAAGGGCGTAATGGCGTGCCAGTTTAAAGGGGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACC‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGCTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCaa                                                                          <  2:360073/236‑1 (MQ=255)
                                                                                                                                                          tgatttgtgTTTTGAACGGTCTATGACGCACCCGGGCGTACGGGGCTGCCAGTTTAAGGTCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACC‑C‑CGCGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTCCCGCCCAAGACGCCTCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCAAAc                                                                         <  2:157509/231‑1 (MQ=255)
                                                                                                                                                                                     gCACCCAGGCGTATTGGGGTGACAGTTTAAAGGGTG‑TTTACTTAGCTT‑CC‑GGGCCGGGGTT‑CATAAACC‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGaaa                                              <  2:45986/239‑1 (MQ=255)
                                                                                                                                                                                      ccaccgggggttaTGGCGTTCCAGTTTAAAGGCGG‑TTTCTTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACC‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAAt                                             <  2:118516/227‑1 (MQ=255)
                                                                                                                                                                                      cacacaGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAAt                                             <  1:131272/239‑1 (MQ=255)
                                                                                                                                                                                       acacaGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATa                                                                                                                                                                                                                                 <  2:442508/63‑1 (MQ=255)
                                                                                                                                                                                       acacaGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATa                                                                                                                                                                                                                                 >  1:442508/1‑63 (MQ=255)
                                                                                                                                                                                       acacaGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑GTTCGTTAGCTT‑CC‑GTGCCCGGGTT‑CATAAACA‑A‑CACCCGGTCCGCAAACCTTGCCTCA‑ATCAGTTTACGCCCCAGACGCAGCCACGGGTCATTGGGGAAAACGTTACTAAAATGCCACAAGGATTCGCCCT‑GGGTTTTTAAAGCGTTCACACCCCCCGGA‑ACGGAAAGGCCCCACCCCGTACATCCAATGCCAccccccgaaac                                            >  2:518889/1‑231 (MQ=255)
                                                                                                                                                                                           aGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGGGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGGTTACGCCCAAGACCCCGCGCCGGGTCATTGGTGAAAACGTTACTGATATTCCACAGAGGTTCGCCCT‑GGGGTTTTAAAGCGTTCACCCACCCCGGG‑TGGCAAAGGCCCCACCCCCTAAATACAATGACACCCCCGCAAACGACa                                        >  2:118276/1‑239 (MQ=255)
                                                                                                                                                                                             gCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGGGTT‑CATAAACA‑C‑CCCGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCACCGCCGGTTCATTGGGGAAAACGTTACTGATATGCCCCCGAGTTTTGCCCT‑GGGGTTTTAAACGGTTCACCCGCGCCGGAGTGGCAAT‑GCCCCACACCGTAAATGCAAAGcccccccccgaaactcaaaa                                      >  2:379547/1‑225 (MQ=255)
                                                                                                                                                                                                tAATGGCGTGCCAGTTTAAAGGCGT‑TTTCCTTAGGTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGGGTTCGGCAAACGTTGCCTCA‑ATCAGGTTACGCCCAAGCCGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACt                                    <  1:652886/238‑1 (MQ=255)
                                                                                                                                                                                                       gTGCCAGTTTAAAGGGGG‑TTTCGTTTGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CCCGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTcc‑ttgcc                            <  2:498884/238‑1 (MQ=255)
                                                                                                                                                                                                         gCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGGTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATTCCACCCCCCGAATCGACAAAATcc‑ttgcctt                          >  1:715000/1‑238 (MQ=255)
                                                                                                                                                                                                                        g‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGGTTACGCCCAAGACACAGCCCCGGTTCATTGGGGGAAACGTTACTGTTATGCCACAGAGTTTCGCCCT‑GGGGTTTTAACTCGATCACCCGCGCCGAG‑TGGGAATGGGCCCACGGCGTAACTGCAAAGCCACCCCCGCAATCGACAAACCCCTTTCCCTTCCTGACCCCacgta            >  2:833067/1‑234 (MQ=255)
                                                                                                                                                                                                                             tttctttctttcctgtccGGAGTT‑CAAAAACA‑C‑CACCCGTGCGGCAAACGTGGCCCCA‑ATCAGTTTACGCCCCAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCC‑TTGCCTTGCTGATCCCAGATTCgg         <  2:543868/222‑1 (MQ=255)
                                                                                                                                                                                                                                 tgtct‑CC‑GTGCCGGGTTT‑CA‑AAACACC‑CAACCGGTCGGGAAAGGTGTCCTCA‑ATCAGTTTACGCGCAAGTAGCAGGTCCGTTTCAATGGTGAAAACGTTACTGATATGCCACAGAGGTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCC‑TTGCCTTGCTGATCCCAGATTCGGCTGCCCt  <  2:916056/235‑1 (MQ=255)
                                                                                                                                                                                                                                 aGCTT‑CC‑GTGCCGGAGTT‑CATAAACC‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCC‑TTGCCTTGCTGATCCCAGATTCGGCTGCCCt  <  1:138094/239‑1 (MQ=255)
                                                                                                                                                                                                                                    tt‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCAc                                                                                                   >  1:16028/1‑141 (MQ=255)
                                                                                                                                                                                                                                    tt‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCAc                                                                                                   <  2:16028/141‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       
GATCATCCATCACCGCTTTCACTGCATGGAGATCGTTAAAGGGAACGTGGATGATGTCTGCCGGTTTCGGCCCAAAGCCGTCGGAATATTTTGGCTGCCCACCCACCGAAACGGTAAA‑CAGCGAGCGACCATGAAAAGCGTTATGGAAGGCAATAATTTTGGTTTTGAACGGGCTATGACGCACACAGGCGTAATGGCGTGCCAGTTTAAAGGCGG‑TTTCGTTAGCTT‑CC‑GTGCCGGAGTT‑CATAAACA‑C‑CACGCGTTCGGCAAACGTTGCCTCA‑ATCAGTTTACGCCCAAGACGCAGCGCCGGTTCATTGGTGAAAACGTTACTGATATGCCACAGAGTTTCGCCCT‑GGGTTTTTAACGCGTTCACCAGCGCAGGA‑TGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCC‑TTGCCTTGCTGATCCCAGATTCGGCTGCCCT  >  NC_000913/3489623‑3490084

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: ‑