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

AGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTG  >  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)
                            ggTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑Gc                                                                                                                                                                                                                     >  1:37444/1‑251 (MQ=255)
                            ggTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGCTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑Gc                                                                                                                                                                                                                     >  2:252836/1‑251 (MQ=255)
                                  cgctccGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTGACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTTTGTTGATCATCGTACGCAAGTGACCAACGCTGGCGATGGTGTCTTTGATGCCAACGATATTCTCGGTTTGCATAGCCAGAGGCGGCACGGGTGCCGGGGTTAATTCCTGACCCGTTTGATGCGGAAAGGTGTACTGTATCTCCGGTATGGttgcgttacgg                                                                                                                                                                                                                 <  2:206141/246‑11 (MQ=255)
                                      ccGTCGCCGCTCAGCAGCATCGTATTCAGCACATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGAt                                                                                                                                                                                                           >  1:130455/1‑251 (MQ=255)
                                         tCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGACTACCGCAAACGACCGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCACCGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATctg                                                                                                                                                                                                        >  2:29015/1‑251 (MQ=255)
                                          cgccgcTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATctg                                                                                                                                                                                                        <  1:182519/250‑1 (MQ=255)
                                               ctcagcagCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGt                                                                                                                                                                                                   >  1:240094/3‑250 (MQ=255)
                                                    cagcaTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAgt                                                                                                                                                                                             >  1:412651/1‑251 (MQ=255)
                                                       caTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAATGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAgtcgt                                                                                                                                                                                          >  1:229110/1‑251 (MQ=255)
                                                        aTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGCTTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAgtcgtc                                                                                                                                                                                         >  2:166624/1‑251 (MQ=255)
                                                         tCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAgtcgtc                                                                                                                                                                                         >  1:113113/1‑250 (MQ=255)
                                                         tCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGGTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAgtcgtc                                                                                                                                                                                         >  1:283696/1‑250 (MQ=255)
                                                              ttCAGCAAATGATCATCGTAACCGCAGAATACCGCAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAgg                                                                                                                                                                                                                             >  2:281754/1‑210 (MQ=255)
                                                              ttCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAgg                                                                                                                                                                                                                             <  1:281754/210‑1 (MQ=255)
                                                               tCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGAGTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGAtt                                                                                                                                                                                   <  1:194425/250‑1 (MQ=255)
                                                                   cAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGTCTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCgtg                                                                                                                                                                              <  2:112301/251‑1 (MQ=255)
                                                                       tgctcatcGTAACCGCAGAATACCGAAAACGAAGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCTCCGGTAGGGAGAC‑GCTACGGGCGATCTGCTGGTAAGAGTCGTCAAGATTTCTTGGTGc                                                                                                                                                                          >  2:106363/4‑251 (MQ=255)
                                                                           catcGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTAATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTGAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGAc                                                                                                                                                                       <  2:162536/250‑1 (MQ=255)
                                                                               gTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑g                                                                                                                                                                                                                      <  1:10139/199‑1 (MQ=255)
                                                                               gTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑g                                                                                                                                                                                                                      >  2:10139/1‑199 (MQ=255)
                                                                                 aaCCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTGTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCag                                                                                                                                                                <  1:78791/251‑1 (MQ=255)
                                                                                    cGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCagt                                                                                                                                                               >  2:104241/1‑249 (MQ=255)
                                                                                        gAATTCCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCTACTTTCCAGTAGTAgg                                                                                                                                                         <  2:294243/251‑1 (MQ=255)
                                                                                                aaaaCGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATg                                                                                                                                                 >  2:301890/1‑251 (MQ=255)
                                                                                                        gggCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGAt                                                                                                                                         <  2:349446/251‑1 (MQ=255)
                                                                                                        aggCGTACCGACTTAACTGTGATGATCTTCGTACGCAAGTGCCCAACGCTGTCGATGGAGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCCCGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGCAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGCTTTCGTGGTGCGACTTTCCAGTAGTCGGGGTTGATGGGGACGAt                                                                                                                                         >  2:422620/2‑251 (MQ=255)
                                                                                                              aCCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTCCGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATc                                                                                                                                   <  1:183892/251‑1 (MQ=255)
                                                                                                                     tAACTGTGTTGATCATCGTACGCAAGTTACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCGCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGcc                                                                                                                             <  2:372053/250‑1 (MQ=255)
                                                                                                                       aCTGTGTTGAGCATCGTACGCAAGTGACCAACGCTGTCGGTGGTGTCTTTGATGCCAACGATATTCTCGTTTAGCAGAGCCAGACGCGTCTCGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACTGC‑ACGGGCGATATTCTGGTAATAGGCGTCAAGATTTCGTGCGGCGCCTCTCCAGTAGCAGGGGTTGAGGTCGAGGATGCCATCAGCGCCGTa                                                                                                                          <  2:40857/251‑1 (MQ=255)
                                                                                                                          tggtTGATCATCGTACTCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAg                                                                                                                         <  1:389704/247‑1 (MQ=255)
                                                                                                                                atcatcGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTgcgc                                                                                                                 <  2:218249/251‑1 (MQ=255)
                                                                                                                                  catcGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTACCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCa                                                                                                                <  1:166624/250‑1 (MQ=255)
                                                                                                                                      gTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCACCGATATTCTCGTTTTGCAGAGCCAGACGCGTCAAGGTGTCCGGGGTTAATTCCTGACCCGTCAGATTCGGACAGTTGTACAGGATCACCGGTCCGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTACCTGAAGTAGGGGTTGATGCCGACGTGACCATCAGCGCAGTAGGCTTGCGCAtgctt                                                                                                           >  1:410002/1‑250 (MQ=255)
                                                                                                                                        aCGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTgcg                                                                                                         <  1:257118/251‑1 (MQ=255)
                                                                                                                                            aaGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACAATCAGCGCCGTAGGCTTGCGCATGCTGCGcc                                                                                                       >  2:249607/1‑249 (MQ=255)
                                                                                                                                                       gCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCtt                                                                                          <  1:31667/251‑1 (MQ=255)
                                                                                                                                                                 gtgtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGa                                                                                <  2:332804/251‑1 (MQ=255)
                                                                                                                                                                  tgtCTTTGATGCCAACGATATTCTCGGTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTACTGACCCGTCAGATTCGGAAAGTTGAACAGGATAACCAGTGGGGAAAC‑TCTCCGGGCGATCTGCTAGGAATAGGCGTAAAAATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGTCGACGATACCATCAGCGCCGTCGTCTTGAGCATGCTGCGCCAGATTGatcactccgacagtgta                                                                                >  1:204278/1‑234 (MQ=255)
                                                                                                                                                                    tCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGacgg                                                                             >  1:444546/1‑248 (MQ=255)
                                                                                                                                                                     cgttGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGGTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAgg                                                                             <  2:141318/248‑1 (MQ=255)
                                                                                                                                                                            gCCAACGATATTCTCGTTTGGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCGCCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGGCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCg                                                                      <  1:422620/250‑1 (MQ=255)
                                                                                                                                                                                 cGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGcc                                                                 >  1:400275/1‑250 (MQ=255)
                                                                                                                                                                                          cGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCcaatcaat                                                          <  2:37444/248‑1 (MQ=255)
                                                                                                                                                                                           gTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTTACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCg                                                      <  1:431908/251‑1 (MQ=255)
                                                                                                                                                                                               tGCAGCGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCTTCAGATTCGGAAAGTTGTCCAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATCCCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGTTTCGTCAGTGGAAGGTTAACCGACGCCTATCAATACCGGCAc                                                  >  2:455771/1‑251 (MQ=255)
                                                                                                                                                                                                    agCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATa                                                                                                                                        >  1:78362/1‑160 (MQ=255)
                                                                                                                                                                                                    agCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATa                                                                                                                                        <  2:78362/160‑1 (MQ=255)
                                                                                                                                                                                                     gCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGCTTCGGAAAGTTGTACAGGATCCCCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGccc                                            >  2:394117/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGtt                               >  2:33374/1‑251 (MQ=255)
                                                                                                                                                                                                                         cGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGACGACAATGGTTACAGctt                        >  1:317450/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCg                <  2:44281/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGt             <  2:412651/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTg            >  2:141366/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cccGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCAt        >  1:351763/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGCCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATgcg     >  1:196285/1‑251 (MQ=255)
                                                                                                                                                                                                                                               aGATTCGGAAAGTGGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGGCAAGATTGCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGGGAGGGCCATGCGCTg  <  2:229110/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
AGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGAC‑GCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTG  >  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: ‑