Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,239,855 A→G 17.6% E326G (GAG→GGG)  dadX → alanine racemase, catabolic, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,239,8550AG17.6% 112.4 / 2.1 57E326G (GAG→GGG) dadXalanine racemase, catabolic, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base A (15/31);  new base G (7/4);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 8.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGCGATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGG  >  NC_000913/1239653‑1240091
                                                                                                                                                                                                           |                                                                                                                                                                                                                                                
gcgcgcGATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTCAGCCGGGTGGGCCGTGGGGCAAGGAGGTCACAAATGGATG‑TGTcggcgcccgcgccggcaa                                                                                                                                                                                                             >  1:894877/1‑228 (MQ=255)
 cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGTGGGGGAAGGAGGTCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACg                                                                                                                                                                                                            >  1:365544/1‑239 (MQ=255)
 cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCCGTTGGGGCAGGGGATCA‑AAAATGTTGATTGTGCCGCCCCGGCCGGaaag                                                                                                                                                                                                            >  1:776991/1‑237 (MQ=255)
 cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTTGGCTGTGGGGCAAGGAGGTCA‑AAAATGATGATGTCTCCGCCGCTGCCGGGACg                                                                                                                                                                                                            >  1:235213/1‑239 (MQ=255)
 cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTAAGCCGGTTTGGCTGTGGGGGAAGGGGGTCA‑CAATTTTTGATTTGCCCGCCCCTGCCGGGACg                                                                                                                                                                                                            >  1:1228248/1‑239 (MQ=255)
 cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGAAGCCGGTTGAGCCTGGGGGCACGGGGGTCA‑AAATTGATGATGTGGCCGCCGCGGCCGGaaac                                                                                                                                                                                                            >  1:158986/1‑237 (MQ=255)
   aggagacagcAGGGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGt                                                                                                                                                                                                          <  1:1041484/230‑1 (MQ=255)
    gcgATTAACCGGCAATCGGCATTGGCGCCGCAGGGTACGCCCGCGGTTATCCGCGCCACGCG‑GCTGCCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGt                                                                                                                                                                                                          <  1:1267717/238‑1 (MQ=255)
      gATGAACCGCGAATGCGCATTTTCCCCCCGGGGTACGCCGACGGTTACCCGGCCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTggg                                                                                                                                                                                                       <  1:1301635/239‑1 (MQ=255)
         gAACAGCGAATCGGCATGGCCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCt                                                                                                                                                                                                     <  1:967194/238‑1 (MQ=255)
          aaCAGGGAATCGGCATTGGCGCCGGAGGGTACGCAGACGGTTAGCCGCGCCACGCG‑CCTACCGGTACCCCGGTTTTAGTGGACGGCGTGCGCCCCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCt                                                                                                                                                                                                     <  1:844603/237‑1 (MQ=255)
                 attCGGCATTGTCCCCGCTGGGTCCGGCGCGGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGGTGGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTg                                                                                                                                                                                             <  1:541813/236‑1 (MQ=255)
                         tggTCGCCGCAGGGTTCGCCGGCGGTTACCCGCGCCACGGG‑CCGACCGGTACCCGTGTTTTAGGGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTgcgcg                                                                                                                                                                                   <  1:536608/237‑1 (MQ=255)
                          tggcgccgcAGGGTACGGCGACGGTTATCCGGGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTgcgcg                                                                                                                                                                                   <  1:1173024/235‑1 (MQ=255)
                                         gcgcGACGGTTACCGGCCGCAGGCG‑CGTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTc                                                                                                                                                                    <  1:298652/235‑1 (MQ=255)
                                            ccggggttgcccgcgcc‑cgcgcCCTACGGGTCCCCCTGGTTTAGGGGGCGGCGGGGCGACCATGACGGTGGGGACCGTCTCGGTGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTTTTGGTACGCCGGTTGAGCTGGGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCgg                                                                                                                                                                <  2:22239/229‑1 (MQ=255)
                                                   tacccGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGGGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGa                                                                                                                                                          <  1:842908/235‑1 (MQ=255)
                                                     tCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGAGGTGCCCGCCGCTGCCGGGACGGTGGGCGTTGAGGTTA‑TGTGGGCGCTGGGGCTACGCGTCCCGGGTGTGGCgg                                                                                                                                                       >  1:371455/1‑239 (MQ=255)
                                                                                cTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGTGGGGCAAGGGGATCA‑AAATTGATGAAGTCGCCGCCGCTGCCGGGACGGTGGGGTATGAGTTGA‑TGTGGGCGCTGGCGCTACCCGTCCCGGTTTGGACGGGGTTAC‑TTGTTGTTAGGCCGGTCggg                                                                                                                            >  1:547168/1‑238 (MQ=255)
                                                                                     ttggTGGGCGGCGTGGGCCCCATGGCGGTGGGGGCCGGCTCGATGGATACGCTAGCGGTCGATTTACGGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGGGCGGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCaa                                                                                                                       <  2:1069939/236‑1 (MQ=255)
                                                                                             cGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGGACGGGGGGCTATGAGTTGA‑TGTGCGGGCTGGCGCTACGCGTCCCGGTTGGGACGGGGTAAC‑TTGTTTTTAGACGGGTCGGCGGCAACGgct                                                                                                                  >  1:824201/1‑236 (MQ=255)
                                                                                                          aTGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGATACGCCGGGTGAGCTGTGGGGCCAGGAGATCC‑AAATTGGTGGTTGCGCCGCCGCTGCCGGGACGGGGGGCTATGAGGTTA‑TGTTCGCGCTGGGGCTACGCGTCCCCGTTGTGACGGGGTAAC‑TTTTTTTAAACCGGATCGGAGGCAACACCATCTGGGGGGaaaaaaa                                                                                                  >  2:533299/1‑239 (MQ=255)
                                                                                                                 tGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGGTGGGCAAGGGGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGGACGGGGGGGTATGTGGTGT‑TTTGCGGGCTGGCGCTACGCGTCCCGGGTGTGGCGGGGGAAA‑TTGTTTTTATGCGGGACGGGGGGCACGGCTTCTGGGTGCaaaaaaaaaaaccc                                                                                           >  2:23635/1‑232 (MQ=255)
                                                                                                                 tGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCCGGTGAGCTGGGGGGCAAGGAGAGCA‑AAATTGTTGTTGTCGCCGCCGCTGCCGGGAAGGAGGGCTAGGAGTTTATTGTCCG‑GCTGGGCCTAGCCGGCCCGGGTGGGACGGGGGGAC‑TTTTTGTTAGCCGGGCTCGAGGCAACCGCTTTTGGGGGGGAAAAAAAaatcc                                                                                            >  2:957579/1‑238 (MQ=255)
                                                                                                                  ggggCCCGTCGCGATGGATATGCTAGGGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATcatcca                                                                                           <  1:796644/238‑1 (MQ=255)
                                                                                                                   gggACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACCCCGGGTGAGCTGTGGGGCAAAGGGATCC‑AAATTGTTGATGTCGCCGCCCCTGCCGGGAGGGGGGGCTACGAGTTTA‑TTGGCGCGCTGGCGTGACGCGGCCCGGTGGGTGCGGGGTTACTTTGTTGTACGCGGGGTGGGGGTCACCGCTCTCTGGGGACaaaaaaacaaaccac                                                                                          >  2:983549/1‑231 (MQ=255)
                                                                                                                   gggACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTAAGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGGTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGGCCCGGGTGTGACGGGGTAAA‑TTTTTGTAAACCGGATCTGAGGCAACGTCTTCTTGGGGCAAAAAAAAcatccat                                                                                          >  1:500059/1‑238 (MQ=255)
                                                                                                                   gggACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGTTGAGGTGTGGGGCCAGGAGATCA‑AAATTGGTGTTGTCGCCGCCGCCGCCGGAACGGGGGGGTATGGGGTGG‑TGGTCACGCTGGCGCCCCCGGCCCCGGTTGGTGCGGGGTAAC‑TTTTTTGAAACCGGGGCGGGGGCACGGCGTTTTGGGGGCaaaaaaaaaataccc                                                                                          >  2:1181027/1‑230 (MQ=255)
                                                                                                                       ccGTCTCGATGGATTTGCTGGCGGTCGGTTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCgg                                                                                      <  1:15699/238‑1 (MQ=255)
                                                                                                                            tcGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTc                                                                                <  1:922797/239‑1 (MQ=255)
                                                                                                                                    tctGCGAGCGGTCGATTTAACGCCTTGCCCGCAGGCGTGTATTGGGACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACt                                                                        <  1:772204/236‑1 (MQ=255)
                                                                                                                                       gCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTGAGCTGTGGGGCAAGGGGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGGGGGCTGTGAGTTGA‑TGTGCGCGCTGGGGCTACGCGTCCCGGGTGGGCAGGGGTGAC‑TTTTTTTAAACGGGAGCGGGGGGCACGGCTCTTGGGGGCAAAAAAAACATACATACCGCTGGGGAGCA‑ACGGtgg                                                                    >  2:905932/1‑237 (MQ=255)
                                                                                                                                           gCGGTCGATTTAACGCCTTGCCCGCAGGGGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAgttgt                                                                <  1:319857/239‑1 (MQ=255)
                                                                                                                                                     tAACGGCTTGCCCGGAGGCGGGTATTGGTACGCCGGTTGAGCTGGGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACa                                                      <  1:14272/239‑1 (MQ=255)
                                                                                                                                                           ctgggcgggaGGCGGGTATTGGGACGCCGGTTGAGCTGGGGGGCACGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGGTGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGGGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  1:1241375/230‑1 (MQ=255)
                                                                                                                                                           cTTGCCCGGGGGGGGGTGGTGGTGCGGCGGTTGGGCTGGGGGGCAAGGAGGTCA‑AAATTGATGATGGCGCCGCCGCGGCCGGGACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  2:894877/239‑1 (MQ=255)
                                                                                                                                                           cTTGCCCGCGGGGGGTTTTTGGTACGGGGTTTGGGCTGGGGGGAAGGGGGATAA‑AAATGGATGATGTCGCCGCCGCTGCCGGAACGGTGGGGTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGGGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  2:371455/239‑1 (MQ=255)
                                                                                                                                                           cTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  1:905932/239‑1 (MQ=255)
                                                                                                                                                           cTTGCCCCGAGGGGGGTTTTGGTACGCGGGTTGAGCTGGGGGGCAAGGGGATAA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  2:824201/239‑1 (MQ=255)
                                                                                                                                                            ttGCCGGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGGCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCa                                                <  1:1205255/238‑1 (MQ=255)
                                                                                                                                                             tGCCCGCGGGCGGGTTTTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCAtt                                              <  1:23635/239‑1 (MQ=255)
                                                                                                                                                                ccGCGGGGGGGTTTTGGGTCGCCGGTTGGGCGGGGGGGCAGGGGGATCA‑AAATTGATGGTGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGc                                            <  2:776991/238‑1 (MQ=255)
                                                                                                                                                                caggaGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGcc                                           <  1:324321/235‑1 (MQ=255)
                                                                                                                                                                  gCAGGGGGGTATTGGTACGCCGGTTGAGGTGTGGGGCAAGTAGATCA‑AAATTGTTGCTGGCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCa                                         <  1:957579/239‑1 (MQ=255)
                                                                                                                                                                                   cGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAATGTCCATTaaa                        >  2:223807/1‑239 (MQ=255)
                                                                                                                                                                                    gCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTTTAAGCCGGATCGGAGGCACCGCCTTTGGGGGGCAAAAAAATCTTCAACCCGG‑TGGGCCGTACACGGAGGTTTTTAAGGTGAAAGAGCAATTGCCCATGGTATAGTCCATAaaaa                       >  2:1371628/1‑239 (MQ=255)
                                                                                                                                                                                    gCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAAAGCCATTGCCCATGATAGTGTACATTaaaa                       >  1:529056/1‑239 (MQ=255)
                                                                                                                                                                                    gCCGGTTGAGCTGTGGGGCAAGGAGATAA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTaaa                        <  1:1244113/238‑1 (MQ=255)
                                                                                                                                                                                     ccGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATAATATGGTCCATTAAAAg                      >  2:1359484/1‑239 (MQ=255)
                                                                                                                                                                                     ccGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATAATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGCAAGGGCCATTGCCCCTGATAGTGGCCATTAAAAg                      >  2:1244113/1‑239 (MQ=255)
                                                                                                                                                                                         tttggtgttgGGGAAAGGAGTACA‑AAATTGGTGGGGTCGCCGCCGCTTCCGGAACGGTGGGTTATGAGTTGA‑TGTGCGCGCTGGCGCTACGGGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATTCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATg                  <  2:547168/231‑1 (MQ=255)
                                                                                                                                                                                           ggggtttgGGGCAAGGGATTTA‑AATTTGATGGGTTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAA‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGa                <  2:235213/233‑1 (MQ=255)
                                                                                                                                                                                           gAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGa                <  1:223807/239‑1 (MQ=255)
                                                                                                                                                                                             gCTGTGGGGCAAGGAGATAA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGGCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCAGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGac               <  1:310620/238‑1 (MQ=255)
                                                                                                                                                                                              cTGTGGGGCAAGGAGATCA‑AAATTGAGGATGGCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGaca              <  1:553599/238‑1 (MQ=255)
                                                                                                                                                                                                gtcgggCAACGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGGGCAAAAAAATCATCCATCCGGCTGGCGAGCA‑ACTGTAGTTGGTAATTGGACAGAACCATTGCCCATGACAAGGGCCATTAAAAGGAAGGgcacaa           >  2:654893/4‑237 (MQ=255)
                                                                                                                                                                                                        aGGAAAACA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTcccc    <  1:654893/238‑1 (MQ=255)
                                                                                                                                                                                                           agaTCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCCGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCgg  <  1:429192/237‑1 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                                
GCGCGCGATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCG‑CCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCA‑AAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGA‑TGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAC‑TTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCA‑ACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGG  >  NC_000913/1239653‑1240091

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