Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,578,287	0	T	C	46.7%	79.2 / 14.1	60	E147G (GAA→GGA)	gntR	d‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (26/6);  new base C (0/28);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.45e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGC  >  NC_000913/3578141‑3578436                                                                                                                                                    |                                                                                                                                                       tgatAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTccac                                                                                                                                                    >  1:480979/1‑151 (MQ=255)    tAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCcacca‑c                                                                                                                                                >  2:253273/1‑151 (MQ=255)       tACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑Cggg                                                                                                                                             >  2:264911/1‑151 (MQ=255)         cGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGaa                                                                                                                                           <  1:590833/151‑1 (MQ=255)         cGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGa                                                                                                                                            >  2:30068/1‑150 (MQ=255)          gTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAAt                                                                                                                                          <  2:774458/151‑1 (MQ=255)            tcgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATAc                                                                                                                                        >  1:930687/1‑151 (MQ=255)                  agaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGcc                                                                                                                                   >  2:9113/1‑150 (MQ=255)                    aCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATCCCCGCCa                                                                                                                                  <  1:383183/149‑1 (MQ=255)                     cGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAC‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATCCCCGCCACtt                                                                                                                               <  2:586643/151‑1 (MQ=255)                         cgcCGAGATAGGCAATGGGGCGACGCCCGCGAGAAATAAGGGCGGTGGTCTTCTGGCGTGCTCCTAAAAATTTAACAAAC‑CCGACGGCAATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCCCCA‑CGGAAATCCCCGCCACTTCaa                                                                                                                            <  2:151332/150‑1 (MQ=255)                           ccGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAatc                                                                                                                          >  2:952856/1‑150 (MQ=255)                                   ggaaaTGGGGCGAGGCCCGCGAGCAAAAAGGGCGGTGGTCATCGGGCGTGCTGCTAAAAAGTAACAAAAC‑CCGACGGCGATATAAGGGCAGGGCGTTTTGCTTCCCATCAGTCCCCCCA‑GGGAAATCCCCGCCACTTCAATCATCTTTAAg                                                                                                                 <  2:777343/148‑1 (MQ=255)                                      agtgtgGCGATCCCCGCGAAAAAAAATGGAGGTGGTCTTCGGGCTTGCTCCTTCAAA‑TTAACAAAACCCGGCGGCAATATCAAGGCAGGGCTCCTTGTTGCCCATCAGTCCCACCC‑CGGAAATCCCCGCCACTTCAATCATCTTTa                                                                                                                   <  2:1064495/144‑1 (MQ=255)                                        tgtgGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTgcg                                                                                                            >  2:731402/1‑151 (MQ=255)                                           ggCGAGCCCCGCGAGAAAAAATGGCGGGGTTCATCTGGCGTGCTGCTTAAAATTTAAAAAAC‑CCGCCGGCGATAACAAGGCAGGGCGTTTTTCTGTCCATAAGTCCCCCCA‑GGGAAATCCCCGCCACTTCAATCATCTTTAAg                                                                                                                 <  2:156241/143‑1 (MQ=255)                                            gCGATGCCCGCGAGAAAAAATGGCGTTGGTCATCTGGCGTGCTGCTTCAAATTTACCAAAC‑CCGCCGGCTATATCAGGGCATGGCTATTTGCTGTCCATCAGTCCCCCCA‑CGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCgcggc                                                                                                        <  2:457120/151‑1 (MQ=255)                                              gacgCCCGGGAGAAATAAGGGCGGGGTTCATCTGGCGTGCTTCTCAAAATTTACAAAAC‑CCGCCGGGAATATCAAGGCTGGGCGTCTTTCGGCCCATCAGTCCCGCGC‑GGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgt                                                                                                      <  2:666274/148‑1 (MQ=255)                                                 gCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGtgtg                                                                                                   >  1:787291/1‑151 (MQ=255)                                                 gCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGtgtg                                                                                                   >  1:1156277/1‑151 (MQ=255)                                                  cccGCGAAAAATAATGGCGGTGGTCCTCTGGCGTGTTTTTAAAAAGTAACAAAAC‑CCGACGGCAATAAAAGGGCTGGGCGATTTGCTGCCCATCAGTCCCCCCC‑GGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTgg                                                                                                  <  2:1040175/151‑1 (MQ=255)                                                     gcgAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGt                                                                                                >  1:41757/1‑150 (MQ=255)                                                          aaaaaggggCGGTGGTCTGCGGGCGTGCTTCTAAAAATTTAAAAAAA‑CCGCCGCGGATAACAAGGCTTGGCGATTTGCTTTCCATCACTCCCACCG‑GGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                                                                          <  2:1162705/144‑1 (MQ=255)                                                           aataatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCg                                                                                         >  1:763646/1‑151 (MQ=255)                                                           aataagGGGGGTGGTCATCGGGCGTGCTGCTCAAAATTTATAAAAC‑CCGACGGCAATATAAAGGCAGGGCGTTTTGCTGTCCATCAGTCCCCCCG‑GGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTAc                                                                                              <  1:978710/146‑1 (MQ=255)                                                              aagggCAGGTGTCATCGGGCGTGCTTCTTAAATTTTAAAAAAC‑CCGCCGGAAATATCAGGGCATGGCGATTTGCTGTCCATCAGTCCCCCCA‑CGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                                                                      <  1:843923/148‑1 (MQ=255)                                                                 gggaggggtctgcgggCGTCTTGCTTAAAATTTAAAAAAC‑CCGCCGGCAATATAAGGCCAGGGCCCCTTGTTGTCCATCAGTCCCCCCG‑GGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                                                                    <  2:288784/136‑1 (MQ=255)                                                                       ggTCTTCGGGCTTGCTCCTTAAAATTTAAAAAAC‑CCGACGGCGATATCAAGGCAGGGCGTCTGGCTGTCCAGCAGTCCCCCCA‑GGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAg                                                                              <  1:867644/150‑1 (MQ=255)                                                                         tcattgggggTGCTGCTTAAATGTTAAAAAAC‑CCGCCGGCGATACCAGGGCAGGGCGTTTTGCTCTCCTTCAGTCCCCCCA‑CGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAgg                                                                             <  1:50150/140‑1 (MQ=255)                                                                          caccgggggttctgctTAAAAGTAATAAAAC‑CCGCGGGCGATATAAGGGAATGGCGGTTTGCTGCCCTTCTGTCCCCCCA‑CGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCg                                                                          <  1:1046413/140‑1 (MQ=255)                                                                           aTCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGt                                                                         >  2:771794/1‑151 (MQ=255)                                                                           aTCTGGCGTGCT‑CTTAAAAATTAACAAAACCCGACGGCAATATAAGGGCATGGCGCTTTGCTGTCCATCAGTCCCACCG‑GGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGcc                                                                           <  2:1052517/149‑1 (MQ=255)                                                                                     ctgttaaaaaGTTAAAAAAC‑CCGACGGCGATATAAAGGCAGGGCGCCTTGCTGTCCATCAGTCCCCCCG‑GGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATAtt                                                                  <  2:584460/141‑1 (MQ=255)                                                                                       gctTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                             >  1:596879/1‑151 (MQ=255)                                                                                          tCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCCgg                                                             >  2:1142348/1‑148 (MQ=255)                                                                                            aaatttaaaaaac‑CCGCCGGCGATATCAGGGCAGGGCGTCTTGCTGCCCATCATTCCCCCCG‑GGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGagag                                                         <  1:992405/141‑1 (MQ=255)                                                                                            aaatttaaaaaac‑CCGACGGCGAAAAAAAGGCAGGGCGTCTTGCTGCCCATCATTCCCCCCG‑GGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGc                                                        <  2:1118513/142‑1 (MQ=255)                                                                                                ttAACAAAC‑CCGACGGCAATATAAGGGCAGGGCGTTTTGCTGCCCATCATTCCCCCCA‑CGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                     <  1:956055/150‑1 (MQ=255)                                                                                                 tATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                                    >  2:654307/1‑150 (MQ=255)                                                                                                   aaaaaacccgcgggCGATAACAAGGCAGGGCGTTTTGCTGCCCATCATTCCCCCCG‑GGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAt                                                  <  2:36468/139‑1 (MQ=255)                                                                                                            gACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGt                                         >  1:310819/1‑151 (MQ=255)                                                                                                               gggaatatAAGGGCAGGGCGACTTGCTTCCCACCGGTCCCCCCA‑CGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAgg                                            <  1:1007614/141‑1 (MQ=255)                                                                                                                gCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCtt                                     >  1:493528/1‑151 (MQ=255)                                                                                                                 cGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCtt                                     >  1:788814/1‑150 (MQ=255)                                                                                                                   atatCAAGGCATGGCGCTTTGCTGCCCATCATTCCCCCCG‑GGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCtt                                     <  1:601863/148‑1 (MQ=255)                                                                                                                     ataagggCAGGGCTTCTTGCTCCCCATCATTCCCCCCA‑CGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTg                                    <  2:370984/142‑1 (MQ=255)                                                                                                                     atCAAGGCAGGCCGTCTTGCTGCCCATCATTCCCCCCA‑CGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTcc                                <  1:95802/151‑1 (MQ=255)                                                                                                                        agggCAGGGCGTCTTGCTGTCCTTCAGTCCCCCCA‑GGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAtt                             <  1:1127728/149‑1 (MQ=255)                                                                                                                        aaGGCAGGGCGTTTGGTTATCCATCAGTCCCCCCA‑GGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCt                                  <  2:1135367/146‑1 (MQ=255)                                                                                                                          ggcggggcgctttGCTGTCAATCAGTCCCCCCG‑CGGAATTCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCa                               <  2:1128984/136‑1 (MQ=255)                                                                                                                                gCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGttt                     >  2:648298/1‑151 (MQ=255)                                                                                                                                 cGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGttt                     >  2:387718/1‑150 (MQ=255)                                                                                                                                  gACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                   >  1:713613/1‑151 (MQ=255)                                                                                                                                     ttGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc                >  2:114404/1‑151 (MQ=255)                                                                                                                                      ttgctccccaccgTTCCCCCAGCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                   <  1:404942/136‑1 (MQ=255)                                                                                                                                      tGCTCTCCATCAGTCCCCCCA‑CGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTta                    <  1:435055/146‑1 (MQ=255)                                                                                                                                         tGTCCACCAGTCCCCCCC‑CGGGAAGCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                   <  2:99149/144‑1 (MQ=255)                                                                                                                                            ccATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTg         >  1:732813/1‑151 (MQ=255)                                                                                                                                                 aGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGcc     >  1:1116047/1‑150 (MQ=255)                                                                                                                                                 aGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTAAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGcc     >  2:926366/1‑150 (MQ=255)                                                                                                                                                    tCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGc  >  1:482498/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                       TGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGC  >  NC_000913/3578141‑3578436

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: ‑