Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,318,215 G→A L184L (CTG→TTG)  phnK ← carbon‑phosphorus lyase subunit PhnK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,318,2150GA93.8% 162.7 / ‑3.8 65L184L (CTG→TTG) phnKcarbon‑phosphorus lyase subunit PhnK
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (44/16);  total (45/20)
Fisher's exact test for biased strand distribution p-value = 7.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGC  >  NC_000913/4317968‑4318457
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
aGGCCGGATAAGGCGGTTACGCCGCATCCGGCACCGTCTCAATGCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCg                                                                                                                                                                                                                                                 <  1:300509/251‑1 (MQ=255)
  gctggATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCgtg                                                                                                                                                                                                                                               <  1:238272/248‑3 (MQ=255)
   ccGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGc                                                                                                                                                                                                                                              >  1:31770/1‑251 (MQ=255)
     ggATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCt                                                                                                                                                                                                                                            <  1:445075/251‑1 (MQ=255)
             cGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGCGCGACGCCTAAATCATGGGTGACAATCACCACCGCGGGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAATCGGGCCTCCACCGac                                                                                                                                                                                                                                    >  1:9110/1‑251 (MQ=255)
                  acgcagtatctggaaccgagtCAATTCGGCAAATCCGAGGCCACCAGCAGCTGTGTATAAGGATGATGCGGGTCGTCGAGTACGTGGTCGGTTAACCCACTCCCCACCCATTGCCCCTGCTTCATCCCCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCCCCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCTGGTCGAGCAAGCGGGCCTGCACCGACACATc                                                                                                                                                                                                                               <  2:363813/231‑1 (MQ=255)
                      cgcATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGCTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCTACACATCCAGc                                                                                                                                                                                                                           >  1:50301/1‑251 (MQ=255)
                      agcATCCGGCACCGGCACACTTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGAAGATGCGGGTCGTCGAGCACGCGGCCGGGTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACCCCTAAATCTTGGGGCGCAAGCACGACCGCGAGGTTCAGCTCCCCCACCACGCCCCGCCGCGGGTCGAGCCAGCGGGCCCGCAcctaccactccagc                                                                                                                                                                                                                           >  2:72524/2‑239 (MQ=255)
                        cATCCGGCACCGGCTCAATTCTGCAAAACCGTTGACACCAGCAGCTCTGCATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCTCTTGCCCCTTCTTCCTCACCAGCAAACGGACTGCCAGCAGGCGGGAGACGCCTCAATCATGCGTGACAAACACCACCACGCAGGGCAGCTCCACCACCACGCCGCCACACCGGTCGAGCACGCGGGCCCGCCCCGACCCATCCAcccc                                                                                                                                                                                                                         >  2:288183/1‑251 (MQ=255)
                               cACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCACACGGTCCGCCCGCAGGCGGGCGACGCCTAAATCATGTGTGACAATCACCACCGCGAGGTTCAGCACCACCACCAGCCCGCGCAGCAGGTCGAGCAAGCGGGCCATAACCGGCACAACCAGCCCACcggtc                                                                                                                                                                                                                  >  2:380825/1‑251 (MQ=255)
                                                aaaaCCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGCCGCCTAAATCCTGGGTGACAATCACCACCACGAGGTTCAGCCCCACCACCAGGCCACGCAGCCGGGCGAACAACCGGGCCGGCACCGGCCCATCCACCCCGCCGGTCGGTTCATCCCTAAccac                                                                                                                                                                                                 >  2:398032/1‑251 (MQ=255)
                                                    ccGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGATCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAACCACCAGc                                                                                                                                                                                             >  1:297589/1‑251 (MQ=255)
                                                     cGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCt                                                                                                                                                                                            >  1:305595/1‑251 (MQ=255)
                                                       aTGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTc                                                                                                                                                                                          >  1:96147/1‑251 (MQ=255)
                                                         gACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCgg                                                                                                                                                                                        >  1:364094/1‑251 (MQ=255)
                                                                    cTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGATCAGCTCCACCACCAGGCCCCGCAGCCGGTCGCGCAAGCGGGCCTGCACCGACACATCCCACCCGCCGGCCGGGTCATCCCTAAACACCCGCTTCGGATGcagccccc                                                                                                                                                                             >  2:337365/1‑244 (MQ=255)
                                                                                   atgCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGCCACATCCAGCCCGCCGGTCGGTTCATCCATAAACCCCAGCTTCGGATGCGTCCCCCGGTTGCGGGCaaacc                                                                                                                                                              >  2:214032/1‑248 (MQ=255)
                                                                                    cggctctccctacctccgcgGTCGGTTAACCAACTCTCCACCAGTGGCCCTGGCTTCTCCATCACAAAACGCCCCGCCGGCAGGGGGGCGACCCCAAAATCTGGGGGGTCAACCACCCCCGCGAGTTTCAGTCCCCCCACCAGGCCCAGCAGCAGTTCGAGCAAGCGGGCCTCCACCTACCCTTCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTGCGGATGCGTCACCAGGTTGCGGGCAATCTg                                                                                                                                                             <  2:399698/235‑1 (MQ=255)
                                                                                      cGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCCTACACCCCAGCTTCGGATGCGTCCCCAGGTTGCGGGCCATCTGc                                                                                                                                                            >  2:317097/1‑250 (MQ=255)
                                                                                       gggTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGGGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCCCCGACACATCCACCCCGCCGGTCGGTTCATCCCTAAACACCACCTTCGGATGCGTCACCAGGGTGCGGGCAATCTGCac                                                                                                                                                          >  2:273893/1‑250 (MQ=255)
                                                                                               gAGCACGTGGTCGGTTAACCCACCCCCCACCCCTTGCCCCGGCTTCATCACCAGCAAACGGCCCGCCGGCAGGCGGGCGACGCCTAAATAATGGGGACCAATCCCCCCCGCGAGGTTCATCTCCACCACCAGGCCGCGCAGCAGGCCGAGCAGGCGGGCCTGCCCCGACACTTCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACgctgct                                                                                                                                                  <  1:54330/251‑1 (MQ=255)
                                                                                                   aCGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGGGACAATCACCACCGCGAGGTTCAACTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAACGGGCCTTCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGGCACCAGGTTGCGGGCCATCTg                                                                                                                                                             >  1:399698/1‑236 (MQ=255)
                                                                                                      cggtcggtTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAACCGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGCCACATCCAGCCCGCCCGTCGCTTCATCCATAACCACCAGCTTCGGATCCGTCACCAGGGTGCGGGCAATCTGAAAACGCTGCGGCACAcc                                                                                                                                           >  2:155677/1‑251 (MQ=255)
                                                                                                          gggttcaccccctccccccacctctccccTGCTACATCCCCGGCAACGGGGCCGCCGGGAGGCGGGGGCCGCCTAAATCATGGGTGACAATCCCCACCTGGAGGTTCAGCCCCACCACCAGGCCGCGCAGCGGGCCGAGCAGGCGTGCCTGCCCCGACACATCCAGCCCGCCGGTCGGTTCACCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATAccgccg                                                                                                                                       <  2:96070/226‑1 (MQ=255)
                                                                                                          cggtTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGCCACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCCTCGGATGCGTCACCCGGTTGCGGGCCATCTGCAACCGCTGCTGCATAccaccg                                                                                                                                       >  2:121657/1‑251 (MQ=255)
                                                                                                          cggatAACCCACTCTCCACCACTTGCCCCGGCTTCATCACCAGCAAACGGTCCGCCAGCACGCGGACGACGCCCAACTCATGGGCCACAATCACCACCGCCCGGATCAGCTCCACCACCAGGCCGCGCAGCAGGTCGACCAAGCGGGCCTGCACCGACACATCCAACCCGCCGGTCGGTTCACCCATACACACCAGATTCGGATTCGCCCCCAGGTTTCGGGATATCTGCAAACCCTCCTGCacccccccg                                                                                                                                       >  2:407213/5‑243 (MQ=255)
                                                                                                               aaCCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCCTAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCCATCTGCAACCGCTGCTGCATACCGCCggcaa                                                                                                                                   >  2:385969/1‑247 (MQ=255)
                                                                                                                  cctctcccctccacgttccccTGCTTCAACCCCAGCAAACGGTCCCCATGCAGGCGGGCGACTCTAAAATCATGGTTGACAATCACCACCGCGAGGTTCAGCTCCCCCCCCAGGCGGTGCAGCGGGTCGAGCAAGCGGGCCTGCCCCGACACGTCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAggt                                                                                                                               <  1:69869/234‑1 (MQ=255)
                                                                                                                    tctctCCACCACGTGCCCCTGCTTCATCCCCAGCAAACGGTCCGCCAGCGGGCGGGCGACGCCTAAATCAGGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACAGCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCTGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAgg                                                                                                                                <  1:429984/247‑1 (MQ=255)
                                                                                                                    cctctCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCCCCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAggtgg                                                                                                                             <  1:209919/250‑1 (MQ=255)
                                                                                                                    aCTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCCCCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAACGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCTTCCATAAACACCAGCTTCGGATGCGTAACCCGGTTGCGGGCCATCTGCAAACGCTGCTGCCTACCGCCGGGAAAggggg                                                                                                                             >  2:429984/1‑248 (MQ=255)
                                                                                                                    aCTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAggtgg                                                                                                                             >  1:241341/1‑251 (MQ=255)
                                                                                                                       cccccccccggcccccTGTTTCACCCCCAGCAACCGGCCCGCCAGCGGGCGGGCGACGCCTAAATCAGGGGTGACAATCCCCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGAAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCg                                                                                                                          <  2:448384/239‑1 (MQ=255)
                                                                                                                                   ccccTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGtcgtcgg                                                                                                              >  1:10915/1‑250 (MQ=255)
                                                                                                                                          ttCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGGTCATCCATAAACACCAGCTTCGGATGCGGCACCAGGTTGCGGGCAATCTGCAAAAGCTGCTGCATACCGCCGGCAAAGGGTGTCGGCAGGGCGGCGATCCGGtt                                                                                                       >  2:97182/1‑251 (MQ=255)
                                                                                                                                            catcaCCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCTTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTgg                                                                                                     <  1:195346/251‑1 (MQ=255)
                                                                                                                                                aCCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGGGGTCGGCAGGTCGTCGATCCGGTTGGCggg                                                                                                 >  1:279309/1‑251 (MQ=255)
                                                                                                                                                 ccAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGATCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACAAGCTTCGGATGCGTAACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAATGGTGGACGGCAGGTCGTCGATCCGGTTGGCggg                                                                                                 >  1:238172/1‑250 (MQ=255)
                                                                                                                                                   aGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCCCCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATACACACCCGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCCGCTGCATACCGCCGGAAAAGGTGGTCGGCCGGTCGTCGACCCGGTTGGCGGGAAt                                                                                              >  2:213446/1‑251 (MQ=255)
                                                                                                                                                       aaCGGTCCGCCAGCAGGTGGGCGCCGCCTAATTCATGGGTGACAATCACCCCCGCGAGGTTCACCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTcc                                                                                          <  1:127046/251‑1 (MQ=255)
                                                                                                                                                       aaCGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTcc                                                                                          >  1:151568/1‑251 (MQ=255)
                                                                                                                                                          ggTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAcc                                                                                       >  1:150876/1‑251 (MQ=255)
                                                                                                                                                             ccGCCAGCAGGCGGGCGACGCCTAAATCAGGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctct                                                                                    <  1:150782/251‑1 (MQ=255)
                                                                                                                                                               gCCAGCAGGCGGGCGACGCCTAACTCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACAACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCCCCGACACATCCAGCCCGCCGGTCGGTTCCTCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCACACGCTGCTGCCTACAGCCGGAAAAGGTGGTCGGCAGGGCGTCGATCCGGTTGGCGGGacactccaactcgtc                                                                                  >  2:485426/1‑237 (MQ=255)
                                                                                                                                                                   aacagGCGGCCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCCGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCACGCCCGCCGGTCGGTTAATCCATCAACACCAGCTTAGGATGCGTCACCAGGTTGCGGGCAATCTGCACACGCTGCTGCATACCGCCGTCAGAGGTGGTCGGCAGGTCGTCGCACCGGTTGGCGGGAGTCGCCAccccctacac                                                                               >  2:222977/3‑243 (MQ=255)
                                                                                                                                                                               aCGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCCGGTCGTCGATCCGGTTGGCGGGGATCTCCACCTCTTCCAGCCACTTCTgcgc                                                                   >  2:191297/1‑251 (MQ=255)
                                                                                                                                                                                cGCCTAAATCATGGGTGACAATCACCACCGCGAGTTCCAGCTCCCCCACCAGGCCGCGCAGCAGTTCGAGCAGGCGGGCCTGCACCGACCCATCCACCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGcgcg                                                                  <  1:84164/251‑1 (MQ=255)
                                                                                                                                                                                           tGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCAGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAat                                                       >  1:350780/1‑251 (MQ=255)
                                                                                                                                                                                           ggggTGACAATCCCCCCCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAat                                                       <  1:118832/250‑1 (MQ=255)
                                                                                                                                                                                            gggTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAat                                                       >  1:150192/1‑250 (MQ=255)
                                                                                                                                                                                                      ccccccccGCGAGGTTCAGCTCCACCCCCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCCCCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTaa                                             <  1:128748/245‑1 (MQ=255)
                                                                                                                                                                                                           accGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACg                                        >  1:214513/1‑251 (MQ=255)
                                                                                                                                                                                                            ccGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACCCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGt                                       >  1:348560/1‑251 (MQ=255)
                                                                                                                                                                                                              gcgAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGc                                     >  1:142709/1‑251 (MQ=255)
                                                                                                                                                                                                               cgAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGcc                                    >  1:219181/1‑251 (MQ=255)
                                                                                                                                                                                                                 gaggttaagacccaccaccaGGCCGCGGAGAAGGTCGAGCAAGGGGGCCTGCACCGACACTCCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGccc                                   <  1:320265/240‑1 (MQ=255)
                                                                                                                                                                                                                  ggtttCAGCCCCCCCCCCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCCCCGCCACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGcccc                                  <  1:103360/248‑1 (MQ=255)
                                                                                                                                                                                                                          cTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGcc                           >  1:21463/1‑250 (MQ=255)
                                                                                                                                                                                                                           tCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCAt                         >  1:236901/1‑251 (MQ=255)
                                                                                                                                                                                                                           cccaccaccaGGCCGCGCCGCAGGTCGAGCAAGCGGGCCTGCCCCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCAt                         <  1:381328/250‑1 (MQ=255)
                                                                                                                                                                                                                                    aGGCCGCGCAGCAGGTCGTGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGCAAAGGTGGTCGGCAGGTCGTCGGTCCGGTTGGCGGGAATCTCCACCTCTTCCAACCACTTCTGCGAGGTGGCCCGAATATAGCCGTAATGACGTGCCCCGGTCGCTATCACCCGCTc                >  2:356958/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGGCGATCCGGTTGGCGGGAAGCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCCCGAATATCGCCGGAATGACGGGCCCCGGTCGCCATCacccccccc               >  2:430738/1‑243 (MQ=255)
                                                                                                                                                                                                                                             agcagGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGGGCCCCGGGCGCCATCAGCCCCCCGCCGATAtt       >  2:341809/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gcagGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTg      >  2:27384/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   gTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGGGCCCCGGTCGCCATCAGCCGCTCGCCGATATTgccgc  >  1:373746/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
AGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGC  >  NC_000913/4317968‑4318457

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