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

ATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGC  >  NC_000913/3578140‑3578430                                                                                                                                                     |                                                                                                                                                atgatAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTcca                                                                                                                                              <  1:175399/151‑1 (MQ=255)         aCGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACggg                                                                                                                                       >  1:1076649/1‑150 (MQ=255)           gatcgtcgAGACGTGCGCCGAGTTAGGCAATGTGTGGATGCCGGCGAGCAATTATGGCAGTGGTCATCTGGCGTGCTGACTCAAAGCTATCAAAA‑CCGACGCCGATATCAAGACTGGGCGTCTTGCTGTCCATCAGTTCGACGACCGGGAt                                                                                                                                    >  1:406019/3‑151 (MQ=255)            ttCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAAt                                                                                                                                    >  1:288692/1‑150 (MQ=255)             tcgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATa                                                                                                                                   <  2:254540/150‑1 (MQ=255)             tcgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATAc                                                                                                                                  >  1:510371/1‑151 (MQ=255)               gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATAccc                                                                                                                                >  2:395380/1‑151 (MQ=255)                tcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCg                                                                                                                               >  1:363821/1‑151 (MQ=255)                 cgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGc                                                                                                                              >  1:942797/1‑151 (MQ=255)                 cgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGc                                                                                                                              >  1:36480/1‑151 (MQ=255)                            ccGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatc                                                                                                                    >  2:223368/1‑150 (MQ=255)                            ccGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatc                                                                                                                    >  2:537957/1‑150 (MQ=255)                             cGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatcat                                                                                                                  >  1:1007568/1‑151 (MQ=255)                             cGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatcat                                                                                                                  >  2:729100/1‑151 (MQ=255)                                     gCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                                                          >  1:693879/1‑151 (MQ=255)                                      cAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                                                          >  2:47072/1‑150 (MQ=255)                                      cAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                                                          >  1:657013/1‑150 (MQ=255)                                                  gCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGCTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGCGCTCGGCgtgt                                                                                              >  2:58846/1‑150 (MQ=255)                                                        gAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACg                                                                                       >  1:278005/1‑151 (MQ=255)                                                        gAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACg                                                                                       >  2:115335/1‑151 (MQ=255)                                                         aaaaaaaatGGTAGTGGTCGTCGGGGGTGCTGCTAAAAATTAAAAAAAC‑CCGGCGGCAATATCAGGGCAGGGCGCTTTGCTTTCCATCATTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTAc                                                                                        <  2:576051/143‑1 (MQ=255)                                                            aataatGGCGGTGTTCTTCTGGCGTGCTGCTTAAAAGTTAAAAAAC‑CCGACGGCAATATCAGGGCAGGGCGCCTTGCTGTCCATCATTCCCCCCACGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTa                                                                                         <  2:755609/145‑1 (MQ=255)                                                            aataatGGATGTGGTCATCGGGCGTGCTCCTTAAAATTTAACAACA‑CCGAGGGCAATATCAGGGCAGGGCTACTTGCTCCCCATCTGTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                                                                    <  2:914637/150‑1 (MQ=255)                                                               aatGGCAGGGGTCTTCGGTCGTCCTGCTACAAATTAAAAAAAC‑CCGCCGGTATTATCAGGGTAGGGCGTTTGGCGTCCCATCATTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                 <  2:820580/150‑1 (MQ=255)                                                               aagggCGGTGGTCATGGGGCGTGCTGCTT‑AAAATTAACAAAACCCGGCGGCAATACCAAGGCGTGGCGCTTTGTTGTCCATCAGTCCCCCCACGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                 <  1:415621/147‑1 (MQ=255)                                                                   gCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAgg                                                                             >  2:571493/1‑150 (MQ=255)                                                                      gTGGTCATCTGGCGTGCTGCTACAAATTTAACAAAC‑CCGCCGGCAATATCAAGGCAGGGCGTCTTGTTTTCCATCCGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGa                                                                            <  1:250197/148‑1 (MQ=255)                                                                      gTGGTCATCGGGCGTGCTGCTAAAAAGTTATAAAAC‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGCCCATCAGTCCCACCACGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                          <  2:517347/150‑1 (MQ=255)                                                                      gTGGTAATCGGGCGTGTTGCTAAAAAGTTAAAAAAC‑CCGACGGCGATATAAAGGCAGGGCGATTTGCTGTCAATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                         <  2:1153176/151‑1 (MQ=255)                                                                        gttcattgggCGTGCTGCTT‑AAAATTAACAAAACCCGACGGCAATACAAGGGCGGGGCTCTTTGTTGTCCACCGGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                          <  1:537140/140‑1 (MQ=255)                                                                        ggTCATGGGGCGTGCTGCTAAAAATTTAAAAAAC‑CCGACGGCGATATCAAGGCAGGGCGATTTGCTTCCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                         <  1:1065647/149‑1 (MQ=255)                                                                          tCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGcc                                                                     >  1:532220/1‑151 (MQ=255)                                                                          tCATCGGGCTTGCTGCTTAAAATTTAACAAAA‑CCGACGGCGATATAAAGGCAGGGCGTCTTGCTTTCCATCAGTCCCCCCGGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGc                                                                      <  2:1112731/150‑1 (MQ=255)                                                                            aTCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGt                                                                   >  1:599541/1‑151 (MQ=255)                                                                            aTCGGGCGGGCTGTTTAAAAGTTATAAAAA‑CCGCCGGCGATAAAAAGGCGTGGCGTCTTGCTTTCCAGCATTCCCACCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCg                                                                    <  1:833103/150‑1 (MQ=255)                                                                                ggctgccgactaaaaaGTTATCAAAC‑CCGCCGGCGATTTAAAGGCAGGCCGTTTTTTTTCCCATCAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGc                                                                      <  1:229043/131‑1 (MQ=255)                                                                                            cAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATAGCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGCTCGGTGAGGATCAGGCCGTCGATATTCCAGgag                                                     >  1:597337/1‑149 (MQ=255)                                                                                            aaaaTTTATAAAAC‑CCGCCGGCGATACAAAGGCGGGCCGTCTTGTTTTCCATCATTCCCCCCGGGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                       <  2:788605/146‑1 (MQ=255)                                                                                                 ttATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                              >  1:856525/1‑151 (MQ=255)                                                                                                   aaaaaac‑CCGACGGCGATACAAAGGAAGGGCGACTTGCTGCCCATCATTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                              <  2:1084519/146‑1 (MQ=255)                                                                                                   aTCAAAC‑CCGCCGGCGATACCAAGGCAGGGCGTTTTGCTGCCCACAATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCAt                                                <  2:118748/147‑1 (MQ=255)                                                                                                     cAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTc                                          >  2:205051/1‑151 (MQ=255)                                                                                                      aaaa‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCCCCACGGGAATACCCGCCACCTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAtt                                           <  1:330045/149‑1 (MQ=255)                                                                                                         acccGACGGCGATAAAAAGGCATGGCGCCTTGCTGTCCATCAGTCCCCCCGCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAg                                       <  2:808652/149‑1 (MQ=255)                                                                                                             gACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGt                                   >  1:100227/1‑151 (MQ=255)                                                                                                               cGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTc                                 >  2:291938/1‑151 (MQ=255)                                                                                                                ggCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCt                                >  1:679702/1‑151 (MQ=255)                                                                                                                 gCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCtt                               >  2:1088147/1‑151 (MQ=255)                                                                                                                  cGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTg                              >  2:219109/1‑151 (MQ=255)                                                                                                                    atatAAAGGCAGGGCGACTTGCTGTCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCt                            <  2:562988/151‑1 (MQ=255)                                                                                                                     tataAAGGCAGGGCGCTTTGCTGCCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc                             <  2:277803/149‑1 (MQ=255)                                                                                                                     tatCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTc                           >  2:17692/1‑151 (MQ=255)                                                                                                                     taacAAGGCAGGGCGTTTTGCTGTCCATCAGTCCCCCCGCGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc                             <  1:531835/146‑1 (MQ=255)                                                                                                                         aggtcacggcaatttGCTGTCCAGCCGCCCCCCCACGCCATTACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAg                                       <  1:537947/122‑1 (MQ=255)                                                                                                                         aaGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAtt                       >  2:422013/1‑151 (MQ=255)                                                                                                                             ccgcccgatttattGTCCATAATTCCCCCCACGGAAATTCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTc                     <  2:465082/136‑1 (MQ=255)                                                                                                                                   gtcTTGCTGTCCATCAGTCCCCCCGCGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTta              <  2:701432/148‑1 (MQ=255)                                                                                                                                        gCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCg        >  1:961722/1‑151 (MQ=255)                                                                                                                                        gCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCg        >  1:267958/1‑151 (MQ=255)                                                                                                                                          tctCCATCAGTCCCCCCGCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAcc         <  2:940515/146‑1 (MQ=255)                                                                                                                                            tCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGt    >  1:1126247/1‑151 (MQ=255)                                                                                                                                               aTCAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTgc  <  1:189242/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                ATGATAGTACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGC  >  NC_000913/3578140‑3578430

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