Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913815,8400TC26.6% 116.5 / 10.3 66I270V (ATT→GTT) ybhKputative transferase YbhK
Reads supporting (aligned to +/- strand):  ref base T (34/13);  new base C (0/17);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 8.68e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGTATAAAAACTGCCAGG  >  NC_000913/815598‑816080
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
cccttaaccggCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTGGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACGGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGctct                                                                                                                                                                                                                                           <  2:172915/241‑1 (MQ=255)
                   ttCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCAc                                                                                                                                                                                                                       >  1:307798/1‑251 (MQ=255)
                                ttggAATTTGCGGGGTTTGTGGCCCGGATAAGGGGTTCCCGCCGCACCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTCCCTTCGGATTCCCGTTACCCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCGGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCcgac                                                                                                                                                                                                           <  1:214904/247‑1 (MQ=255)
                                  tGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGACTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGAt                                                                                                                                                                                                        >  1:569534/1‑251 (MQ=255)
                                           gcgGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGCCATGACGATACGGAATATCGCTGGCCTCCAGTACCCCCTGGATCACAATCCGCCCTTTCACCGCCGCGACATCCACTTTTGGCCCGACGATGAcaccaac                                                                                                                                                                                               >  2:128507/1‑245 (MQ=255)
                                                  gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGAcc                                                                                                                                                                                        >  2:530167/1‑251 (MQ=255)
                                                   taagCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCt                                                                                                                                                                                       >  1:456027/4‑251 (MQ=255)
                                                         ggATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTAc                                                                                                                                                                                 >  1:237952/1‑251 (MQ=255)
                                                          gATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTAcc                                                                                                                                                                                >  1:117184/1‑251 (MQ=255)
                                                                 ggtccccgccgcATCCGGCATAAACACCGCCCCCTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTCCCGTTAACCTAAAGCCGGTAACGCCTTTCCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCGGGCCTCCATTACCCCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATAc                                                                                                                                                                         <  1:8129/245‑1 (MQ=255)
                                                                  gTTCACGCCGCATCCGGCATAAACAACGCGACCTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTACCCTAAAGCCTGAAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCAGGACGATACGGAATATCGCGGGCCTCCAGTACCTCCTGGATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACt                                                                                                                                                                        <  1:187996/251‑1 (MQ=255)
                                                                        gcctcATCCAGCATAAACAACGCGCACTTTGTATTCAATAGGCACCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTCACAACTGGCGGTCATGACGATACGGTAAATCGCTGGCCTCCAGTACCTCCTGGTTCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGAAGTCGACAGCATCCATGACCTCTTTACCAACATACttcccc                                                                                                                                                                   >  2:48509/5‑245 (MQ=255)
                                                                         ccgcaccgcgaaTAAACACCGCTCACTTTCTCTTCAAACGGACCCCGAAATTTTCCTTCGATTTTCCGTAACCCAAAAGCCTTAACCCCCTTTACCAGCGCTTTAGGTAACACCGGCCGTTCATGACGATCCGGATTATCGCTGCCCCCAGGTACCCCTGGGACCCCATTCCCCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCAt                                                                                                                                                                 <  2:139383/240‑1 (MQ=255)
                                                                         ccccccccgtaaTAAACAACGCGTCTTTGGTCTCCAATGGGAAGCTGGATTTTTCCTCCGGATTCCCTTAAACCAAACGCCTGTAACGCTTTTCCCCCCGCGTAATGAAACAACGGGCGGTCAGGAGCATCCGAAGTATCGCGGGCCCCCATTCCCTCCTGGACCCCAACCCCCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCAt                                                                                                                                                                 <  2:237952/240‑1 (MQ=255)
                                                                            cATCCGGCATAAACACCGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTACCCTAAAGCTTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTCCCTCCTGGATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAAt                                                                                                                                                              <  1:530167/251‑1 (MQ=255)
                                                                              cccGGCATAAACAACGCGCTCTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCTTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTCCCTCCTGGATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTg                                                                                                                                                            <  1:168271/250‑1 (MQ=255)
                                                                              cccGGCATAAACAACGCGCCCTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTACCCTAAACCCTGTAACGCCTTTTCCAGCGCGTTATTTAACACCTGGCGGTCAGGACGATACGGAATATCGCTGGCCTCCAGTCCCTCCTGGTTCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTg                                                                                                                                                            <  1:25057/250‑1 (MQ=255)
                                                                                cGGCATAATCAACGCGCACTTTGTTTTAAATAGGACGCAGGAATTTTCCTTCGGTTTTCCGTTAACCTAAACCCTGTCACGCTTTTTCCAGCGCGTTATGTACCAACTGGCGGTCATGACTAGTCGGAATATCGCTGCCCCCCAGTCCCCCCTTGTTCCCATCCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGcc                                                                                                                                                          <  1:85050/251‑1 (MQ=255)
                                                                                    aTAAACAACGCGCACTTTGTCTCCAATAGGAAGCCGGAATTTCCCTTCGGATTCCCGTTAACCTAAAGCCTGTAACCCCTTTTCCACCGCTTTATGTAACAACTGGCGGTCATGCCGATACGGAATATCGCTGGCCTCCAGTACCCCCTGGATCACAACCCGCTCTTCCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAgct                                                                                                                                                      <  1:23329/251‑1 (MQ=255)
                                                                                      aaaCATCGCGCACTTTGTCTTCAATAGGAAGCCGGATTTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTCCCAGCCCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTCCCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAgctt                                                                                                                                                     <  1:363121/250‑1 (MQ=255)
                                                                                       aacaacGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGATGGCCTCCAGTACCTCCTGGTTCCCAATCCGCTCTTTCCCCGCCGGGACATCCACTTTTGGCCCGACGATGACCGCATCAATTCCCTTTTTACCCACCTACTGCTCCATAATTGCCAgcttgc                                                                                                                                                   >  2:222751/1‑251 (MQ=255)
                                                                                          aacGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCAACAATGACCTTTTTACCAACATCCTGCTCCATAATTGCCAGCTTGCttt                                                                                                                                                >  1:461370/1‑251 (MQ=255)
                                                                                           ccgcgCACTTTGTCTTCAATAGGAAGCCGAAATTTTCCTTCGGATTTCCTTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTCCCTCCTGGATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTc                                                                                                                                               <  1:438620/250‑1 (MQ=255)
                                                                                           acGCGCACTTAGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTTGCCCGACGATGACCGCATCACTGACCTTTTTACCCACATACTGCTCCATAATTGCCAGCTTGCTTTc                                                                                                                                               >  2:519305/1‑251 (MQ=255)
                                                                                              cgctttttGTCTTCAATGGGAGGCCGGAATTTTCCTTGGGTGTTCCGTGACCCAAAAGCCTTTAACGCTTTTGCCAGCGTGTTATTTAACAACTGGCGGTCATGACGATACGGAATACCGCTGGCCTCCAGTACCCCCTGGATCAAAACCCGCCCTTCCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAg                                                                                                                                            <  2:572595/246‑1 (MQ=255)
                                                                                               cccactttgatttaaaTAGGAAGCAGGATTTTTCCTGCGGATTTCCTTTAACCAAAAGCCTTAAACGCCTTTTCCAGCTCGTTTGGAACCAGCGCGCGTTAAGGGCGATACGAACGATCGCTGGCCCCCAGTCCCCCTTGGATCAAAATCCGTTCTTGCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATGCTGCTCCATAATTGCCAGCTTGCTTTCCAg                                                                                                                                            <  2:550242/237‑1 (MQ=255)
                                                                                                cccTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTCCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTCCCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCt                                                                                                                                          <  1:579970/249‑1 (MQ=255)
                                                                                                cACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGAc                                                                                                                                                                                                                             >  1:10887/1‑168 (MQ=255)
                                                                                                cACTTTGTCATCAATAGGAATCCGGAATTTTCCTTCGGATTTCCTTTAACCAAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGTTCATGACGATACGTAATATCGCTGGCCTCCAGTCCCCCCTGGATCACACCCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCt                                                                                                                                          <  1:520379/251‑1 (MQ=255)
                                                                                                 aCTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAACGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACAt                                                                                                                                                                           >  2:36048/1‑217 (MQ=255)
                                                                                                 aCTTTGTCTTAAATAGGAAGCCGGATTTTTCCTTCGGATTTCCGTTACCCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTCCCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACAt                                                                                                                                                                           <  1:36048/217‑1 (MQ=255)
                                                                                                           tCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCCTGACGATACGGAATATCGCTGGCCTCCAGAAACTCCTGGCTCACAAACCGCTCATTCACCGCAGAGCCATCCACTTTTTGACCGACGATGACCCCACCAATTCACTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCCGCTCCAAAtttt                                                                                                                                 >  2:337793/1‑248 (MQ=255)
                                                                                                            aaaTAGGAGGCCGAAATTTCCCTTCGGTTTCCCGTTACCCAAAACCCTTTACCGCCTTTTCCAGCGCGTTATGTAACAACTGGGGGTCATGACGATACGGAATATCGCTGGCCCCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGcc                                                                                                                               <  1:533136/250‑1 (MQ=255)
                                                                                                              aTAGGAAGCCGGCATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCACCGCGTTCTGTAACAACTGGCGGTCATAACGATAAGGAATATCGCTGGCCTCCAGCACCAACTGGCTCACAATCCGCTCTTTCACCCCCGGGACATCCCCTTGTGGCCCGACGATAGCCGCATCAATGACCTTTTTCACAACATCCTGCTCCATAATGCCCATCTTGCTTTCCCGCTTCAAATTAGcagc                                                                                                                             >  1:550535/1‑248 (MQ=255)
                                                                                                                   aaGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGt                                                                                                                         >  1:23083/1‑250 (MQ=255)
                                                                                                                    aGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGCTGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTa                                                                                                                        >  1:373510/1‑250 (MQ=255)
                                                                                                                    aGCCGAAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTACCGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCCCCAGTACCTCCTGGATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTaa                                                                                                                       <  1:574082/251‑1 (MQ=255)
                                                                                                                        gaaaTTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACCCCTTTTCCAGCGCGTTATGTACCAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTCCCCCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAaactc                                                                                                                   <  1:357950/249‑1 (MQ=255)
                                                                                                                           aTTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGAACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATAATGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGGAcactcaac                                                                                                                >  2:359495/1‑251 (MQ=255)
                                                                                                                              ttCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTACAGCTTCAAATTAGCCGACGGTAaactcaacgc                                                                                                              >  1:128347/1‑248 (MQ=255)
                                                                                                                                    cGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCCCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCCGGAAAACTCAACTCAcccccc                                                                                                       >  2:554565/1‑251 (MQ=255)
                                                                                                                                       aTTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGCCGATGACCGCATCACTGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGcccc                                                                                                       >  1:127071/1‑248 (MQ=255)
                                                                                                                                        tttCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGAGCGTTATGTAACAACTGGCAGTCATGACGATACGGAATATCGCTGTCCTCCAGTAGATCCTGGATCACAATCCGCTCTTTAACCGCCGGGACAACCAATTTTGGCCCGACGGTTACCGCAAAAATGACCTTTTTACCGACAGACTGCACAAAAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAAGTCACGCCACAgct                                                                                                   >  1:540498/1‑249 (MQ=255)
                                                                                                                                          ccctttAACCTACAGCCGGTACCGCCTTTCCCAGCGCTTTATTTAACAACTGGCGGTCATGACTATACGGAATATCGCTGGCCCCCAGTACCCCCTGGACCCCAACCCCCCTTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGAtt                                                                                                  <  1:159434/246‑1 (MQ=255)
                                                                                                                                           ccGTTACCCTAAAGCCTGTAACGCTTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCAGCAGGTAAACTCAACTCACGCCCCAGATTGc                                                                                                <  1:178909/251‑1 (MQ=255)
                                                                                                                                             gTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCg                                                                                              >  1:73048/1‑251 (MQ=255)
                                                                                                                                                aaCCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGACTATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGATTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGat                                                                                            >  1:71111/1‑250 (MQ=255)
                                                                                                                                                            gTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTNCCGCCCCAGATTGCCGATATAAACCATCGgc                                                                               >  2:287146/1‑251 (MQ=255)
                                                                                                                                                             gtaaccccttgcccAGCGCGTTATTAAACACCGGGCGGTCATACGGATACGGAAAACCGCTGCCCCCCATTCCCTCTGGGATCACATTCCCCTTTTACACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCGTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTGCAAATTACCCGCAGGTAAACGCAACTCACGCCCCAGATTGCCGATATAAACCATCgg                                                                                <  1:481761/238‑1 (MQ=255)
                                                                                                                                                              aaCGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCAGCAGGTAAACTCAACTCACGCCCCAGCTTGCCGATATAAACCATCGgcgc                                                                             >  1:381685/1‑251 (MQ=255)
                                                                                                                                                               aCGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATc                                                                                  >  1:243623/1‑245 (MQ=255)
                                                                                                                                                                actcTTTTCCCAGCGGGTTATTTAACAACGGGGTGTCATGACGAGACGGAATCTCGCTGGCCCCCATGCCCTCCTGGATCACAACCCGCCTTTTCACCTCCGAGACAGCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTGTTACCAACATACTTCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATc                                                                                  <  2:243623/242‑1 (MQ=255)
                                                                                                                                                                 gCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGcc                                                                                                                                                                                                               >  1:49127/1‑118 (MQ=255)
                                                                                                                                                                      tttCCAGCGCGTTATGTAACAACTGGCGTTAATGACGATACGGAATATCGCTGGCCTCCAGTACCCCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGt                                                                       <  1:173562/250‑1 (MQ=255)
                                                                                                                                                                                    aTGTAAAAACGGGCGGTCATGACGATACGGAATATCCCTGGCCTCCAGTACCCCCTGGATCCCAACCCGCTCTTCCCCCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCt                                                         <  1:207223/251‑1 (MQ=255)
                                                                                                                                                                                                         gACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTcagcagc                                     >  1:550427/1‑251 (MQ=255)
                                                                                                                                                                                                            gATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCg                                                                                              >  1:451604/1‑191 (MQ=255)
                                                                                                                                                                                                            gATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTcagcagc                                     >  2:92126/1‑248 (MQ=255)
                                                                                                                                                                                                            ccaTACGTAACATCCCTGGCCCCCACTACCTCCGGGACCACAACCCGTCCTTTCCCCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTcagcagcag                                   <  1:20166/249‑1 (MQ=255)
                                                                                                                                                                                                                  ggAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTgg                             >  1:293033/1‑251 (MQ=255)
                                                                                                                                                                                                                       atCGCTGGCCTCCCGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAAt                                                             >  1:399448/1‑214 (MQ=255)
                                                                                                                                                                                                                             ggCCTCCAGTACGTCCTGGATCACAAACCGCTCTTTCACAGCCGAGACATCCACTTTTGGACCGACGATGACCGCATCAATGACCTTTTTACCAACATCCTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGCAAACTCAACTCAAGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGCCGAAATGCCTGGTCGATTTCCTTCAGCAGCAGCATTGGCATCCGGcggg                  >  2:63831/1‑248 (MQ=255)
                                                                                                                                                                                                                               ccTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGt                 >  2:505659/1‑250 (MQ=255)
                                                                                                                                                                                                                                 cccAGTACCCCCTGGATCACAACCCGCTTCTTCACCGCCGAGACATCCACTTTTGGCCCTACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGtata              <  1:45282/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cctcctGGATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGTATAAAAACTg       <  1:296340/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGTATAAAAACTGCCAgg  >  1:238161/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
CGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAGTACCTCCTGGATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGTATAAAAACTGCCAGG  >  NC_000913/815598‑816080

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