Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,842,914	0	A	G	44.3%	70.1 / 25.4	61	M182V (ATG→GTG)	gdhA	glutamate dehydrogenase, NADP‑specific
Reads supporting (aligned to +/- strand):  ref base A (1/33);  new base G (27/0);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.77e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGT  >  NC_000913/1842769‑1843057                                                                                                                                                  |                                                                                                                                                  cGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGatga‑tg                                                                                                                                             <  1:911511/151‑1 (MQ=255)  gATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAtga‑tga                                                                                                                                            <  2:765768/151‑1 (MQ=255)      cGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGaaaaa                                                                                                                                        <  2:20431/151‑1 (MQ=255)        aaaGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGG‑TGaaaaa                                                                                                                                        >  2:317736/1‑149 (MQ=255)          aGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGG‑TGaaaaaaagc                                                                                                                                    >  1:541873/1‑147 (MQ=255)               aaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGGAAAAACACTcc                                                                                                                                 >  2:365236/1‑149 (MQ=255)                aaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCcaac                                                                                                                              <  1:962087/151‑1 (MQ=255)                 aGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCcaaca                                                                                                                             <  1:394017/151‑1 (MQ=255)                   cGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGGAAAAGAACTccc                                                                                                                                >  1:604545/1‑145 (MQ=255)                     aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TTGAAAAAATCTCcacc                                                                                                                              >  2:826708/1‑144 (MQ=255)                     aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGA‑GGAAAAAAGACTCcaacaaaa                                                                                                                          >  1:393451/1‑148 (MQ=255)                     aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATAc                                                                                                                         <  2:374374/151‑1 (MQ=255)                       ggTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCg                                                                                                                       <  2:844515/151‑1 (MQ=255)                          gAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCt                                                                                                                    <  1:99595/151‑1 (MQ=255)                               gATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGA‑TGAAAAAGCTCTCCCACAAAACCCCCCGc                                                                                                                  >  1:517498/1‑148 (MQ=255)                                aTGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCt                                                                                                              <  2:868777/151‑1 (MQ=255)                                 tGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGA‑TGAAAAAACTCTCCAACAAAAACGCCTGCGGCtt                                                                                                             >  2:254765/1‑151 (MQ=255)                                  gCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCtt                                                                                                             <  2:201251/150‑1 (MQ=255)                                      ttttGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGA‑TGAAAAAAATCTCCCACAATAACGCCCGCGTCTTCAcc                                                                                                         >  2:1011235/1‑150 (MQ=255)                                        ttGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCgg                                                                                                       <  2:707068/150‑1 (MQ=255)                                        ttGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCAAACAACACCGCCTGCGTCTTCACCgg                                                                                                       <  2:264286/150‑1 (MQ=255)                                             aGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAACACTCCCCCAATAACGCCCGCGTCTTTTCCGGG‑GAggg                                                                                                >  1:688644/1‑151 (MQ=255)                                             aGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGA‑TGAAAAAACTCCCCCACAAAAACCCCCGCGGCTTCAcccggaaagg                                                                                                 >  2:750278/1‑143 (MQ=255)                                              ggCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGG‑AGGAAAAACTCTCCCACAAAACCGCCTGCGGCTTCCCCCGT‑AAg                                                                                                  >  2:196505/1‑148 (MQ=255)                                              ggCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAggg                                                                                                <  2:583527/150‑1 (MQ=255)                                              ggCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGc                                                                                               <  1:96527/151‑1 (MQ=255)                                                 gcTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGG‑TGGAAAAAAACTCCAACAATAACGCCTGCCGCTTTACCGGT‑AAagggc                                                                                              >  2:711056/1‑147 (MQ=255)                                                 gcTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCtt                                                                                            <  2:913698/151‑1 (MQ=255)                                                   tgatgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAAATCTCCCACAAAACCCCCCGCGGCTTTACCCGT‑AAGGGCCttt                                                                                           >  2:359154/1‑150 (MQ=255)                                                   tgatgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGA‑TGAAAAAACTCTCCCACAAAACCGCCCGCGTCTTTACCCGT‑AAggggcctt                                                                                           >  2:149282/1‑145 (MQ=255)                                                    gatgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCttt                                                                                           <  1:760971/149‑1 (MQ=255)                                                    gatgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTc                                                                                          <  1:867947/150‑1 (MQ=255)                                                    gatgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTAGGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCa                                                                                         <  2:205615/151‑1 (MQ=255)                                                      tgaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGGAAAAACTCTCCCACAATACCGCCCGCGGCTTCACCCGT‑AAGGGGCCTTTAtt                                                                                       >  2:641650/1‑151 (MQ=255)                                                       gaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCAttt                                                                                      <  2:857168/151‑1 (MQ=255)                                                           gAACTGTATCGCCCCCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATCCCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTggct                                                                                  <  1:179010/151‑2 (MQ=255)                                                           gAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTgg                                                                                    <  1:684739/149‑1 (MQ=255)                                                            aaCTGTATCGCCCCCTGGGCGCGGATACCGACTTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTggcgg                                                                                 <  2:318869/151‑1 (MQ=255)                                                              cTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGG‑TGAAAAAGATCTTCAAAAAAAACCCCCGCGGCTTCCCCCGG‑CAAGGCCCTTCATTTtgcgg                                                                                 >  1:331709/1‑149 (MQ=255)                                                               tGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAACTCTCCCACAAAAACGCCCGCGGCTTCCCCGGTAAAGGGCCTTT‑ATTTtgcgg                                                                                 >  2:286381/1‑148 (MQ=255)                                                                    cGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAAATCTCCCACAATACCGCCTGCGGCTTCCCCGGG‑AAGGGGCTTTTCTTTtgcgggcggctt                                                                          >  2:132493/1‑143 (MQ=255)                                                                    cGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCtt                                                                          <  2:789452/150‑1 (MQ=255)                                                                        aCCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAACTCTCCCACAAAACCGCCCGCGGCTTCCCCCGG‑AAGGGGCTTTTATTTTGCGGCGGGCTtttt                                                                       >  2:581030/1‑149 (MQ=255)                                                                         ccTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCg                                                                     <  1:36955/150‑1 (MQ=255)                                                                          cTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGA‑TGAAAAAACTCTCCCACAAAAACCCCCGCGTCTTTCCCCGG‑AAGGGGCTTTTATTTTGCGGCAGGCTtttt                                                                       >  1:104082/1‑147 (MQ=255)                                                                               cgcgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGA‑TGAAAAAAGTCTCCCACAAAAACGCCCGCGTCTTTACCGGG‑GAGGGCCTTTTATTTTGCGGCCGTCTttttttccc                                                                  >  1:812173/1‑142 (MQ=255)                                                                                gcgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGaa                                                              <  1:317710/150‑1 (MQ=255)                                                                                       ccGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCg                                                       <  2:490614/150‑1 (MQ=255)                                                                                        cGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAACTCTCCCACAATAACCCCCGCGTCTTCCCCCGG‑AAGGGCCTTTCATTTTGCGGCAGTCTTTTTTCCCCCGGAGCTAcccg                                                      >  2:151514/1‑148 (MQ=255)                                                                                            gTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTa                                                   <  2:115134/149‑1 (MQ=255)                                                                                               ccGGAAGGTGATATCGGGGTTGGTGGTCGTGAGGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATCCCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGt                                               <  2:359688/150‑1 (MQ=255)                                                                                                cGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGA‑TTAAAAAAAACTCCAACAAAACCACCCGCGTCTTCACCGGG‑AAAGGGCTTTCCTTTTGCGGCAGGCTTATTCTCCCCGGAGCTACCGGCCACgggc                                              >  2:373859/1‑148 (MQ=255)                                                                                                        gATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTAt                                      <  2:29167/150‑1 (MQ=255)                                                                                                         atatCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATTGAAAAAG‑TCTCCAACAATAACGCCTGCGGGTTCACCCGT‑AAGGGCCCTTTATTTTGCGGCCGTCTTTTTTCCCCCGGAGCTACCGGGCACGGGCTGGGTTAttt                                    >  2:548817/1‑151 (MQ=255)                                                                                                           atCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAAATCTCCAACAAAACCCCCCGCGTCTTTACCCGG‑AAGGGGCCTTTATTTTGGGGGGGGCTTTTTTGCCCGGGAGCTACCGGCTACGGGCTGGGtttttt                                    >  2:737669/1‑149 (MQ=255)                                                                                                               gggTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAACTCTCCCACAAAAACGCCCGCGTCTTCTCCCGG‑GAGGGGCTTTTATTTTGGGGCGGGCTTTTTTTCCCCGGAGCCACCGGCCACGGGGTGGGtttttt                                    >  2:869539/1‑145 (MQ=255)                                                                                                                ggTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGaa                             <  2:896030/151‑1 (MQ=255)                                                                                                                ggTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGaa                             <  2:907318/151‑1 (MQ=255)                                                                                                                      tggtCGTGAAGTCGGCTTTATGGCGGGGGTGA‑TGAAAAAACTTTTCACCAAAAACCCCCCCGTCTTCTCCCGT‑GAGGGGGTTTCTTTTTGCGGGCGGCTTATTTTCCCCGCAGCCACCGGCTACAGGCTGGTTTTTTTTACAGAAGCa                          >  1:897528/1‑148 (MQ=255)                                                                                                                             gAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTaaaa                 <  2:994656/150‑1 (MQ=255)                                                                                                                                           ggCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGt  <  2:604600/151‑1 (MQ=255)                                                                                                                                           ggCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGt  <  2:393479/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGT‑AAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGT  >  NC_000913/1842769‑1843057

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