Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913803,8690GT26.7% 49.8 / 18.4 46V123L (GTG→TTG) ybhJputative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base G (3/30);  new base T (12/0);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 1.58e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCT  >  NC_000913/803628‑804116
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                            
gtctTCGCATAATACGTCCGGACATATGGTTAAACTTAAACTTATGTTTGATTCATTAGCATCTCACTATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCAGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGTCGGCACTATAACCGGTTCTGACCATGTTTTTGGTTTAGCGGCGTCACAGTGTTCTGTCTGGTATTTTTATGCCttct                                                                                                                                                                                                                                                      >  2:290490/2‑248 (MQ=255)
ctctTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGGATTGTACAGACCGCTAAAGCGTCCGGTATCGAACGGTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACGCTGCGCCGTCGGCGGCACTATGAACGGTGATGACCGTGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTGTGTGCCTcc                                                                                                                                                                                                                                                       <  2:182606/249‑1 (MQ=255)
   ttCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTccccct                                                                                                                                                                                                                                                   >  1:114365/1‑246 (MQ=255)
     cGCATAATACGTCCGGAAATATGGATAAACTGAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCatat                                                                                                                                                                                                                                                <  2:483220/250‑1 (MQ=255)
       cATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTTATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATGGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGc                                                                                                                                                                                                                                             <  2:127093/251‑1 (MQ=255)
                tCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGTTTTTTTTTCCTCCCCATATTTCgggcctccc                                                                                                                                                                                                                                    >  1:514890/1‑244 (MQ=255)
                 ccGGAAATATGGATAAACTTAAAATTGAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCACAGCGTTTTGGCGGAAGTTTTTTGCACCCGCATATTTCGGGcatgcc                                                                                                                                                                                                                                   >  1:211653/1‑247 (MQ=255)
                      aaTATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCa                                                                                                                                                                                                                                <  1:13619/248‑1 (MQ=255)
                           ggATAAACTTAAAATTAAGTTTGATTCTTTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGc                                                                                                                                                                                                                        <  2:137736/251‑1 (MQ=255)
                            gATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCGCACATTTCGGTcatacatatatagga                                                                                                                                                                                                                          >  1:417286/1‑240 (MQ=255)
                            gATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTGGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTcatccat                                                                                                                                                                                                                                  <  2:417286/240‑1 (MQ=255)
                                        aTTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGTCTCCCCGTACTTAGGTCATCCCTTAAAATATGTGTGAGATGATggg                                                                                                                                                                                                           >  2:394097/1‑250 (MQ=255)
                                               ttGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCTGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAggcgg                                                                                                                                                                                                     <  1:351112/250‑1 (MQ=255)
                                               ttGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAggcgg                                                                                                                                                                                                     <  2:114365/250‑1 (MQ=255)
                                                               ccACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCg                                                                                                                                                                                  <  1:524177/250‑1 (MQ=255)
                                                                   gATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTTCGGTCATCCCTCAAAATATGCGGGGGGTGATTGGAGGCGGGGGCAAAAAGATCCTCGGGtt                                                                                                                                                                              >  1:439920/1‑250 (MQ=255)
                                                                      atTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGGCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCGGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATGGTCGGCCCATCTTTATGGCGGTATTTTTGGGCCTCCGCATAGTTCGTTCATCCATGAAGATATGCGGGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAg                                                                                                                                                                            <  1:39772/250‑1 (MQ=255)
                                                                         aCCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGAATGGAACGTTTCCCGATGCCGTATGTGATGACCAACTGCCATAACTCACTCTGAGCCGTCGGCGGCACTATTACCGGTGATGACAATGTTTTTGGTTTATCGGCGGCCAAGCGTTATGGCGGTATTTTTGTGCACCCGCAAATTGAGGTAAATCCTTAATAATTGCGTGAGAAGGTTGCAGGCGGCGGACAAAAtatcccacagacagac                                                                                                                                                                           >  1:68063/1‑234 (MQ=255)
                                                                            tttGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCCTTTTTCGGGCATCCCTCAATATATGCGTGGGATGATGGCAGGCGGCGGCAAAATGATCCCCCGGTCAGACAGccc                                                                                                                                                                     >  1:428997/1‑250 (MQ=255)
                                                                               ggtgGTATTTTCCAGTCCGCAAAAGGGTGCGGTATGAAGCGGTGCCCGCTGCCGTATGTGCGGACCATCTGCCATAACTCCCTGTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTAGTTTGGTGCCTCCGCAGATTGCGGTAATCCATCAATATATGAGTGAGATGATGGCAGGCGGCGGCAAAAGGATCCTCGGGTCAGACAggaaca                                                                                                                                                                   <  2:385336/248‑6 (MQ=255)
                                                                                     aTTGTACAGCCCGATAAAGCGTCCGGAATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGTCTCCCCCTTATTCGGTCCTCCATCAATATATTCGTGAGATGGTGGCAGGGGGCGGCCACAATATCCTCGGGTCAGACACCCCCACCCgctt                                                                                                                                                            >  1:440880/1‑248 (MQ=255)
                                                                                        gTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCCCATATTTCGGTCATCCCTCAATAATTGGGTGGGAGGATGGCAGGGGGGGGCAAAATGATCCTACGGTCAGACAGCCACACCCGTTACgg                                                                                                                                                         >  1:304458/1‑251 (MQ=255)
                                                                                          aCAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGt                                                                                                                                                        <  1:67243/250‑1 (MQ=255)
                                                                                                    aaaGCGTCCGGTATGGAACGTTTCCCGCGGCCGTATGTGCTGACCACCTGCCATAACTCACGCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGAc                                                                                                                                             <  2:439920/251‑1 (MQ=255)
                                                                                                            cGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTTTTTTTTTGGCTCCGCAAATTTAGGGCCTCCATCCATATATGCGTGGGATGATGGGAGGGGGGGGCAAAAAGATCCTCCGGGCAGACAGCCACCCCCCGTACGGGGCATTTGGGACAAAGGccg                                                                                                                                     >  1:472102/1‑249 (MQ=255)
                                                                                                                                    gTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTTTTTTTGTGCCTCCGCATATTGCGGGCCTCCCTCCATATATTCGTGAGATGATGGCAGGGGGCGGCAAAAAGATACTCGGGGCAGACAGACACACCCGGTAAGGGGCATTAGGGACAATGGCAGGCGGTGAGGGcggcggggggtggg                                                                                                             >  1:63980/1‑243 (MQ=255)
                                                                                                                                             gACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAg                                                                                                     <  1:394097/250‑1 (MQ=255)
                                                                                                                                                        aTAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGac                                                                                         <  1:460543/251‑1 (MQ=255)
                                                                                                                                                                      gcCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCCCCCCCTATTTCGGTCATCCCTCAATATATGCGTGGGGTGGTGGCAGGGGGGGGCAAAATGATCCCCGGGTCAGACAGCCCCACCCGTTACGGGGGATTAGGGGCAAAGGCAGTCGGGGGGGGGGGCGGGGGGGTTGTAAAACAACGGCTTTATGACACCTGGGGTATCgc                                                                           >  1:261796/1‑250 (MQ=255)
                                                                                                                                                                       ccGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGAc                                                                          <  1:320422/251‑1 (MQ=255)
                                                                                                                                                                         gTCGGCGGCACTATTATCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTa                                                                        <  1:410962/251‑1 (MQ=255)
                                                                                                                                                                                gcACTATTAACGGTGAGGACCATGTTTTTGGTTTAGCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGc                                                                <  2:434889/251‑1 (MQ=255)
                                                                                                                                                                                    tATTAACGGCGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTCCCCCTATTTCGGGCATCCATCAATATATGCGTGAGATGATGGGAGGGGGCGGCAAAAAGATTCTCGGGGCAGACAGCCACACCCGTTACGGTGGATTAGGGACAATGGCCGGCGGGGGGGGGGGCGGGGGGGTGGGAAAAAAACTGCTTAATAACAACCGGGGTATCGACTATTCCggggggg                                                             >  1:134695/1‑247 (MQ=255)
                                                                                                                                                                                       tAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTgatt                                                          <  1:11749/250‑4 (MQ=255)
                                                                                                                                                                                           ggTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGt                                                     <  2:211719/251‑1 (MQ=255)
                                                                                                                                                                                               atgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATAGGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTGAATGACACCTGGGATATCGACTATCCGGGCGGGGTTGCGGTGCAt                                                 <  2:440880/251‑1 (MQ=255)
                                                                                                                                                                                                   ccATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTaa                                                                                             <  1:267308/204‑1 (MQ=255)
                                                                                                                                                                                                     aggTTTTTGGTTGATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGAcc                                           <  2:514890/249‑1 (MQ=255)
                                                                                                                                                                                                              gTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTGCGGTCATCCCTCAATATATGCGGGAGGTGATGGCAGGCGGCGGCAAAATGATCCTCGGGGCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGGCGGTGGGGGCGGCGGTGGGGTGGGAAAACAGCTGCTTAATGACACCTGGGATATCCACTATCCCGGCGGGGGTGCGGGGCCTCTGACCGGaaaaacc                                  >  1:72460/1‑248 (MQ=255)
                                                                                                                                                                                                                ttATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGTATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTg                                 <  1:502364/250‑1 (MQ=255)
                                                                                                                                                                                                                 tATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTgc                                <  1:118214/250‑1 (MQ=255)
                                                                                                                                                                                                                    cggcggCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGcc                             <  1:164172/250‑1 (MQ=255)
                                                                                                                                                                                                                         gcCCAGCGTTATGGCGGTATGTTTGTGCCTCCGCATATTGCGGTCATACATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGTATATCGACTATCCGGGCGTGGTTGCGGGGCATCTGACCTTAAAACCTGCGCCGTAtgt                       <  2:472102/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gaTTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTgg      <  2:134695/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aTTTTTGTGCCTCCGCGTATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAAGGATCCTCGGGTGAGGCAGGCACTCCCGTTACGGTGCAGTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACTCCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTTACCGGAAAACCTGCGTCGTATGTGGGGCCACAGGATGTGgc     <  1:387691/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ttGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGGTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTTCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCt  <  1:470901/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                            
CTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCT  >  NC_000913/803628‑804116

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