Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	47.5%	75.8 / 12.7	61	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (20/12);  new base C (0/29);  total (20/41)
Fisher's exact test for biased strand distribution p-value = 3.78e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786706‑787005                                                                                                                                                      |                                                                                                                                                       cGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAc                                                                                                                                                       <  1:519954/151‑1 (MQ=255)               aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAcc                                                                                                                                         <  2:492913/151‑1 (MQ=255)               aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAc                                                                                                                                          <  2:576367/150‑1 (MQ=255)                      gAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccgccg                                                                                                                                     >  2:512013/1‑148 (MQ=255)                       aaTCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCgctg                                                                                                                                  >  1:639826/1‑150 (MQ=255)                            ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                                            >  2:930194/1‑151 (MQ=255)                            ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                                            >  2:166678/1‑151 (MQ=255)                            ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                                            >  1:465388/1‑151 (MQ=255)                            ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                                            >  2:392790/1‑151 (MQ=255)                                 ggCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCg                                                                                                                       >  1:185194/1‑151 (MQ=255)                                         aatccgcccggggccGCATCAGGCATTGGGCTCAATTGTTAGCAACAAAAAACCCGACTCCCTTCCAGCCGGTTTTCTCATTTAAAACGAATTCCTTTTCCTTCGCTTTCCCCGGTTTGGAAACCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                               <  2:1020450/137‑1 (MQ=255)                                           gTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAg                                                                                                             >  2:874792/1‑151 (MQ=255)                                            ccgcgGAGGGTCGCATCAGGCATTTTGCTCCATTTTTACAAACAAAAAAGCGGACTCCCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTCCTTCGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGc                                                                                                            <  1:166639/150‑1 (MQ=255)                                              ccgcAGCGTCGCATCAGGCATTGTGCTCCATTTTTATCAACAAAAAAGCCTCCTCACTGGCTTTCGGTTTTCCCTTTTAAAACGAATGGCTTTTCCTTCCCTTTCCCCTGGTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt                                                                                                          <  1:274261/150‑1 (MQ=255)                                               cggggggTCGCATAAGTCAATTTGCCCAATTTTTAGCAAAAAAAAAGCCGCCTCCCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt                                                                                                          <  2:52115/144‑1 (MQ=255)                                               cgcAGCGTCGCATCAGCCAATGGGCTCCATTGTTAGCAAAAAAAAGGCCGACTCATTGGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt                                                                                                          <  2:51648/150‑1 (MQ=255)                                                    gcccGCATCGGGAATGGTGCCCCATTGTAAGAAACAAAAAACCGGTTCCCCTTGCAGTCGCGTTTCTCATTTAAACCTAATACCTTTTCCTTCCTTTCCCCCTGGTTTGCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     <  1:304582/147‑1 (MQ=255)                                                    ctgcgcctgaGGCAGTGTGCCCCATTGTAAGCAAAAAAAAAGCCCGCTCACGTGCAGTCGTCTTTCTCTTTTAAAACAAATGACGTTTCCTTCCCTTTACCCTGGTTGGCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     <  1:613770/141‑1 (MQ=255)                                                    cGTCGCATCGGGCATTGTGCTCCATTTTTAAAAACAAAAAACCCGACTCACTTGCATTCGGTTTCCGCATTTAAACCAAATGACTTTTCTTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     <  1:178781/150‑1 (MQ=255)                                                        gCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                >  2:997192/1‑151 (MQ=255)                                                        ccATCAGGCATTGTCCTCCATTGTTAACAAAAAAAAGGCCGACTCACTTGCAGTCGTTTTTCTCATTTTAAACGAATGCCTTTTACTTCGCTTTCCCCTGTTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                <  1:91535/150‑1 (MQ=255)                                                        ccATCAGGCAAGGTGCTCCATTGTTAGAAAAAAAAAAGCCGACTCATTTGCAGTCGGTTTTCTCATTTTAAACGAATGCCTTTTACTCCGCTTTCCCCTGTTTGGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                <  1:685218/150‑1 (MQ=255)                                                         caacgggCAATGTGCTCAATTGTTAAAAACAAAAAGCCCGACTCACTTGCGGGCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTACCCCTGGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                <  1:1046144/145‑1 (MQ=255)                                                         cATCGGGCAATGTCCTCCATTTTTAGAAAAAAAAACCCCGCCTCCCTTGCAGTCTGTTTTCCCATTTTAACCGAATGCCTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                <  2:258155/150‑1 (MQ=255)                                                         cATCAGGCATTTTGCTCATTTTTTAGCAACAAAAAAGCCGACTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGCCTTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                <  2:25792/150‑1 (MQ=255)                                                         aaTCAGGAAAGGTGCTCCATTGTTAGCAACAAAAAACCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACAAATGCCGTTTACTCCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                <  1:565621/149‑1 (MQ=255)                                                          ctCAGGCAATTTCCTCCATTTTTACCAAAAAAAAACCCGCCTCACTTGCAGCCGGCTTTCTCATTTAAAACAAATGCCTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              <  1:585532/150‑1 (MQ=255)                                                               gCAATTTGCTCCATTGTTAAAAAAAAAAAAGCCGCCCCCCTTGCAGTCGGCTTTCTCATTTTAAACGATTGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                                          <  2:162182/150‑1 (MQ=255)                                                               gCAAGGTGCTCCTTTGTTAAAAACAAAAAAGCCGACTCACTTGCAGGCGTTTTTCTCATTTTAACCGAATGACGTTTCCTTCGCTTTACCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                                          <  1:958536/150‑1 (MQ=255)                                                               acAATGTGCTCCATTTTAAGAAAAAAAAAACCCGACTCCCTTGCAGTCGGCTTTCTCATTTAAAACGAATACCTTTTCCTCCCCTTTCCCCGGTTTCGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAg                                                                                            <  1:273184/147‑1 (MQ=255)                                                                aaatggccTCCATTTTAAAAAACAAAAAAGCCGACTCCTTTGCTGTCGGCTTTCTCATTTAAACCGATTGACTTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                                         <  2:850316/143‑1 (MQ=255)                                                                       ccccATTGTAAACAACAAAAAACCCGTCTCCCTGCCTGCCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCCCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                                 <  2:190865/149‑1 (MQ=255)                                                                         ccATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGtt                                                                               >  1:154761/1‑151 (MQ=255)                                                                          cATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                              >  1:855841/1‑151 (MQ=255)                                                                           aTTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                             >  1:865170/1‑151 (MQ=255)                                                                               ttAGCAAAAAAAAAGCCGACTCTCTTGCAGTCGGTTTTCCCATTTAAAACAAATGACTTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAAt                                                                         <  1:109958/151‑1 (MQ=255)                                                                                      caaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                   >  1:807450/1‑150 (MQ=255)                                                                                                  ctctctTGCGGTCGGTTTTCTCATTTAAAACAAATGACGTTTACTTCTCTTTCCCCTGGTTGTCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                      <  2:952934/151‑1 (MQ=255)                                                                                                  cTCCCTTGCAGTCGTCTTTCCCATTTTAAACGAATGACGTTTCCTTCGCTTTACCCTGGTTCGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                      <  1:364750/151‑1 (MQ=255)                                                                                                  cTCCCTTGCAGTCGGTTTTCTCATTTTAAACGTATGACGTTTCCTTCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                      <  2:381585/151‑1 (MQ=255)                                                                                                    cttttcctggcGGCTTTCTAATTATAAGGAAATGACTTTTCCTTCGTTTACCCCTGGTGTGAATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     <  2:710732/140‑1 (MQ=255)                                                                                                    cACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCa                                                    >  2:216861/1‑151 (MQ=255)                                                                                                    cACTTGCAGTCGGCTTTCTCATTCTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCa                                                    >  2:522972/1‑151 (MQ=255)                                                                                                        tgcatccGGTTTTCTCATTTAAAACAAATGCCGTTTTCTCCGCTTTCCCCGGTTTTCAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                <  2:986234/145‑1 (MQ=255)                                                                                                         ccAGGCGGCTTTCTAATTAAAAGCGAATGACGTTTCCTCCGCTTCCCCCGGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                <  1:185802/149‑1 (MQ=255)                                                                                                              cGGTTTTCTCATTTAAAAAGAATGCCTTTTACTTCGTTTTACCCTGGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGc                                           <  2:302873/150‑1 (MQ=255)                                                                                                                           ttAAACGAGTGCCGTTTACTCCGCTTTACCCTGGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                             <  1:843945/151‑1 (MQ=255)                                                                                                                           ttAAACGAATGACTTTTCCTCCCTTTTCCCCTGTTTTGAACCCGCCGCTTCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                             <  2:510567/151‑1 (MQ=255)                                                                                                                           ttAAACGAATGACGTTTACTTCGCTTACCCGTGGTTAGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                             <  1:372692/151‑1 (MQ=255)                                                                                                                            tAAACAAATGACTTTTCCTTCG‑TTTACCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                           <  2:717544/151‑1 (MQ=255)                                                                                                                                cGAATGACGTTTACTTCGTTTTCCCCTGGTTTGCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGTGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc                        <  2:593301/151‑1 (MQ=255)                                                                                                                                     gACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGaa                   >  2:523826/1‑151 (MQ=255)                                                                                                                                      aCGTTTCCTTCGTTTTCCCCTTGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGaa                   <  1:883523/150‑1 (MQ=255)                                                                                                                                       cttttCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtc                 <  2:864356/149‑1 (MQ=255)                                                                                                                                       cGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt                  >  2:820352/1‑150 (MQ=255)                                                                                                                                         tttCCTTCCCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtct                <  2:301584/150‑1 (MQ=255)                                                                                                                                               tccttttCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc         <  2:931198/147‑1 (MQ=255)                                                                                                                                               tCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTcttctt          >  2:886067/1‑150 (MQ=255)                                                                                                                                               tCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTcttctt          >  1:165247/1‑150 (MQ=255)                                                                                                                                                  cTTTACCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGc       <  2:154797/150‑1 (MQ=255)                                                                                                                                                    ttACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCa    >  1:777573/1‑151 (MQ=255)                                                                                                                                                      aCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  1:1045104/1‑151 (MQ=255)                                                                                                                                                      aCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  2:473718/1‑151 (MQ=255)                                                                                                                                                      aCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGCCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCCCAAGAATCTCTTCTTCGCTCacg  >  2:907695/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786706‑787005

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