Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,234,243	0	A	C	42.6%	78.0 / 28.4	61	P162P (CCA→CCC)	pgi	glucosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (8/27);  new base C (26/0);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.63e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGACGTAATGCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTC  >  NC_000913/4234096‑4234387                                                                                                                                                    |                                                                                                                                                  cAAAGACGTAATGCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATAGCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATAc                                                                                                                                                <  1:627804/151‑1 (MQ=255) cAAAGACGTAATGCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGCTCTTACCTCCGGCCaca                                                                                                                                                 >  1:867189/1‑148 (MQ=255)   aaGACGTAATGCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACAt                                                                                                                                              <  1:739098/151‑1 (MQ=255)    aGACGTAATGCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCCGCCCCTACa                                                                                                                                               >  1:406846/1‑149 (MQ=255)     gACGTAATGCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGGTCTGCCCCCCGCCCCCACATgg                                                                                                                                            >  2:197855/1‑151 (MQ=255)         tAATGCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTTACCTCCGCCCCCACAAGGGGAc                                                                                                                                        >  1:156390/1‑151 (MQ=255)           aTGCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGGCCCCGGCCCCAACATgggtgcc                                                                                                                                       >  2:239472/1‑145 (MQ=255)             gCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGaa                                                                                                                                    <  2:444673/151‑1 (MQ=255)             gCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGaa                                                                                                                                    <  2:627279/151‑1 (MQ=255)                ggAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTTACCCCGGCCCCTACAAGGTGACCGaa                                                                                                                                    >  1:1029080/1‑148 (MQ=255)                      cAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCg                                                                                                                            <  2:415478/150‑1 (MQ=255)                       aaCGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑Tc                                                                                                                         <  1:414997/151‑1 (MQ=255)                           cgGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTa                                                                                                                     <  2:440508/151‑1 (MQ=255)                                cTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTTACCTCGGCCCATACATAGTGACCAAAAGTCTGCG‑TCCGGACaaa                                                                                                                 >  2:684544/1‑150 (MQ=255)                                 tGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCCTACATGGGGACCGAAGCCCTGCG‑TCCGTCCaaaa                                                                                                                >  1:233785/1‑150 (MQ=255)                                 tGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCCGCCCCTACATGGTGACCGAAGCTCTGCG‑GCCGTACaaaa                                                                                                                >  1:615893/1‑150 (MQ=255)                                  ggAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACaaaaa                                                                                                               <  1:58013/150‑1 (MQ=255)                                  ggAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCCGCCCCCCCCTGGGGACCGCAGCCCTTCG‑CCCGTACaaaaa                                                                                                               >  1:985609/1‑150 (MQ=255)                                   gagaAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCCCCGCCCATACAAGGGGACCGAAGCTCTTCG‑TCCGTACAAAAacc                                                                                                             >  1:191883/1‑151 (MQ=255)                                     gaagaTGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAaccac                                                                                                           <  1:204466/151‑1 (MQ=255)                                      aagaTGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCCGCCCCTACATGGGGGCCGAAACTCTGCG‑TCCCTGCaaaaaac                                                                                                             >  2:539338/1‑146 (MQ=255)                                      aagaTGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGCTCTCACCCCGGCCCCTACATGGGGACCCCAGATCTGCG‑GCCCTCCAAAAaacacc                                                                                                          >  1:756058/1‑151 (MQ=255)                                       agaTGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTTACCCCCCGCCCTACATAGTTACCACAACGCTGCC‑TCCCCACaaaaaaaacc                                                                                                          >  1:448147/1‑145 (MQ=255)                                       agaTGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCCTACAAGGTGACCGAAGCTCTCCG‑GCCGCACAAAAaacacc                                                                                                          >  1:155297/1‑150 (MQ=255)                                       agaTGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCCGCCCCCACAAAGGGGCCGAAGCGCCGCG‑TTCCTACaaaaa                                                                                                               >  2:97177/1‑145 (MQ=255)                                        gaTGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGGCCCCCGCCCCCAAAAGGGGACCCAAGCCCTGCG‑GCCCGCCaaaaaac                                                                                                             >  2:995549/1‑144 (MQ=255)                                         aTGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGGTTTTACCCCGGCCCATACATGGGGACCGAAGCGATGCG‑TCCGCACAAAAaccacc                                                                                                          >  2:470777/1‑148 (MQ=255)                                           gAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGaa                                                                                                      <  2:3412/150‑1 (MQ=255)                                             aaaCCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTTTGACCTCGGCCCCTACACGGGGGCCGGAGCGCTGCG‑TCCGTACaaaaaaacaccca                                                                                                       >  2:119497/1‑139 (MQ=255)                                              aaCCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATg                                                                                                  <  1:631842/151‑1 (MQ=255)                                               aCCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTTTTACCTCCGCCCCCACATGGTGACCGAAACTCTTCG‑GCCCTTCAAAAAACACCCCAACAttc                                                                                                 >  1:844171/1‑149 (MQ=255)                                                 cTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTGTGACCCCCCCCCCCAAATGGTGGCCGAAACTCCTCG‑TCCGCTCAAAAAACACCCTAACAAGCa                                                                                                >  2:184453/1‑150 (MQ=255)                                                 cTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCCTAAATGGGGGACGAAGCCCTGCG‑TCCGGACAAAAaacaccccgaaatgc                                                                                                 >  2:283796/1‑140 (MQ=255)                                                         aaGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTTAACTCGGCCCCTACATGGTGACCGAAGCTCTGCG‑TCCGTTCAAAAACCACCCGAACAAGCACTTTGttt                                                                                        >  2:12122/1‑150 (MQ=255)                                                         aaGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCTGCCCCTACAAGGGGACCGCAGGTCTGCG‑TCCGTACAAAAAACACCTGAACATACACTtttttt                                                                                        >  1:839356/1‑150 (MQ=255)                                                         aaGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGttt                                                                                        <  2:707495/150‑1 (MQ=255)                                                         aaGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTc                                                                                       <  1:952588/151‑1 (MQ=255)                                                         aaGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTc                                                                                       <  1:662398/151‑1 (MQ=255)                                                              attattTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACg                                                                                  <  2:10549/151‑1 (MQ=255)                                                                  tttCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCCCCGCCCCCACATTGTGACCGAAGCCCTGCG‑TCCGGACAAAAACCACCCTAACATGCACCTTTTTTTTTACGTCgg                                                                              >  1:886039/1‑150 (MQ=255)                                                                      cGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTTTGACCCCCGCCCACAAATAGTGGACGAAACACTTCG‑GCCGGACAAAAACAACCCGAAAATACACTTTTTTTTTTACGTCGGTggg                                                                          >  2:71899/1‑151 (MQ=255)                                                                        gTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGCCCTCCGCCCCAACCTTGGGAACGAAACTCTGCG‑TCCCGGCAAAAAAAACCATAACATGGACTTTTTTTTTAACGTTGATggg                                                                          >  1:674590/1‑149 (MQ=255)                                                                             tGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTTAACCCGGCCCACACAAGGTGAGCGAAGACCTTCG‑GCCGGACAAAAAACACCTGAACAAGCACCTTTTTTTTTACGTCGATggggatca                                                                     >  1:972779/1‑144 (MQ=255)                                                                               gAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATc                                                                 <  1:750656/151‑1 (MQ=255)                                                                               gAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGGTCTGAACCCGGCCCCTACATGGGGACCGCAGCGCTGCG‑TCCGTTCAAAAAACACCCGAAAAAATACTTTTTTTTTAACGTCGATGGGGGTcac                                                                    >  2:320358/1‑148 (MQ=255)                                                                                 aaGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATcgc                                                               <  1:918761/151‑1 (MQ=255)                                                                                   ggTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTTTTACCCCCGCCCCTACAAGGGGAACCAAGCCCTGCG‑GCCGTACAAAAAACACCCTAACATACACTTTTTTTTTTACGTCCGTGGGGGTCCCATCCCgg                                                             >  2:130663/1‑151 (MQ=255)                                                                                    gTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTTACCTCCGCCCCTACATGGGGACCGAAGCTCTTCG‑TCCGCACAAAAAACACCCCAAAAAGCACTTTTTTTTTTACGGCGGTGGGGCTCACATCCCgg                                                             >  1:57619/1‑150 (MQ=255)                                                                                       ataCCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAg                                                          <  2:798003/150‑1 (MQ=255)                                                                                       ataCCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAg                                                          <  2:448125/150‑1 (MQ=255)                                                                                       ataCCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGt                                                         <  1:958582/151‑1 (MQ=255)                                                                                                   cAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGt                                             <  1:865699/151‑1 (MQ=255)                                                                                                   cAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGt                                             <  1:409217/151‑1 (MQ=255)                                                                                                   cAATCACTGACGTAGTGAACATCGGGATCGGCGGCTTTTTCCTCCGCCCCCACATGGTGAGCGAAGCTCCCCGTTCCGT‑CAAAAAACACCCCAACAAGCACTTTTTTTCTTACCTCCCTGGGGCACCCACCACGGGAGGGCTGgaaaaa                                               >  2:651167/1‑149 (MQ=255)                                                                                                     aTCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAg                                              <  2:209332/148‑1 (MQ=255)                                                                                                     aTCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTaa                                           <  1:333141/151‑1 (MQ=255)                                                                                                       cACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTaaa                                          <  1:283812/150‑1 (MQ=255)                                                                                                                aGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAAc                                 <  2:577426/150‑1 (MQ=255)                                                                                                                                ggcggTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGt                  >  2:419245/1‑149 (MQ=255)                                                                                                                                      tCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCt           >  2:842270/1‑150 (MQ=255)                                                                                                                                              cGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCtt   <  1:1041549/150‑1 (MQ=255)                                                                                                                                               ggCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTc  >  2:748906/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CAAAGACGTAATGCCGGAAGTCAACGCGGTGCTGGAGAAGATGAAAACCTTCTCAGAAGCGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCG‑TCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGCGGAAGTGCTGAAAAAAGTAAACCCGGAAACCACGCTGTTCTTGGTAGCATCTAAAACCTTC  >  NC_000913/4234096‑4234387

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