Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,578,301	0	A	C	55.7%	61.2 / 99.0	62	G142G (GGT→GGG)	gntR	d‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (25/2);  new base C (0/34);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 3.99e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-11

TCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTG  >  NC_000913/3578152‑3578444                                                                                                                                                        |                                                                                                                                                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)                                                                                                                                         cccccca‑cGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTgcg               <  2:366199/141‑1 (MQ=255)                                                                                                                                          ccccca‑gGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTg                 <  1:337271/139‑1 (MQ=255)                                                                                                                                          ccccca‑cGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAg           <  1:1133501/145‑1 (MQ=255)                                                                                                                                           cccccgcGGGAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGcc             <  2:211642/140‑1 (MQ=255)                                                                                                                                           cacca‑cGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATg       >  2:655129/1‑151 (MQ=255)                                                                                                                                             cca‑gGGGAATGCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGt     <  2:189750/151‑1 (MQ=255)                                                                                                                                                  gggAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTg  >  2:217211/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                TCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCA‑CGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTG  >  NC_000913/3578152‑3578444

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