Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,524,596 C→T D139N (GAT→AAT)  yjhH ← putative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,524,5960CT98.4% 208.3 / NA 64D139N (GAT→AAT) yjhHputative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/31);  minor base G (0/1);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTG  >  NC_000913/4524353‑4524842
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
agagTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCCGCAGCATCGTATTCCGCCAATGATCATCGCAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGGTCATCGTACGCCAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCCACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGaa                                                                                                                                                                                                                                                    >  1:218249/1‑249 (MQ=255)
  agTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGtt                                                                                                                                                                                                                                                <  2:352354/251‑1 (MQ=255)
  agTTCCGGAGCAAAGTTAGCGCTGGCGGTTATAGCAACGTCGCCGCTCAGCAGCATCGTATTCAGCAAAAGATCATCGTAACCGCAGAATACCGAACACGGCGGGCGTAACGACTTAACTGTGTTGGTCATCGTACGCAAGTGACCAACGCTGTCGATGGGGTCTTTGATGCCAACGGTATTCTCGTTTTTCAGAGCCAGACGCGTCCCGGTTTCCGGGGTTAAGTCCTTGCCCGTCAAATTCGGAAAGtt                                                                                                                                                                                                                                                >  2:314551/1‑251 (MQ=255)
   gTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTg                                                                                                                                                                                                                                               <  1:390215/251‑1 (MQ=255)
           gCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAgggt                                                                                                                                                                                                                                       >  1:340984/1‑249 (MQ=255)
                   aGCGCTGGCGTTTAGCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACTACGGGCGTACCGACTTAACTGTGATGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTGTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACTGTTTCCGGGGTTAAGTCCTGACCCGTCAGTTTCGGAAAGTTGTACAGGATCACCGGt                                                                                                                                                                                                                                <  2:199794/250‑1 (MQ=255)
                   aGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTGTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTa                                                                                                                                                                                                                               <  2:150405/251‑1 (MQ=255)
                          gCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAc                                                                                                                                                                                                                        <  1:314551/251‑1 (MQ=255)
                            ggTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGc                                                                                                                                                                                                                     >  1:37444/1‑251 (MQ=255)
                            ggTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGCTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGc                                                                                                                                                                                                                     >  2:252836/1‑251 (MQ=255)
                                  cgctccGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTGACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTTTGTTGATCATCGTACGCAAGTGACCAACGCTGGCGATGGTGTCTTTGATGCCAACGATATTCTCGGTTTGCATAGCCAGAGGCGGCACGGGTGCCGGGGTTAATTCCTGACCCGTTTGATGCGGAAAGGTGTACTGTATCTCCGGTATGGttgcgttacgg                                                                                                                                                                                                                 <  2:206141/246‑11 (MQ=255)
                                      ccGTCGCCGCTCAGCAGCATCGTATTCAGCACATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGAt                                                                                                                                                                                                           >  1:130455/1‑251 (MQ=255)
                                         tCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGACTACCGCAAACGACCGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCACCGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATctg                                                                                                                                                                                                        >  2:29015/1‑251 (MQ=255)
                                          cgccgcTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATctg                                                                                                                                                                                                        <  1:182519/250‑1 (MQ=255)
                                               ctcagcagCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGt                                                                                                                                                                                                   >  1:240094/3‑250 (MQ=255)
                                                    cagcaTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgt                                                                                                                                                                                             >  1:412651/1‑251 (MQ=255)
                                                       caTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAATGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgtcgt                                                                                                                                                                                          >  1:229110/1‑251 (MQ=255)
                                                        aTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGCTTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgtcgtc                                                                                                                                                                                         >  2:166624/1‑251 (MQ=255)
                                                         tCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgtcgtc                                                                                                                                                                                         >  1:113113/1‑250 (MQ=255)
                                                         tCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGGTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgtcgtc                                                                                                                                                                                         >  1:283696/1‑250 (MQ=255)
                                                              ttCAGCAAATGATCATCGTAACCGCAGAATACCGCAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAgg                                                                                                                                                                                                                             >  2:281754/1‑210 (MQ=255)
                                                              ttCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAgg                                                                                                                                                                                                                             <  1:281754/210‑1 (MQ=255)
                                                               tCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGAGTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGAtt                                                                                                                                                                                   <  1:194425/250‑1 (MQ=255)
                                                                   cAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGTCTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCgtg                                                                                                                                                                              <  2:112301/251‑1 (MQ=255)
                                                                       tgctcatcGTAACCGCAGAATACCGAAAACGAAGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCTCCGGTAGGGAGACGCTACGGGCGATCTGCTGGTAAGAGTCGTCAAGATTTCTTGGTGc                                                                                                                                                                          >  2:106363/4‑251 (MQ=255)
                                                                           catcGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTAATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTGAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGAc                                                                                                                                                                       <  2:162536/250‑1 (MQ=255)
                                                                               gTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACg                                                                                                                                                                                                                      <  1:10139/199‑1 (MQ=255)
                                                                               gTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACg                                                                                                                                                                                                                      >  2:10139/1‑199 (MQ=255)
                                                                                 aaCCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTGTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCag                                                                                                                                                                <  1:78791/251‑1 (MQ=255)
                                                                                    cGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCagt                                                                                                                                                               >  2:104241/1‑249 (MQ=255)
                                                                                        gAATTCCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCTACTTTCCAGTAGTAgg                                                                                                                                                         <  2:294243/251‑1 (MQ=255)
                                                                                                aaaaCGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATg                                                                                                                                                 >  2:301890/1‑251 (MQ=255)
                                                                                                        gggCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGAt                                                                                                                                         <  2:349446/251‑1 (MQ=255)
                                                                                                        aggCGTACCGACTTAACTGTGATGATCTTCGTACGCAAGTGCCCAACGCTGTCGATGGAGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCCCGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGCAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGCTTTCGTGGTGCGACTTTCCAGTAGTCGGGGTTGATGGGGACGAt                                                                                                                                         >  2:422620/2‑251 (MQ=255)
                                                                                                              aCCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTCCGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATc                                                                                                                                   <  1:183892/251‑1 (MQ=255)
                                                                                                                     tAACTGTGTTGATCATCGTACGCAAGTTACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCGCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGcc                                                                                                                             <  2:372053/250‑1 (MQ=255)
                                                                                                                       aCTGTGTTGAGCATCGTACGCAAGTGACCAACGCTGTCGGTGGTGTCTTTGATGCCAACGATATTCTCGTTTAGCAGAGCCAGACGCGTCTCGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACTGCACGGGCGATATTCTGGTAATAGGCGTCAAGATTTCGTGCGGCGCCTCTCCAGTAGCAGGGGTTGAGGTCGAGGATGCCATCAGCGCCGTa                                                                                                                          <  2:40857/251‑1 (MQ=255)
                                                                                                                          tggtTGATCATCGTACTCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAg                                                                                                                         <  1:389704/247‑1 (MQ=255)
                                                                                                                                atcatcGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTgcgc                                                                                                                 <  2:218249/251‑1 (MQ=255)
                                                                                                                                  catcGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTACCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCa                                                                                                                <  1:166624/250‑1 (MQ=255)
                                                                                                                                      gTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCACCGATATTCTCGTTTTGCAGAGCCAGACGCGTCAAGGTGTCCGGGGTTAATTCCTGACCCGTCAGATTCGGACAGTTGTACAGGATCACCGGTCCGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTACCTGAAGTAGGGGTTGATGCCGACGTGACCATCAGCGCAGTAGGCTTGCGCAtgctt                                                                                                           >  1:410002/1‑250 (MQ=255)
                                                                                                                                        aCGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTgcg                                                                                                         <  1:257118/251‑1 (MQ=255)
                                                                                                                                            aaGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACAATCAGCGCCGTAGGCTTGCGCATGCTGCGcc                                                                                                       >  2:249607/1‑249 (MQ=255)
                                                                                                                                                       gCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCtt                                                                                          <  1:31667/251‑1 (MQ=255)
                                                                                                                                                                 gtgtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGa                                                                                <  2:332804/251‑1 (MQ=255)
                                                                                                                                                                  tgtCTTTGATGCCAACGATATTCTCGGTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTACTGACCCGTCAGATTCGGAAAGTTGAACAGGATAACCAGTGGGGAAACTCTCCGGGCGATCTGCTAGGAATAGGCGTAAAAATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGTCGACGATACCATCAGCGCCGTCGTCTTGAGCATGCTGCGCCAGATTGatcactccgacagtgta                                                                                >  1:204278/1‑234 (MQ=255)
                                                                                                                                                                    tCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGacgg                                                                             >  1:444546/1‑248 (MQ=255)
                                                                                                                                                                     cgttGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGGTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAgg                                                                             <  2:141318/248‑1 (MQ=255)
                                                                                                                                                                            gCCAACGATATTCTCGTTTGGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCGCCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGGCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCg                                                                      <  1:422620/250‑1 (MQ=255)
                                                                                                                                                                                 cGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGcc                                                                 >  1:400275/1‑250 (MQ=255)
                                                                                                                                                                                          cGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCcaatcaat                                                          <  2:37444/248‑1 (MQ=255)
                                                                                                                                                                                           gTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTTACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCg                                                      <  1:431908/251‑1 (MQ=255)
                                                                                                                                                                                               tGCAGCGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCTTCAGATTCGGAAAGTTGTCCAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATCCCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGTTTCGTCAGTGGAAGGTTAACCGACGCCTATCAATACCGGCAc                                                  >  2:455771/1‑251 (MQ=255)
                                                                                                                                                                                                    agCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATa                                                                                                                                        >  1:78362/1‑160 (MQ=255)
                                                                                                                                                                                                    agCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATa                                                                                                                                        <  2:78362/160‑1 (MQ=255)
                                                                                                                                                                                                     gCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGCTTCGGAAAGTTGTACAGGATCCCCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGccc                                            >  2:394117/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGtt                               >  2:33374/1‑251 (MQ=255)
                                                                                                                                                                                                                         cGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGACGACAATGGTTACAGctt                        >  1:317450/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCg                <  2:44281/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGt             <  2:412651/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTg            >  2:141366/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cccGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCAt        >  1:351763/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGCCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATgcg     >  1:196285/1‑251 (MQ=255)
                                                                                                                                                                                                                                               aGATTCGGAAAGTGGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGGCAAGATTGCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGGGAGGGCCATGCGCTg  <  2:229110/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
AGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTG  >  NC_000913/4524353‑4524842

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