Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913803,8690GT35.7% 46.9 / 16.4 42V123L (GTG→TTG) ybhJputative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base G (9/18);  new base T (15/0);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTG  >  NC_000913/803623‑804111
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGcc                                                                                                                                                                                                                                                      >  1:659839/1‑251 (MQ=255)
  tctctcTTCGCATAAGACGTCAGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACTATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCGGCCGTATGTGCTGACCAACTTCCATAACTCACTCTGCGCCGTCGGCGGGACTATTAACGGTGATGACCATGTGTTTGGTTTATCGGCGGCCAAGCTTTATGGCGGTTTTTTTGTGCCt                                                                                                                                                                                                                                                     <  2:633923/250‑1 (MQ=255)
         tCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGGCTCCGCatat                                                                                                                                                                                                                                            >  1:757829/1‑251 (MQ=255)
                        ggAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGGCTCCGCCTATTTCGGGgatccatc                                                                                                                                                                                                                              >  1:532365/1‑250 (MQ=255)
                                               tAAGTTTGATTCATTAGCCTCCCACGATATTACCTTTGTTGGTATTGGACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCGCTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCa                                                                                                                                                                                                       <  2:753/250‑1 (MQ=255)
                                                 aGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTGCGGGCATCCATCAATATATGTGTGAGATGATGGCAggg                                                                                                                                                                                                    >  1:464852/1‑250 (MQ=255)
                                                     tGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTTGCTCCCCCCATTGTGGGCATCCAACAATATATGAGTggggtgttggggggggggg                                                                                                                                                                                                >  2:705310/1‑233 (MQ=255)
                                                          cATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGTTTTTTGTGCCCTCGCATATTTGGGGCATCCATCAATATAAGCGTgaggggatgggaaggggcgggga                                                                                                                                                                                            >  1:467143/1‑230 (MQ=255)
                                                          cATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCCCACCTTTGGGTCAACCCTCAAAATATGCGTgggagtgttgcggggggggggga                                                                                                                                                                                            >  2:322315/1‑228 (MQ=255)
                                                           aTTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCCCATATTGGGGGTATTCATCAATATATGGGTGAGATGATGGGAgggggggggaaa                                                                                                                                                                                           >  1:460087/1‑240 (MQ=255)
                                                             tAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGGCCCCGCCTATTGTGGCGCTCCCTCAATATATGCGTGGGATGATAGTAGGCGGCGGCAAAATg                                                                                                                                                                                       >  1:75638/1‑251 (MQ=255)
                                                                         atatTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGCGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCaaa                                                                                                                                                                                          <  1:28015/236‑1 (MQ=255)
                                                                          tatTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCCCCGCCTATTTCGGTACTACATCAATATTTGGGTGTGATGATGGGAGGGGGCGGCAAAAAGATTCTCggggca                                                                                                                                                                           >  1:489403/1‑247 (MQ=255)
                                                                              aCCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTTCCCATATTTCGGGCATCCATCAATATATTCGTGGGGTGATGGCCGGGGGCGGCAAAATGATACTccgggcaaac                                                                                                                                                                        >  1:741150/1‑240 (MQ=255)
                                                                                   tgttgGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAGCGTTTCCCGCGGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCcac                                                                                                                                                                  <  2:741150/250‑1 (MQ=255)
                                                                                                   aCCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTCCGCCTATTTCGGGCATACATCAAAATATGCGTGGGATGATGGCAGGGGGCGGGAAAAAGATCCTCCGGGCAGACAGCCACACCCGTTAAGGGGCAtt                                                                                                                                                 >  1:218121/1‑251 (MQ=255)
                                                                                                                 cGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAg                                                                                                                                  <  1:102596/251‑1 (MQ=255)
                                                                                                                         aaCGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAg                                                                                                                          <  1:705310/251‑1 (MQ=255)
                                                                                                                              ttCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGggcg                                                                                                                      <  1:783511/250‑1 (MQ=255)
                                                                                                                                cccGCTGCCTGATGTGCTGACCACCTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATGGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGgcggcg                                                                                                                   <  2:375087/251‑1 (MQ=255)
                                                                                                                                        cgaaggtgcagcccAACTGCCATTACTCCCCCTGCGCCGTCGGCGGCACTTTAAACGGTGCTGACCGTGTTTTTGTTTTATCGGCGGCCCAGAGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATTTATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACTGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTg                                                                                                           <  2:250013/239‑1 (MQ=255)
                                                                                                                                           aTGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTa                                                                                                        <  1:489245/251‑1 (MQ=255)
                                                                                                                                               gCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGACCAGCGTTATGGCGGATGTTTTTTGCCTCCCCAAATTGCGGTCATCCATCAAAATATACGTGAGATGATGTCAGGCGGGGGCAAAATTATCCTAGGGTCAGACACCCACACCCGTTACCGTGCATTAAGGACAATGGAAGAGGGTGAGGGGGGCGGGGAGGTGGGaaaa                                                                                                     >  1:91090/1‑250 (MQ=255)
                                                                                                                                                  gACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTTCCTCCGCATATTTCGGGCATCCACCCATATATTCGTGAGATGATTGCCGGCGGCGGCCAAATGATCCTCCGGTCAGACAGCCACACCCGGTACGGTGCCTTAGGGGCAATGGCAGGCGGTGGGGGCGGCGGGGGGGTGGTAAAACAgc                                                                                                 >  1:800068/1‑251 (MQ=255)
                                                                                                                                                            cATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGTCTCCCCATAATTTGGTAATCCATCCAAATATGCGTGGGATTATGGTAGGGGGGGGCCAAAATATTCTCGGGTCAGACAGGCACACCCGGTTCGGGGAATTAAGGGAAATGGGCGGCGGGGAGGGCGGCGGTGAGtttgatcgggagagcaccagttt                                                                                       >  1:330426/1‑230 (MQ=255)
                                                                                                                                                            cATAACCCACTCTGCGCCGGCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTAGCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGtt                                                                                                            <  2:330426/230‑1 (MQ=255)
                                                                                                                                                                  ccactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTGATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACGCCCGTTAGGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTg                                                                                <  2:272733/250‑1 (MQ=255)
                                                                                                                                                                          cgcCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTTCCCCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCg                                                                        <  2:460087/251‑1 (MQ=255)
                                                                                                                                                                                    ggcACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTGTGCCCCCCCATATTGCGGGCATCCATCAATATATTCGTGGGGTGATGGCAGGGGGCGGGAAAATGATCCTCGGGGCAGACAAGCAAACCCGGTACGGGGCATTAGGGACAAATGCAGTAGGGGAGGGGGGCGGCGAGGTGGTAAAAAAGATGGTTAATGGCAACTGGGGTATCGACTATCCggg                                                            >  1:506230/1‑251 (MQ=255)
                                                                                                                                                                                              aCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGGGCCTCCGCATATTTCGGGCCTCCCTCAATATATGCCTGGGATGATGGCAGGCGGCGGGAAAATGATCCTCGGGGCAGACAGCCACCCCCGTTACGGGGCATTAGGGACAAAGGGAGGCGGGGAGGGCGGGGGGGAGGTGGGAAAAAAGCCGCTTAATGAAACCCGGGGATTTGGCTATCCGGGCGTGGTTggg                                                  >  1:327652/1‑249 (MQ=255)
                                                                                                                                                                                              aCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGAATTTTTTTGGCTCCGCACTTTTGGGTCATCCCTCCAAATATGAGTGAGGAGATGGGAGGGGGCGGGAAAATTATCCCCGGGGCAAAAAGCCACACCCCTTACCGGGGATTAGGGGAAATGGTAGACGGTGGGGGGGGGGGGGAGTGGGTAAAACAAACGGTTAATTACAACTGGGATATTGACTATACCGGGGTGGGTGCg                                                  >  1:839530/1‑251 (MQ=255)
                                                                                                                                                                                                 gTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGGCTCCGCATATTTCGGTCATCCATCAATATATGCGGGAGGAGATGGCAGGCGGCGGCAAAATGATTCTCGGGTCAGACAGCCACACCCGTTACGGTGCCTTAGGGAAAAAGGCAGTCGGGGGGGGCGGCGGGGAGGTGGTAAAAAAACTGCTTAATGACACCCGGGATATCGACTATCCCGGCGTGGTTGCgggg                                               >  1:17931/1‑249 (MQ=255)
                                                                                                                                                                                                                       aTCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCCTATTTCGGGCATCCCTCAATATATGGGGGGGATGATGGGAGGCGGCGGCAAAATGATTCTCCGGGCAGACAGCCACACCCGGTACGGGGGATTAAGGGCAATGGGAGGCGGGGAGGGGGGGGGGGGGGTGGGAAAAAAGCTGCTTAATGAAACCCGGGGTATTGACCTTTCCGGGGGGGGTGCGGGGCATCTGACCGGGAAACCCgcgc                         >  1:20441/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGc             <  2:582378/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ttATGGCGGTATTTTTGTGCCTCCGCCTAATTCGGGCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGGCAGAAAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGGGAGGGCGGCGGGGGGGTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGGCTATCCGGGCGGGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCCGATGTGGGGCcac          >  1:552141/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aTGGCGGTATTTTTTTGCCTCCGCATATTGCGGTCCTCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACAccc                                                                                                                                                              >  2:163917/1‑105 (MQ=255)
                                                                                                                                                                                                                                        aTGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACAccc                                                                                                                                                              <  1:163917/105‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGTGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGa      <  2:470325/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGAtg    <  1:571166/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCAGCGCCGTATGTGGGGCCACAGGCtgtg  >  2:110312/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     ttGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTg                            >  1:577479/1‑219 (MQ=255)
                                                                                                                                                                                                                                                     ttGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTg                            <  2:577479/219‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
ATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTG  >  NC_000913/803623‑804111

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