Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913803,8690GT32.5% 60.3 / 12.9 40V123L (GTG→TTG) ybhJputative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base G (4/23);  new base T (12/1);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 3.68e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTAT‑CCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCTGGC  >  NC_000913/803623‑804119
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                              
aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGcc                                                                                                                                                                                                                                                            <  1:438179/251‑1 (MQ=255)
   cactctACGCATAATACGAGCAGAAATATGAATAAACATAAAGTTAAGTTTGATTAATTAGCCTCTCAAGATATTTGCTTTTTTGGTATTGTACAGACAGGTAAGGCTGCCGGTATTGACCGTGTCCCGGTGCCGGATGTGCTGACCAACTGCGATAACTCTCTCTGCGCCGTCGGCGGTACTGTTAGCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCTGTATTTTTGTGCCTcc                                                                                                                                                                                                                                                         <  2:418694/249‑1 (MQ=255)
     ctctTCTCATAATTCGTCCGGAAATATGGATAAACTTAAAATTAAGTTTTATTCATTAGCCTCTCACGATATTACCTTTGTTGTTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTc                                                                                                                                                                                                                                                          <  2:103924/248‑1 (MQ=255)
     cactTCGCATAATACGTCCGGAAATATGTATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACTATATTTCCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCTCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTTATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGTCGGTATTTTTGTGCCTcc                                                                                                                                                                                                                                                         <  1:29086/247‑1 (MQ=255)
            cATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGc                                                                                                                                                                                                                                                <  1:139792/251‑1 (MQ=255)
             ataataCGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTCCGCATAtttcg                                                                                                                                                                                                                                               >  1:454300/1‑248 (MQ=255)
                              atGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAatatat                                                                                                                                                                                                                              <  2:501069/251‑1 (MQ=255)
                                                 aGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCCCCTTTTTCGGTCCTCCAACAATATATGTG‑TGTGA‑TGGTGGGAggg                                                                                                                                                                                                         >  1:77584/1‑250 (MQ=255)
                                                 aGTTTGATTCATTAGCCTCTCACGATAGTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAggc                                                                                                                                                                                                         <  1:277308/251‑1 (MQ=255)
                                                  gTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTCCGCCTATTGGGGTCATCCAACAatatatttgtggggtggtgggagggg                                                                                                                                                                                                          >  1:457929/1‑231 (MQ=255)
                                                   tttGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCCCCTATTGCGGGCATCCCTCAATATATTCG‑TGGGA‑TGATGGCAggcgg                                                                                                                                                                                                       >  1:239345/1‑251 (MQ=255)
                                                           aTTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCCTATTGCGGGCATCCCTCAATATTAGCG‑GGAGA‑TGATTGCAggcggcgggcca                                                                                                                                                                                                >  1:392111/1‑246 (MQ=255)
                                                             tAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAAt                                                                                                                                                                                              <  1:258502/250‑1 (MQ=255)
                                                              aGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGa                                                                                                                                                                                            <  1:418812/251‑1 (MQ=255)
                                                                 ctctCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCGCTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTTTCTTCCCCTATTTCGGGCCTCCATCAAAATATGCG‑GGGGA‑TGGAGGCAGGCGGGGGCaaaaagaatc                                                                                                                                                                                         >  1:135771/1‑245 (MQ=255)
                                                                         atatTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCCTTTTTCGGGCATCCATCCATCTATGCGGTG‑GAGTG‑TGGCCGGCGGCGGCCAAATGATCCTCgggggc                                                                                                                                                                                 >  1:17654/1‑248 (MQ=255)
                                                                                  tggttgGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCca                                                                                                                                                                         <  2:239345/248‑1 (MQ=255)
                                                                                       ggTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTTCCCATATTTCGGTCATTCCTTAATATTTTCG‑GGGGA‑TGATTGGAGGCGGGGGCAAAAAGGTCCTCGGGGTAAAAAGCCAC‑Acccc                                                                                                                                                                  >  1:139458/1‑250 (MQ=255)
                                                                                        gTATTGTACACACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGt                                                                                                                                                                 <  1:138511/251‑1 (MQ=255)
                                                                                          aTTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGTTGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGtt                                                                                                                                                                <  1:291428/250‑1 (MQ=255)
                                                                                               aCAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTTCCTCCGCATATTTCGGGCATCCATCCATATTTTCG‑GGAGA‑TGATGGCAGGCGGGGGGAAAATGATCCTCGGGTCCGACAGCCCA‑ACCCCTTACgggg                                                                                                                                                          >  1:135727/1‑249 (MQ=255)
                                                                                                         aagtcttCCGTTATGGAACGTTTCCCGCTGCCGTATGGGCTGACCATCTTCCATAACTCACTCTTCGCCGTCGGCTGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTGCGGTCATCCAGCAATATATGCG‑TGAGA‑TGATGGCAGTCGGCGGCAAAATGATCCTCGGATCAGACAGCCCC‑ACCCTTTACGGTGCATTAggg                                                                                                                                                  <  2:279938/243‑1 (MQ=255)
                                                                                                                                     tGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTTCGGGCATCCATCCATATATGCG‑TGAGA‑TGATTGCAGGCGGCGGCAAAATGATCCTCCGGGCAGACAGCCAC‑ACCCGTTACGGGGCATTTGGGACAATGGCAGGCGGGGAGGGcggcgggg                                                                                                                      >  1:333817/1‑247 (MQ=255)
                                                                                                                                             gtgcgtaCCAACTGCCATAACTCACGCTGCGCCGGCGGCGGCACTAGTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCAGCAAAATGATCCTCGGGTCAGACAGCCAC‑TCCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTaa                                                                                                             <  2:333971/244‑1 (MQ=255)
                                                                                                                                                     cAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCCCCTATTTCGGGCATCCCTCAAAATATGCG‑TGGGA‑TGAATGCAGGGGGCGGGAAAATGATACCCGGGGCCGACCGCCAC‑ACCCCGTACGGGGCAATAGGGGCAATGGCAGTCGGGGAGGGGGGGGGGGGGGTGGGAAAACAgctgg                                                                                                    >  1:349703/1‑250 (MQ=255)
                                                                                                                                                           ccATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACATCTGCTTAATg                                                                                              <  1:54242/251‑1 (MQ=255)
                                                                                                                                                            cATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAAtta                                                                                             <  1:86834/251‑3 (MQ=255)
                                                                                                                                                                   cactcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCCTATTTCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAAAGATCCTCGGGTCAGAAAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGGGAGGGGGGCGGTGAGGTGGGAAAACAGCTGCTTAATGACACCTgg                                                                                      >  1:206680/1‑251 (MQ=255)
                                                                                                                                                                       ctGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTGTGCCTCCGCATATTTCGGGCATTCATCAATATATGCG‑GGAGA‑TGATGGGAGGGGGCGGGAAAATGAACCTCGGGGGAGAAAGCCAC‑ACCCGGTACGGGGGATTTGGGGCAATGGGAGGCGGTGAGGGCGGGGGGGGGGTGGGAAAAAAGATGCTTAATGAAAACTggggg                                                                                   >  1:236941/1‑248 (MQ=255)
                                                                                                                                                                            ccGTCGTCGGCACTATTAACTGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCTTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTTGGATATCGAc                                                                             <  2:139458/251‑1 (MQ=255)
                                                                                                                                                                                                   gatgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTAT‑CCGGGCGTGGTTGCGGTGc                                                      <  1:211179/250‑1 (MQ=255)
                                                                                                                                                                                                       tccATGTTTTTGGTGTAGCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTAT‑CCGGGCGTGGTTGCGGTGCATCt                                                  <  2:352472/249‑1 (MQ=255)
                                                                                                                                                                                                            gTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCCCATATTTCGGTCATCCCTCAATATATGCG‑TGAGG‑TGATGGCAGGGGGCGGCAAAAAGAACCCCCGGTCAAACAGCCAC‑ACCCGCTACGGGTCATTAGGGAAAAAGGGAGGCGGGGAGGGGGGGGGGGGGGTTGTAAAAAAGCCGCTTAATGACACCCGGGGTATTGACTATTCCggggggggtgcggggcattttgcccg                                             >  1:462733/1‑228 (MQ=255)
                                                                                                                                                                                                                ttGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTTCGGGCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGGGGCGGCAAAATGATCCTCGGGGCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGGGAGGGCGGCGGTGAGGTGGTAAAACAGCTGCTTAATGAAAACCGGGGTATCGACTAT‑CCGGGCGTGGTTGCGGGGCATCTGACCGGaaa                                         >  1:384220/1‑250 (MQ=255)
                                                                                                                                                                                                                                ccAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTAT‑CCGGGCGTGGTTGCGGTGCAGCTGACCGGAAAACCTGCGCCGTATGTgg                        <  2:315551/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGCGTTATGGCGGTATTTTTGTTCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTAT‑CCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTggg                       <  2:463841/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTTCGGCCATCCATCCATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATACCCGGGTCAGACAGCCACCACCC‑TTACGGAGCATTAGGGAAAATGGCAGGCGGCGAGGGCGGGGGCGAGTTGGGAAAAAAGCTGCCTAATGAAACCCGGGGTATCGACCAT‑TCGGGCGCGGGTGCGGTGCATCTTGCCGGAAAACCTTCCCCGTATGCggg                       >  1:394731/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTAT‑CCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCcac                  <  1:99717/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cgttATTTTTGTGCCTCCGCATAGTGCGTTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTAT‑CCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGAtg            <  2:427092/248‑1 (MQ=255)
                                                                                                                                                                                                                                             ggTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTAT‑CCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGAtgt           <  2:454300/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      gtgCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTAT‑CCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCTGGc  >  1:313746/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                              
ATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCG‑TGAGA‑TGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCAC‑ACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTAT‑CCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGCTGGC  >  NC_000913/803623‑804119

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