Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,578,301	0	A	C	45.0%	92.6 / 12.0	60	G142G (GGT→GGG)	gntR	d‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (29/4);  new base C (0/27);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 3.62e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATAAC  >  NC_000913/3578154‑3578449                                                                                                                                                       |                                                                                                                                                     gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATAccc                                                                                                                                                   >  2:709169/1‑151 (MQ=255)      gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCAc                                                                                                                                              >  1:164750/1‑151 (MQ=255)       aCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACt                                                                                                                                             >  2:338295/1‑151 (MQ=255)            cgcCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAat                                                                                                                                        >  2:939369/1‑151 (MQ=255)             gcCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAatc                                                                                                                                       >  1:852135/1‑151 (MQ=255)              ccGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAatc                                                                                                                                       >  2:716173/1‑150 (MQ=255)                  gaTAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCtt                                                                                                                                  >  2:704381/1‑151 (MQ=255)                      ggCATTGGGGCGATCCCCGCGAGAAATAATGGCAGTGGTCATCTGGGGTGCTGCTTAAAATTTATAAAAC‑CCGACGGCGATATAAAGGCAGGGCGCCTTGCTGTCCATCAGTC‑CCCCCAGGGA‑AATCCCCGCCACTTCAATCATCttt                                                                                                                                 <  1:575838/148‑1 (MQ=255)                         attggggCGATCCCCGCAAGAAAAAAGGGGGGTGGTCTTCTGGCGTCCT‑CCTAAAATTTAAAAAAACCCGCCGGCAATATAAAGGCTGGCCGACTTGCTGTCCACCA‑TT‑CCCCCCCGGGAAATCCCCGCCACTTCAATCATCtt                                                                                                                                  <  2:253637/139‑1 (MQ=255)                              ggCGATCCCCGCGACAAAAAAGGGAGGTGTTCATCGGGGGTGCTGCTCAAAATTTAAAAAAC‑CCGCGGGCAATAAAAGGGCAGGGCGTTTTGTTTCCCATCTTTC‑CCCCCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTg                                                                                                                           <  1:285659/146‑1 (MQ=255)                                     cccgggaaaaaaaagggCGGTGCTCATCGGGGGTTCTGCTAAAAAGTTATAAAAC‑CCGAGGGCGATATCAAGGCAGGGCGTTTTGCTGCCCATCAGTC‑CCCCCACGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGtgtg                                                                                                                <  1:483309/135‑1 (MQ=255)                                      ccGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTggg                                                                                                              >  1:601316/1‑151 (MQ=255)                                           aGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGt                                                                                                         >  2:846088/1‑151 (MQ=255)                                            gCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACg                                                                                                          >  2:444255/1‑149 (MQ=255)                                              aataatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGAGCGCGGCGTGTGGGTACGTTc                                                                                                       >  2:685817/1‑150 (MQ=255)                                              aataagGGAAGTGGTCTTTGGGGGTCCTGCTTAAAATTTACAAAAC‑CCGCGGGCAATATAAAGGCATGGCGATTTGCTGTCCATCAGTC‑CCACCACGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTa                                                                                                            <  1:989096/145‑1 (MQ=255)                                                    ggCGGTGTTCATCTGGGGTGTTTCTTAAAATTTACAAAAC‑CCGACGGCGATATCAGGGCATGGCGTTTTGCTGCCCATCAGTC‑CCACCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                                                                                 <  1:219385/150‑1 (MQ=255)                                                    ggCAGTGGTCATGGGGGGTGCTTCTTAAAATTTAAAAAAC‑CCGACGGCGATATCAGGGCAGGGCGTTTTGCTGTCCATCAGTC‑CCCCCAGGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGa                                                                                                  <  1:991587/149‑1 (MQ=255)                                                      cAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGa                                                                                                  >  2:771481/1‑147 (MQ=255)                                                         tGGTCATCGGGCGTGCTCCTAAAATTTTAAAAAAC‑CCGGCGGCGATATAAAGGCATGGCTATTTGCTGTCCATCA‑TT‑CCCCCCCGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                                            <  2:837689/150‑1 (MQ=255)                                                             cATGTGGCGGGCTGCTTAAAATTTAAAAAAC‑CCGCGGGCAATATAAAGGCAGGGCGACTTGCTGCCCATCTGTC‑CCCCCGGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGcc                                                                                        <  1:134730/150‑1 (MQ=255)                                                             cATCGGGGGTGCTTCTTAAAAGTTACAAAAC‑CCGAGGGCAATAGCAGGGCAGGGCGTTTTGCTGTCCATCAGTC‑CCCCCACGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGc                                                                                         <  1:520443/149‑1 (MQ=255)                                                               tCTGGCGTGCTTCTAAAAATTTAACAAAC‑CCGCCGGCGATACAAGGGCAGGGCCTTTTGCTTCCCATCAGTC‑CCCCCGGGGA‑AACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGcc                                                                                        <  2:574897/148‑1 (MQ=255)                                                                  ggcggcctgctTAAAATTAACCAAAC‑CCGAGGGGAATATAAAGGCATGGCTTTTTGCTGTCCATCAGTC‑CCCCCGCGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTc                                                                                     <  2:739724/142‑1 (MQ=255)                                                                  ggCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat                                                                                  >  1:746227/1‑151 (MQ=255)                                                                       gctgctTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTc                                                                              >  2:986301/1‑150 (MQ=255)                                                                            taaaaaTTTAACAAAA‑CCGGCGGCAATATCAAGGCAGGGCGATTTTCTGTCCATCAGTC‑CCACCACGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCa                                                                            <  2:832912/144‑1 (MQ=255)                                                                               aaaGTTATCAAAC‑CCGACGGCGATATAAGGGCAGGGCGTCTTGTGGCCCACCATTC‑CCCCCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGc                                                                     <  1:923837/151‑1 (MQ=255)                                                                                aaGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGc                                                                     >  1:904134/1‑150 (MQ=255)                                                                                     ataaaac‑CCGCCGGCGATACCAGGGCGTGGCGATTTGCTTTCCATCAGTC‑CCCCCACGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAt                                                               <  1:918619/148‑1 (MQ=255)                                                                                     ataaaac‑CCGACGGCAATACCAAGGCATGGCGTCTTGCTTTCCATCATTC‑CCCCCCCGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                                  <  1:612606/145‑1 (MQ=255)                                                                                     aacAAAC‑CCGACGGAAATATCAGGGCATGCCTCTTTGCTGTCCATCGGTC‑CCCCCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                                  <  2:965931/146‑1 (MQ=255)                                                                                     aTCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                                                 >  2:119221/1‑149 (MQ=255)                                                                                        aaac‑CCGACGGAAATAAAAAGGCGTGGCGTCTTGCTCTCCATCAGCC‑CCACCAGGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAt                                                               <  2:55951/148‑1 (MQ=255)                                                                                               gACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGTTATTCCAGGAGAGCATGGATTCGAGGCGt                                                      >  2:36867/1‑151 (MQ=255)                                                                                                  ggCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCt                                                   >  2:170192/1‑151 (MQ=255)                                                                                                  ggCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCt                                                   >  1:57623/1‑151 (MQ=255)                                                                                                     gATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc                                                >  1:415180/1‑151 (MQ=255)                                                                                                      atatCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCt                                               >  2:877513/1‑151 (MQ=255)                                                                                                       tatCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTc                                              >  1:525446/1‑151 (MQ=255)                                                                                                        atCAAGGCATGGCGCTTTGCTCTCAATGAGTT‑CCACCACGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc                                                <  1:70427/148‑1 (MQ=255)                                                                                                           aaGGCATGGCGACTTGCTGCCCATCAGTC‑CCACCACGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTc                                              <  1:882661/147‑1 (MQ=255)                                                                                                                  ggcgttttGCTGTCCA‑CGGTTCCCACCGCGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAtt                                          <  2:813008/138‑1 (MQ=255)                                                                                                                   gCGTCTTGCTGCCCATCAGTC‑CCCCCACGGA‑AAACCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGttt                                  <  2:231325/151‑1 (MQ=255)                                                                                                                   gCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGttt                                  >  2:121930/1‑151 (MQ=255)                                                                                                                    cgccTTGCTGTCCATCAGTC‑CCCCCGCGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTta                                 <  1:135383/148‑1 (MQ=255)                                                                                                                     gtcTTGCTGTCCATCAGTC‑CCCCCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTta                                 <  1:424015/148‑1 (MQ=255)                                                                                                                     attttgttgCCCATCGGTC‑CCACCACGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTta                                 <  2:713709/147‑1 (MQ=255)                                                                                                                      aCTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                                >  2:427454/1‑150 (MQ=255)                                                                                                                        ttGCTGCCCACCAGTC‑CCCCCCGGGA‑AACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTta                                 <  2:264832/147‑1 (MQ=255)                                                                                                                          gCTGTCCATCAGTC‑CCCCCAGGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAcc                            <  1:417436/150‑1 (MQ=255)                                                                                                                               ccATCAGTC‑CCCCCGCGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAg                        <  2:169976/149‑1 (MQ=255)                                                                                                                                cctcctccccccccacGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAg                        <  2:680240/137‑1 (MQ=255)                                                                                                                                    aGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGcc                  >  1:134579/1‑150 (MQ=255)                                                                                                                                              acGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCt        >  1:432169/1‑151 (MQ=255)                                                                                                                                               cGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTg       >  1:425907/1‑151 (MQ=255)                                                                                                                                                ggg‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCt        >  2:921732/1‑149 (MQ=255)                                                                                                                                                    aaTCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTg       <  1:47576/147‑1 (MQ=255)                                                                                                                                                    aaTCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATa    <  1:939394/150‑1 (MQ=255)                                                                                                                                                     aTACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATAAc  >  1:770330/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                     GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTT‑CCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATAAC  >  NC_000913/3578154‑3578449

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