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

GCGCGCGATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGG  >  NC_000913/1239653‑1240091
                                                                                                                                                                                                           |                                                                                                                                                                                                                                                
gcgcgcGATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTCAGCCGGGTGGGCCGTGGGGCAAGGAGGTCACAAATGGATGTGTcggcgcccgcgccggcaa                                                                                                                                                                                                             >  1:894877/1‑228 (MQ=255)
 cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGTGGGGGAAGGAGGTCAAAATTGATGATGTCGCCGCCGCTGCCGGAACg                                                                                                                                                                                                            >  1:365544/1‑239 (MQ=255)
 cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCCGTTGGGGCAGGGGATCAAAAATGTTGATTGTGCCGCCCCGGCCGGaaag                                                                                                                                                                                                            >  1:776991/1‑237 (MQ=255)
 cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTTGGCTGTGGGGCAAGGAGGTCAAAAATGATGATGTCTCCGCCGCTGCCGGGACg                                                                                                                                                                                                            >  1:235213/1‑239 (MQ=255)
 cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTAAGCCGGTTTGGCTGTGGGGGAAGGGGGTCACAATTTTTGATTTGCCCGCCCCTGCCGGGACg                                                                                                                                                                                                            >  1:1228248/1‑239 (MQ=255)
 cgcgcgATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGAAGCCGGTTGAGCCTGGGGGCACGGGGGTCAAAATTGATGATGTGGCCGCCGCGGCCGGaaac                                                                                                                                                                                                            >  1:158986/1‑237 (MQ=255)
   aggagacagcAGGGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGt                                                                                                                                                                                                          <  1:1041484/230‑1 (MQ=255)
    gcgATTAACCGGCAATCGGCATTGGCGCCGCAGGGTACGCCCGCGGTTATCCGCGCCACGCGGCTGCCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGt                                                                                                                                                                                                          <  1:1267717/238‑1 (MQ=255)
      gATGAACCGCGAATGCGCATTTTCCCCCCGGGGTACGCCGACGGTTACCCGGCCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTggg                                                                                                                                                                                                       <  1:1301635/239‑1 (MQ=255)
         gAACAGCGAATCGGCATGGCCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCt                                                                                                                                                                                                     <  1:967194/238‑1 (MQ=255)
          aaCAGGGAATCGGCATTGGCGCCGGAGGGTACGCAGACGGTTAGCCGCGCCACGCGCCTACCGGTACCCCGGTTTTAGTGGACGGCGTGCGCCCCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCt                                                                                                                                                                                                     <  1:844603/237‑1 (MQ=255)
                 attCGGCATTGTCCCCGCTGGGTCCGGCGCGGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGGTGGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTg                                                                                                                                                                                             <  1:541813/236‑1 (MQ=255)
                         tggTCGCCGCAGGGTTCGCCGGCGGTTACCCGCGCCACGGGCCGACCGGTACCCGTGTTTTAGGGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTgcgcg                                                                                                                                                                                   <  1:536608/237‑1 (MQ=255)
                          tggcgccgcAGGGTACGGCGACGGTTATCCGGGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTgcgcg                                                                                                                                                                                   <  1:1173024/235‑1 (MQ=255)
                                         gcgcGACGGTTACCGGCCGCAGGCGCGTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTc                                                                                                                                                                    <  1:298652/235‑1 (MQ=255)
                                            ccggggttgcccgcgcccgcgcCCTACGGGTCCCCCTGGTTTAGGGGGCGGCGGGGCGACCATGACGGTGGGGACCGTCTCGGTGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTTTTGGTACGCCGGTTGAGCTGGGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCgg                                                                                                                                                                <  2:22239/229‑1 (MQ=255)
                                                   tacccGCGCCACGCGCCTACCGGTACCCCTGTTTTAGGGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGa                                                                                                                                                          <  1:842908/235‑1 (MQ=255)
                                                     tCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGAGGTGCCCGCCGCTGCCGGGACGGTGGGCGTTGAGGTTATGTGGGCGCTGGGGCTACGCGTCCCGGGTGTGGCgg                                                                                                                                                       >  1:371455/1‑239 (MQ=255)
                                                                                cTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGTGGGGCAAGGGGATCAAAATTGATGAAGTCGCCGCCGCTGCCGGGACGGTGGGGTATGAGTTGATGTGGGCGCTGGCGCTACCCGTCCCGGTTTGGACGGGGTTACTTGTTGTTAGGCCGGTCggg                                                                                                                            >  1:547168/1‑238 (MQ=255)
                                                                                     ttggTGGGCGGCGTGGGCCCCATGGCGGTGGGGGCCGGCTCGATGGATACGCTAGCGGTCGATTTACGGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGGGCGGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCaa                                                                                                                       <  2:1069939/236‑1 (MQ=255)
                                                                                             cGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGGACGGGGGGCTATGAGTTGATGTGCGGGCTGGCGCTACGCGTCCCGGTTGGGACGGGGTAACTTGTTTTTAGACGGGTCGGCGGCAACGgct                                                                                                                  >  1:824201/1‑236 (MQ=255)
                                                                                                          aTGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGATACGCCGGGTGAGCTGTGGGGCCAGGAGATCCAAATTGGTGGTTGCGCCGCCGCTGCCGGGACGGGGGGCTATGAGGTTATGTTCGCGCTGGGGCTACGCGTCCCCGTTGTGACGGGGTAACTTTTTTTAAACCGGATCGGAGGCAACACCATCTGGGGGGaaaaaaa                                                                                                  >  2:533299/1‑239 (MQ=255)
                                                                                                                 tGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGTTGAGCTGGTGGGCAAGGGGATCAAAATTGATGATGTCGCCGCCGCTGCCGGGACGGGGGGGTATGTGGTGTTTTGCGGGCTGGCGCTACGCGTCCCGGGTGTGGCGGGGGAAATTGTTTTTATGCGGGACGGGGGGCACGGCTTCTGGGTGCaaaaaaaaaaaccc                                                                                           >  2:23635/1‑232 (MQ=255)
                                                                                                                 tGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCCGGTGAGCTGGGGGGCAAGGAGAGCAAAATTGTTGTTGTCGCCGCCGCTGCCGGGAAGGAGGGCTAGGAGTTTATTGTCCGGCTGGGCCTAGCCGGCCCGGGTGGGACGGGGGGACTTTTTGTTAGCCGGGCTCGAGGCAACCGCTTTTGGGGGGGAAAAAAAaatcc                                                                                            >  2:957579/1‑238 (MQ=255)
                                                                                                                  ggggCCCGTCGCGATGGATATGCTAGGGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATcatcca                                                                                           <  1:796644/238‑1 (MQ=255)
                                                                                                                   gggACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACCCCGGGTGAGCTGTGGGGCAAAGGGATCCAAATTGTTGATGTCGCCGCCCCTGCCGGGAGGGGGGGCTACGAGTTTATTGGCGCGCTGGCGTGACGCGGCCCGGTGGGTGCGGGGTTACTTTGTTGTACGCGGGGTGGGGGTCACCGCTCTCTGGGGACaaaaaaacaaaccac                                                                                          >  2:983549/1‑231 (MQ=255)
                                                                                                                   gggACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTAAGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGGTATGAGTTGATGTGCGCGCTGGCGCTACGCGGCCCGGGTGTGACGGGGTAAATTTTTGTAAACCGGATCTGAGGCAACGTCTTCTTGGGGCAAAAAAAAcatccat                                                                                          >  1:500059/1‑238 (MQ=255)
                                                                                                                   gggACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGGACGCCGGTTGAGGTGTGGGGCCAGGAGATCAAAATTGGTGTTGTCGCCGCCGCCGCCGGAACGGGGGGGTATGGGGTGGTGGTCACGCTGGCGCCCCCGGCCCCGGTTGGTGCGGGGTAACTTTTTTGAAACCGGGGCGGGGGCACGGCGTTTTGGGGGCaaaaaaaaaataccc                                                                                          >  2:1181027/1‑230 (MQ=255)
                                                                                                                       ccGTCTCGATGGATTTGCTGGCGGTCGGTTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCgg                                                                                      <  1:15699/238‑1 (MQ=255)
                                                                                                                            tcGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTc                                                                                <  1:922797/239‑1 (MQ=255)
                                                                                                                                    tctGCGAGCGGTCGATTTAACGCCTTGCCCGCAGGCGTGTATTGGGACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACt                                                                        <  1:772204/236‑1 (MQ=255)
                                                                                                                                       gCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGGATTGGTACGCCGGGTGAGCTGTGGGGCAAGGGGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGGGGGCTGTGAGTTGATGTGCGCGCTGGGGCTACGCGTCCCGGGTGGGCAGGGGTGACTTTTTTTAAACGGGAGCGGGGGGCACGGCTCTTGGGGGCAAAAAAAACATACATACCGCTGGGGAGCAACGGtgg                                                                    >  2:905932/1‑237 (MQ=255)
                                                                                                                                           gCGGTCGATTTAACGCCTTGCCCGCAGGGGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAgttgt                                                                <  1:319857/239‑1 (MQ=255)
                                                                                                                                                     tAACGGCTTGCCCGGAGGCGGGTATTGGTACGCCGGTTGAGCTGGGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACa                                                      <  1:14272/239‑1 (MQ=255)
                                                                                                                                                           ctgggcgggaGGCGGGTATTGGGACGCCGGTTGAGCTGGGGGGCACGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGGTGGCTATGAGTTGATGTGCGCGCTGGCGCTACGGGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCa                                                <  1:1241375/230‑1 (MQ=255)
                                                                                                                                                           cTTGCCCGGGGGGGGGTGGTGGTGCGGCGGTTGGGCTGGGGGGCAAGGAGGTCAAAATTGATGATGGCGCCGCCGCGGCCGGGACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCa                                                <  2:894877/239‑1 (MQ=255)
                                                                                                                                                           cTTGCCCGCGGGGGGTTTTTGGTACGGGGTTTGGGCTGGGGGGAAGGGGGATAAAAATGGATGATGTCGCCGCCGCTGCCGGAACGGTGGGGTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGGGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCa                                                <  2:371455/239‑1 (MQ=255)
                                                                                                                                                           cTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCa                                                <  1:905932/239‑1 (MQ=255)
                                                                                                                                                           cTTGCCCCGAGGGGGGTTTTGGTACGCGGGTTGAGCTGGGGGGCAAGGGGATAAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCa                                                <  2:824201/239‑1 (MQ=255)
                                                                                                                                                            ttGCCGGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGGCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCa                                                <  1:1205255/238‑1 (MQ=255)
                                                                                                                                                             tGCCCGCGGGCGGGTTTTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCAtt                                              <  1:23635/239‑1 (MQ=255)
                                                                                                                                                                ccGCGGGGGGGTTTTGGGTCGCCGGTTGGGCGGGGGGGCAGGGGGATCAAAATTGATGGTGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGc                                            <  2:776991/238‑1 (MQ=255)
                                                                                                                                                                caggaGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGcc                                           <  1:324321/235‑1 (MQ=255)
                                                                                                                                                                  gCAGGGGGGTATTGGTACGCCGGTTGAGGTGTGGGGCAAGTAGATCAAAATTGTTGCTGGCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCa                                         <  1:957579/239‑1 (MQ=255)
                                                                                                                                                                                   cGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAATGTCCATTaaa                        >  2:223807/1‑239 (MQ=255)
                                                                                                                                                                                    gCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTTTAAGCCGGATCGGAGGCACCGCCTTTGGGGGGCAAAAAAATCTTCAACCCGGTGGGCCGTACACGGAGGTTTTTAAGGTGAAAGAGCAATTGCCCATGGTATAGTCCATAaaaa                       >  2:1371628/1‑239 (MQ=255)
                                                                                                                                                                                    gCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAAAGCCATTGCCCATGATAGTGTACATTaaaa                       >  1:529056/1‑239 (MQ=255)
                                                                                                                                                                                    gCCGGTTGAGCTGTGGGGCAAGGAGATAAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTaaa                        <  1:1244113/238‑1 (MQ=255)
                                                                                                                                                                                     ccGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATAATATGGTCCATTAAAAg                      >  2:1359484/1‑239 (MQ=255)
                                                                                                                                                                                     ccGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATAATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGCAAGGGCCATTGCCCCTGATAGTGGCCATTAAAAg                      >  2:1244113/1‑239 (MQ=255)
                                                                                                                                                                                         tttggtgttgGGGAAAGGAGTACAAAATTGGTGGGGTCGCCGCCGCTTCCGGAACGGTGGGTTATGAGTTGATGTGCGCGCTGGCGCTACGGGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATTCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATg                  <  2:547168/231‑1 (MQ=255)
                                                                                                                                                                                           ggggtttgGGGCAAGGGATTTAAATTTGATGGGTTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAAATTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGa                <  2:235213/233‑1 (MQ=255)
                                                                                                                                                                                           gAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGa                <  1:223807/239‑1 (MQ=255)
                                                                                                                                                                                             gCTGTGGGGCAAGGAGATAAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGGCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCAGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGac               <  1:310620/238‑1 (MQ=255)
                                                                                                                                                                                              cTGTGGGGCAAGGAGATCAAAATTGAGGATGGCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGaca              <  1:553599/238‑1 (MQ=255)
                                                                                                                                                                                                gtcgggCAACGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGGGCAAAAAAATCATCCATCCGGCTGGCGAGCAACTGTAGTTGGTAATTGGACAGAACCATTGCCCATGACAAGGGCCATTAAAAGGAAGGgcacaa           >  2:654893/4‑237 (MQ=255)
                                                                                                                                                                                                        aGGAAAACAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTcccc    <  1:654893/238‑1 (MQ=255)
                                                                                                                                                                                                           agaTCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCCGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCgg  <  1:429192/237‑1 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                                
GCGCGCGATGAACAGCGAATCGGCATTGTCGCCGCAGGGTACGCCGACGGTTATCCGCGCCACGCGCCTACCGGTACCCCTGTTTTAGTGGACGGCGTGCGCACCATGACGGTGGGGACCGTCTCGATGGATATGCTAGCGGTCGATTTAACGCCTTGCCCGCAGGCGGGTATTGGTACGCCGGTTGAGCTGTGGGGCAAGGAGATCAAAATTGATGATGTCGCCGCCGCTGCCGGAACGGTGGGCTATGAGTTGATGTGCGCGCTGGCGCTACGCGTCCCGGTTGTGACGGTGTAACTTGTTGTAAGCCGGATCGGAGGCAACGTCTTCTGGGTGCAAAAAAATCATCCATCCGGCTGGTCAGCAACTGTAGTTGTTAATGTGACAGAGCCATTGCCCATGATAGTGTCCATTAAAAGGATGGACACTATTTCCCCGG  >  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: