Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913803,8690GT23.8% 80.3 / 17.7 42V123L (GTG→TTG) ybhJputative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base G (3/29);  new base T (10/0);  total (13/29)
Fisher's exact test for biased strand distribution p-value = 1.94e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGG  >  NC_000913/803636‑804112
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
ataataCGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑GTTTTTTTTGGCTCTGCCTATTGCg                                                                                                                                                                                                                                     >  1:49260/1‑251 (MQ=255)
  aataCGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑GAATTTTTTGCCTCCGcccatttgggg                                                                                                                                                                                                                                   >  1:112989/1‑241 (MQ=255)
       gTCCGGAAATATGGAAAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTGGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTTCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCat                                                                                                                                                                                                                                           <  2:433156/238‑1 (MQ=255)
            gAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGTCTCCGCATATTTCgggcctcccttc                                                                                                                                                                                                                          >  1:179721/1‑240 (MQ=255)
              aaTATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑GATTTTTTTGGCTCCCCCTATTGCGGGCATCCCTCAat                                                                                                                                                                                                                        >  1:117732/1‑250 (MQ=255)
                 atGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAatatat                                                                                                                                                                                                                    <  2:81086/251‑1 (MQ=255)
                 atGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGAATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCGACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGCCGGCCCAGCGTTATGGCGG‑TATTTTTGTTTCTCCGCGTATTgtcgtaataccttcatatat                                                                                                                                                                                                                    >  2:354803/1‑232 (MQ=255)
                         aCTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTgaga                                                                                                                                                                                                            <  1:335339/251‑1 (MQ=255)
                                        tGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGCCTCCCCCTATTTCCGGCCTCCATCAATATATTTGTGGGATGATGG‑GAgggggcg                                                                                                                                                                                            >  1:240262/1‑251 (MQ=255)
                                             cATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGt                                                                                                                                                                                                                                   <  1:158629/208‑1 (MQ=255)
                                                   ccTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTTGCTCCGCATATTTCGGGGATCCATCAATATTTGCGTGGGGTGGTGG‑CGGGGGGCGGGAAAATGATc                                                                                                                                                                                 >  1:165276/1‑251 (MQ=255)
                                                      ctcgcGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGTTTTTTTTTGCC‑CCCCCTATTTTGGGCCTCCCTCAATATATGGGTGTGAGGATGG‑CGGGCGGCGGCaaaaaaaacacc                                                                                                                                                                              >  2:300179/5‑243 (MQ=255)
                                                        cACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGGGTTGACCAACTGCCATAACTCCCTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTc                                                                                                                                                                              <  2:497561/249‑1 (MQ=255)
                                                        cACGATATTACCTTTGTTGGAATTGTACAGACCGCTAAAGCGTCCGGTAAGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTc                                                                                                                                                                              <  2:149308/249‑1 (MQ=255)
                                                         aCGATATTTCCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCggg                                                                                                                                                                           <  2:117732/251‑1 (MQ=255)
                                                           gATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑GATTTTTTTGCCTCCCCCTATTGCGGTCATCCACCAAAATTTTCCGGAGGTGATGG‑GCGGCGGGGGCAAAATGATCCTCggggg                                                                                                                                                                         >  1:299782/1‑249 (MQ=255)
                                                              atTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGa                                                                                                                                                                      <  1:179909/251‑1 (MQ=255)
                                                                                                cGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCTCTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAGGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATg                                                                                                                                    <  1:43889/251‑1 (MQ=255)
                                                                                                cGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATg                                                                                                                                    <  1:289082/251‑1 (MQ=255)
                                                                                                 gTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATgg                                                                                                                                   <  1:300179/251‑1 (MQ=255)
                                                                                                                             tATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTgg                                                                                                        <  1:468989/250‑1 (MQ=255)
                                                                                                                              aTGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGt                                                                                                       <  1:353233/250‑1 (MQ=255)
                                                                                                                                  cctgcccAACTGCCATACCCCCCTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGGGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAAc                                                                                                  <  2:98316/246‑1 (MQ=255)
                                                                                                                                                ataaccccctcggcgcCGTCGGCTGCACTATTAACGGTGATGCCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATGAGGGGCAATGGCAGTCGGTGAGGGCGGCGGTTAGTTGGTAAAACAGCTGCTTTATGac                                                                                    <  2:447294/239‑1 (MQ=255)
                                                                                                                                                           tGCGCCGTCGGCGNCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAgctgct                                                                                           <  2:219195/233‑1 (MQ=255)
                                                                                                                                                            gcgcCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATc                                                                        <  2:106504/251‑1 (MQ=255)
                                                                                                                                                              gcCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGa                                                                      <  1:245968/251‑1 (MQ=255)
                                                                                                                                                               ccGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑GATTTTTTTGCCCCCCCAAATTGCGGTCATCCCTCAATATATTTGTGGGATGGTGG‑GAGGCGGGGGCAAAATGATCCTCGGGGCAGACAACCCCACCCGTTACGGGGCATTAGGGGAAAAGGCAGTCGGGGGGGGGGGCGGGGGGGTGGGGAAACAGCTGCTTAATGACACCTGGGGTATCGAc                                                                     >  1:324224/1‑251 (MQ=255)
                                                                                                                                                                cGTCGGCGGCACTATTAACGGTGAGGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGGGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACt                                                                    <  1:300868/251‑1 (MQ=255)
                                                                                                                                                                      cggcACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGAGCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGggc                                                                                                                     <  2:265549/196‑1 (MQ=255)
                                                                                                                                                                        gcACTATTAACTGTGATGTCCATGTTTTTGGTTAATCGGCGGCCCAGCGTTCTGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGAGGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGc                                                            <  1:481357/251‑1 (MQ=255)
                                                                                                                                                                               tAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTTTGCCTCCGCCTATTTGGGGCATCCCTCAATATATGCGGGGGATGGTGGGCGGGGGGGGGCAAAA‑GATCCTCGGGGCAGACAGCCACACCCCTTACGGGGGCTTTGGGACAAAGGGAGTCGGGGGGGGGGGGGGGGGGGTGGTAAAAAAGCTGCTTTATTAAAACCGGGGTATCGACTATTCGGGGGGGGTTg                                                     >  1:64046/1‑251 (MQ=255)
                                                                                                                                                                                       tggaCCATGTTTTTTGTTTATCGGCGGCCCAGCGTTATGTCGT‑TATTTTTGTGCCTCCGCATATTGCTGTCATCCAGCAATATGTGCGTGAGATGATGG‑CAGGCGGCTGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCAtt                                                                                                                                              <  2:109842/152‑1 (MQ=255)
                                                                                                                                                                                         gcccATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCt                                           <  1:95052/249‑1 (MQ=255)
                                                                                                                                                                                                  tttGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGaa                                   <  1:327352/250‑1 (MQ=255)
                                                                                                                                                                                                   ttGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGaaa                                  <  1:445535/250‑1 (MQ=255)
                                                                                                                                                                                                     ggTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAAc                                <  1:455363/250‑1 (MQ=255)
                                                                                                                                                                                                         tATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTgc                            <  1:194727/250‑1 (MQ=255)
                                                                                                                                                                                                                           aTGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAg         <  1:192172/251‑1 (MQ=255)
                                                                                                                                                                                                                                gg‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGAtgt    >  1:435330/1‑251 (MQ=255)
                                                                                                                                                                                                                                gg‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGAtgg    >  2:411130/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aTTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTgg  <  1:85165/250‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
ATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGG‑TATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGG‑CAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACACCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTGACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGG  >  NC_000913/803636‑804112

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