Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	343,444	0	A	G	40.0%	99.8 / 16.9	70	G187G (GGA→GGG)	yahK	broad specificity NADPH‑dependent aldehyde reductase, Zn‑containing
Reads supporting (aligned to +/- strand):  ref base A (9/33);  new base G (28/0);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 5.92e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGA  >  NC_000913/343299‑343592                                                                                                                                                  |                                                                                                                                                      atGGTCTGCGTATTCGTCACCCGCAATGGCAGCTGGCGGCGGTGGCTCCTTGGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACata                                                                                                                                                   <  2:925205/150‑1 (MQ=255)   gTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATg                                                                                                                                                 <  2:1049822/150‑1 (MQ=255)   gTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCCGGGGCATATgg                                                                                                                                                >  1:974855/1‑151 (MQ=255)   gTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGAGGGCGTGGTCGGCATCGGCGGTGTGGGGCattttg                                                                                                                                                >  1:1072386/1‑147 (MQ=255)    ttCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCCGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATtgg                                                                                                                                               >  1:65835/1‑151 (MQ=255)       tGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTCGGGCATATtggg                                                                                                                                              >  1:1325206/1‑149 (MQ=255)              tCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGTTAAGCT‑g                                                                                                                                    >  1:989597/1‑151 (MQ=255)               cGTCACCCGCACGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGTTGGGGCAAATtggggtgtaggcgg                                                                                                                                    >  2:225065/1‑141 (MQ=255)               cGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑g                                                                                                                                    <  2:568077/150‑1 (MQ=255)               cGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCCTATGGGGGTTAAACT‑g                                                                                                                                    >  2:28232/1‑150 (MQ=255)                 tCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGc                                                                                                                                  <  1:1244676/150‑1 (MQ=255)                  cACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGccc                                                                                                                                <  2:936767/151‑1 (MQ=255)                   aCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCAATTGGGGGTTGAACC‑GGccc                                                                                                                                >  1:1159388/1‑150 (MQ=255)                     ccGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGAAAAAAAGGGGGCGGGGTCGGCATCGGCGGACTTGGGGATATaggggtataagtttgccac                                                                                                                               >  2:987036/1‑135 (MQ=255)                       gCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGGAAAAAAGTGGGCGTGGTCGGCATCGGCGGGGTGGGGGATATggggggttaactgggcc                                                                                                                                 >  2:939664/1‑133 (MQ=255)                        cAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCATATTGGGGGTAAACG‑GGGCCCcgcg                                                                                                                           >  2:647938/1‑150 (MQ=255)                        cAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCATATGGGGATTAAGCT‑GGccc                                                                                                                                >  1:382425/1‑145 (MQ=255)                          agagCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATg                                                                                                                        <  2:1204313/151‑1 (MQ=255)                           gagCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATgg                                                                                                                       <  1:305813/151‑1 (MQ=255)                               agcTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCAGCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGc                                                                                                                   <  2:119634/151‑1 (MQ=255)                                  tGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTTGGACAAAAGGGGGTTTAACT‑GGGCCACGCGGTggggg                                                                                                                    >  1:337479/1‑147 (MQ=255)                                    gcggcgGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACAtg                                                                                                              <  1:1170203/151‑1 (MQ=255)                                         ggTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACTTgtggtg                                                                                                         <  2:1317918/151‑1 (MQ=255)                                         ggTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATgtggtg                                                                                                         <  1:5183/151‑1 (MQ=255)                                           tGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCATATGGGGGTTAAGCT‑GGCCCACGCGCTGGGGGGACAAGtggggg                                                                                                        >  2:1018941/1‑150 (MQ=255)                                            ggCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGc                                                                                                       <  1:382726/150‑1 (MQ=255)                                             gCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGCGTGGGGGATATGGGGGTTAAAGT‑GGGCCACCCGAAGGGGGCACATggggttgga                                                                                                      >  1:523044/1‑146 (MQ=255)                                                    tgttgtGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATTGGGGTTTAGCT‑GGGCCCCCCGCTGGGGGGCCACGtggtggtcttttcc                                                                                                >  2:1238805/1‑141 (MQ=255)                                                         tgtgCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGCGCGTGGTCGGCATCGGCGGACTGGGGGAAATTGGGGTGTATGT‑CGCCCACGCGATGGGGGCACATGTGGTGGCCTTTaaccatt                                                                                            >  1:547503/1‑142 (MQ=255)                                                               gggATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCa                                                                                    <  2:759131/150‑1 (MQ=255)                                                                 gATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGTTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGGATTTACCCCTTCTGAGGCaaa                                                                                  >  1:1285705/1‑150 (MQ=255)                                                                    caccacGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGAGATATGGGGGGTAAGCC‑GGCCCCCGCGATGGGGGCACATGTGGTGGCATTTTCCACCTCTGAGGCaaaa                                                                                 >  2:393321/1‑148 (MQ=255)                                                                         cGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGTTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAAAGCGaa                                                                           >  1:545704/1‑149 (MQ=255)                                                                         cGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGGGTAAAAC‑GGGCCCCACGGTGGGGGGCCCAGTGGCGGCATTTTCCACCTCTGGGGCaaaa                                                                                 >  1:893654/1‑143 (MQ=255)                                                                          gTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGTTAAGCT‑GGCCCACGCGATGGGGGGACATGTGGGGGCATTTACCACTTCTGAGGCAAAACGCGAAg                                                                          >  1:456987/1‑149 (MQ=255)                                                                          gTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCg                                                                        <  2:1267597/151‑1 (MQ=255)                                                                            aTTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGACATATGGGGGTTTAGCC‑GGCCCACGCGATGGGGGCACATGTGGGGGCATTTACCACTTCTGAGGGAAAAAGCGAAGCgg                                                                       >  2:494925/1‑150 (MQ=255)                                                                              tCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCCGGGGAATATGGGGATTTAACT‑GGGCCCCGCGGTGGGGGCACATGTGGGGGGATTTTTCACCTCTGAGGGAAAAACCGAAgggg                                                                       >  2:240689/1‑148 (MQ=255)                                                                                gccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGAATATGGGGGTTTAACT‑GGCCCCCCCGATGGGGGGACATGTGGGGGGATTTTCCACTTCTGAGGCAAAAAGCGAAGCGGCaaa                                                                   >  2:872022/1‑150 (MQ=255)                                                                                 ccgcTACGTCACTGGCAGGCCGGGCGGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGGCAATGGGGATTACGCTCGGCCCCGCCGTTTGGGGGA‑ACGTGGGGGGATTTATCACCTCTTAGGGAAAAagcggatgggagaa                                                                   >  2:790730/1‑139 (MQ=255)                                                                                 ccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGCCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAg                                                                 <  1:293280/151‑1 (MQ=255)                                                                                  cgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGc                                                                <  1:1332020/151‑1 (MQ=255)                                                                                      aCGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTg                                                            <  1:53253/151‑1 (MQ=255)                                                                                      aCGACACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCt                                                             <  2:998740/150‑1 (MQ=255)                                                                                       cGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGTTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAAAGCGAAGCGGCCAAAGCCCCgg                                                           >  1:993752/1‑151 (MQ=255)                                                                                         tCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGCGGTCGGCATCGGCGGACCCGGGCATATGGGGGTTAAACT‑GGCCCCCGCGCCGGGGGGCCATGTTGGGGGATTTTACACTTCTCGCGCAAAACCCGAAGCGGCGAAAggactgtgg                                                         >  2:1277226/1‑143 (MQ=255)                                                                                          cACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTggggg                                                        <  2:1069123/151‑1 (MQ=255)                                                                                          cACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTgggg                                                         <  2:773112/150‑1 (MQ=255)                                                                                              ggCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGcc                                                      <  2:1159325/149‑1 (MQ=255)                                                                                                cAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGGTTAAGCT‑GGGCCACGCGGTGGGGGCCCATGTGGGGGGATTTACCACCTTTGAGGCAAAAAGCGAAGCGGGAAAAGCCCTGGGGGCCGAt                                                   >  2:545866/1‑150 (MQ=255)                                                                                                      gggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCCTATGGGGGATAAGCT‑GGCCCACGCGATGGGGGCACAAGTGGTGGGATTTACCACTTTTGGGGCAAACAGCGAAGGGGCAAAAGCCCTGGGGGCCGATGAAgt                                              >  2:447235/1‑149 (MQ=255)                                                                                                       ggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgt                                           <  1:845387/151‑1 (MQ=255)                                                                                                            ggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTTGGAGATAATGGGGTTAAGCT‑GGGCCACGCGAAGGGGGGCAAAGTGGGGGGCTTTACCAATTCTGAGGCAAAACAGGAAACGGCAAAAGCCCTGGGGGCAGATGAAGTTGTTaaa                                       >  1:381894/1‑149 (MQ=255)                                                                                                            ggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTa                                         <  1:174982/148‑1 (MQ=255)                                                                                                                aaaaaGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCCTATGGGGATTAAACT‑GGCCCCCCCGCTGGGGGCGCATGTGGTGGCATTTACCACTTCTGAGGGAAAACGCGAAGCGGCAAAAACCCTGGGGGCCGATGAAGGTGTTAACCCACg                                  >  2:18322/1‑151 (MQ=255)                                                                                                                 aaaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGc                                 <  1:327123/151‑1 (MQ=255)                                                                                                                       gggCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGcc                           <  1:110317/151‑1 (MQ=255)                                                                                                                        ggCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCg                          <  2:625545/151‑1 (MQ=255)                                                                                                                        ggCGTGGTCGGCATCGGCGGACTTGGGCCTATTGGGGATAAACG‑GGGCCCCCCGACGGGGGGACATGTGGTGGCATTTACCACTTCTCAGGCAAAAAGCGAAACGGCAAAAACCCTCGGGGGCGATGAAGGTGTTAAATCAAGCAATgtc                           >  2:995261/1‑148 (MQ=255)                                                                                                                         gCGTGGTCGGCATCGGCGGTGTGGGGCATATTGGGGTTAAGCT‑GGGCCCCGCGGTGGGGGGACATGTGGGGGCATTTACCACTTCTGAGGCAAAAAGCGGAGCGGCAAAAGACCTGGGGGCCGGTGGAGTTGTTAACTCACGc                                 >  2:585699/1‑143 (MQ=255)                                                                                                                           gTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATg                       <  2:932677/151‑1 (MQ=255)                                                                                                                           gTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATg                       <  1:1286471/151‑1 (MQ=255)                                                                                                                               tCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGGACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATg                       >  1:825785/1‑147 (MQ=255)                                                                                                                               tCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATgaga                    <  2:267523/150‑1 (MQ=255)                                                                                                                                       ggcggTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTc           <  2:656457/151‑1 (MQ=255)                                                                                                                                         cggTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCAt         <  2:1108123/151‑1 (MQ=255)                                                                                                                                          ggTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATc        <  1:28232/151‑1 (MQ=255)                                                                                                                                             cTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTg      >  1:1239996/1‑150 (MQ=255)                                                                                                                                               gggACATATGGGGATTAAGCA‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGAGGCAAAAGCGCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCGCATCTGa     >  2:868355/1‑149 (MQ=255)                                                                                                                                                ggACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCGCGCAATGCCGATGAGATGGCGGCTCATCTGAaga  <  2:996869/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      ATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCT‑GGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGA  >  NC_000913/343299‑343592

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