Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	40.3%	109.1 / 12.0	62	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (29/8);  new base C (0/25);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.21e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTC  >  NC_000913/786706‑787002                                                                                                                                                      |                                                                                                                                                    cGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAc                                                                                                                                                    >  1:373303/1‑151 (MQ=255)  gTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAcc                                                                                                                                                   >  2:1038135/1‑151 (MQ=255)     agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg                                                                                                                                                >  2:899692/1‑151 (MQ=255)     agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTg                                                                                                                                                <  1:27092/151‑1 (MQ=255)               aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAcc                                                                                                                                      >  1:141139/1‑151 (MQ=255)                       aaTCCCCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCATTGTCCTCCATTGTTAGCAAAAAAAAAGCCGCCTCCCTTGCAGTCGGTTTTCTCATTTAAACCGAAGGACGTTTCCTGCGCTTTCCCCGGGTTTGCAACCGCCgctgc                                                                                                                              <  1:806177/151‑1 (MQ=255)                         tCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCtt                                                                                                                            >  1:762334/1‑151 (MQ=255)                            ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCg                                                                                                                         >  1:714781/1‑151 (MQ=255)                            ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTCCTCCATTGTTAGCAACAAAAAACCCGACTCCCTTGCAGTCGGCTTTCTCATTTTAACCGAATGACTTTTACTTCGCTTTCCCCTGTTTTGAAACCGCCGCTGCttt                                                                                                                           <  1:1007704/149‑1 (MQ=255)                                   ccTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGAt                                                                                                                  >  1:347517/1‑151 (MQ=255)                                       aTAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCCCTGGCAGTCGGCTTTCTCATTTTAAACGAATGCCGTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                               <  1:953942/150‑1 (MQ=255)                                        tAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                             >  2:173296/1‑151 (MQ=255)                                         aaGTCGCGCAGCGTCGCATCAGGCATTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt                                                                                                             <  2:363534/150‑1 (MQ=255)                                          aGTCGCGCAGCGGCGCATCAGGCAATGTGCCCCATTGTTAGCAACAAAAAAGCCGCCTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGCCGTTTACTCCCCTTTACCCTGTTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCa                                                                                                           <  2:110053/151‑1 (MQ=255)                                           gTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACCCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAg                                                                                                          <  2:769493/151‑1 (MQ=255)                                                 cAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAc                                                                                                    >  2:789239/1‑151 (MQ=255)                                                  aGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGTATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc                                                                                                   >  2:969413/1‑151 (MQ=255)                                                   gCGTCGCATCAGGCAATGTGCTCCATTTTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGTCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAc                                                                                                    <  1:832821/149‑1 (MQ=255)                                                   gCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                  <  2:690223/151‑1 (MQ=255)                                                    cGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                  >  2:418169/1‑150 (MQ=255)                                                      tCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                  >  1:148194/1‑148 (MQ=255)                                                         cATCAGGCAATGTCCTCCATTTTTAGAAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGCCGTTTACTTCGCTTTCCCCTGTTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCa                                                                                                 <  1:891745/146‑1 (MQ=255)                                                         cATCAGACAATGTGCCCCATTTTGAGCAAAAAAAAAGCCGACCCACTGGAAGTCGGTTTTCTCTTTTTAAACGAATAACTTTTACTTCGTTTGCCCCTGTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                             <  1:155151/150‑1 (MQ=255)                                                            cAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                           >  1:29425/1‑149 (MQ=255)                                                             aGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                        >  1:654687/1‑151 (MQ=255)                                                              ggCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataata                                                                                          <  2:666210/148‑1 (MQ=255)                                                                 aaTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGttt                                                                                    >  2:722623/1‑151 (MQ=255)                                                                  aTGTGCTCAATGGTTAGCAACAAAAAAGCCGACTCCCTTGCAGTCGGTTTTCTCTTTTAAAACGAATGCCGTTTCCTCCGCTTTCCCCTGTTTTAAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                                      <  1:150725/148‑1 (MQ=255)                                                                     tacactaTTGTAAACAAAAAAAAACCCGTCTCACTTGCTGTAGGTTTTCGCATTTAAAACAATCAACGTTTCCTGCGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAg                                                                                         <  1:882999/136‑1 (MQ=255)                                                                       cTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                              >  2:818097/1‑151 (MQ=255)                                                                            ttGTTAGCAACAAAAAACCCGACTCCCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCgg                                                                              <  1:28132/146‑1 (MQ=255)                                                                             tGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                          >  2:504023/1‑149 (MQ=255)                                                                               ttAGCAACAAAAAAGCCGACTCCCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaa                                                                        <  2:395352/149‑1 (MQ=255)                                                                               ttAGCAACAAAAAAGCCGACTCCCTTGCAGTCGGCTTTCTCTTTTAAAACAAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaaa                                                                       <  2:957120/150‑1 (MQ=255)                                                                                   caacaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtctc                                                                  >  2:44170/1‑151 (MQ=255)                                                                                    aacaaAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtctc                                                                  <  2:63550/150‑1 (MQ=255)                                                                                    aacaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcg                                                                 >  1:978663/1‑151 (MQ=255)                                                                                       aaaaaaGCCGACTCCCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                              <  2:745545/151‑1 (MQ=255)                                                                                       aaaaaaGCCGACTCACTTGCAGTCGTCTTTCTCATTTAAAACAAATGACGTTTACTTCGCTTTCCCCTGTTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                               <  1:770540/150‑1 (MQ=255)                                                                                          aaaGCCTACTCCCTTGCGTTCGTCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTACCCGGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                           <  1:968702/151‑1 (MQ=255)                                                                                          aaaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                           >  1:993766/1‑151 (MQ=255)                                                                                          aaaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                           >  2:890859/1‑151 (MQ=255)                                                                                           aaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACt                                                          >  1:238496/1‑151 (MQ=255)                                                                                              ccGACTCACTGGCAGTCGGCTTTCTCATTTTAACCGAATGACGTTTACTTCGCTTTACCCTGTTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                       <  1:815985/151‑1 (MQ=255)                                                                                                    cACTTGCATTCGCCTTTCTCATTTTAACCGAATGACGTTTCCTTCGTTTTCCCCTGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATacac                                                   <  1:777748/149‑1 (MQ=255)                                                                                                        tGCGGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                              <  1:359684/150‑1 (MQ=255)                                                                                                         gctgTCGGTTTTCTCATTTAAACCGAATGACTTTTCCTTCGCTTTACCCTGGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                             <  1:789600/147‑1 (MQ=255)                                                                                                           aGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                          >  2:317779/1‑151 (MQ=255)                                                                                                               ggCTTTCTCATTTTAACCGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                          <  1:487018/147‑1 (MQ=255)                                                                                                                    tctcATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCgg                                 >  1:188759/1‑151 (MQ=255)                                                                                                                    tctcATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCgg                                 >  2:910483/1‑151 (MQ=255)                                                                                                                     ctcATTTTAAACGAATGACGTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                  <  2:430717/149‑1 (MQ=255)                                                                                                                         ttttAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGa                               <  1:44163/148‑1 (MQ=255)                                                                                                                          tttAAACGAATGACTTTTACTTCGCTTTCCCCGGTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGatat                            <  1:872120/150‑1 (MQ=255)                                                                                                                                cGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc                     >  2:80625/1‑151 (MQ=255)                                                                                                                                   aTGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                    >  1:491323/1‑149 (MQ=255)                                                                                                                                         tttCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt               <  1:728175/148‑1 (MQ=255)                                                                                                                                         tttACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt               <  2:820968/148‑1 (MQ=255)                                                                                                                                              ttCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtctt          <  1:437444/148‑1 (MQ=255)                                                                                                                                                cGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttct        <  1:462662/148‑1 (MQ=255)                                                                                                                                                cGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCg     >  1:419959/1‑151 (MQ=255)                                                                                                                                                   tttACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc  >  2:344228/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTC  >  NC_000913/786706‑787002

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