Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,133,302	0	T	G	42.4%	96.8 / 13.6	59	H219P (CAC→CCC)	yghS	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (7/27);  new base G (25/0);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.14e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGG  >  NC_000913/3133153‑3133449                                                                                                                                                      |                                                                                                                                                     aaTGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg                                                                                                                                                     <  1:654538/151‑1 (MQ=255)    gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTc                                                                                                                                                   <  1:734354/150‑1 (MQ=255)    gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCg                                                                                                                                                  <  2:942733/151‑1 (MQ=255)    gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCg                                                                                                                                                  <  1:590692/151‑1 (MQ=255)       gTTGGGTTGCGGCGCGGTTGTACTATTAAAAGGTTGAGATTGGAGTGCATCCATTCGTTAAATCATCCCATATCCATTTTGGTAATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtggctacc                                                                                                                                               >  1:391189/1‑145 (MQ=255)           ggCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                                                                                           <  1:403497/151‑1 (MQ=255)              tGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAACTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGACTTGCGCGTAGTGACGGTTGCAGTCTTTCATCGGCGGGGTTTCGcgccacc                                                                                                                                           >  1:499499/1‑142 (MQ=255)                cggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTTGcccccccacttt                                                                                                                                      >  1:725270/1‑142 (MQ=255)                        ttGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGCCCATCCCGTTCAAgg                                                                                                                                  >  1:991869/1‑147 (MQ=255)                          gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGGTCttt                                                                                                                            >  2:342814/1‑151 (MQ=255)                          gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGGTCAAGGATCttt                                                                                                                            >  2:378538/1‑151 (MQ=255)                           tACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGGTCaaaggttctt                                                                                                                            >  2:188556/1‑142 (MQ=255)                            aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGTGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCCGTTCAAGGATCttt                                                                                                                            >  1:86731/1‑149 (MQ=255)                             cTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTAAAGATCttt                                                                                                                            >  1:224618/1‑148 (MQ=255)                               aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGcc                                                                                                                       <  2:1006439/151‑1 (MQ=255)                               aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGcc                                                                                                                       >  2:549889/1‑151 (MQ=255)                                 tAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGcc                                                                                                                       >  2:330508/1‑149 (MQ=255)                                  aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGCTCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTgcg                                                                                                                         <  1:75896/146‑1 (MQ=255)                                  aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCAt                                                                                                                     <  2:625545/150‑1 (MQ=255)                                  aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCAt                                                                                                                     <  2:889727/150‑1 (MQ=255)                                  aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCAt                                                                                                                     <  1:561336/150‑1 (MQ=255)                                  aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGCCCCTCCAGTTCTAGGATCTTTTCGCCCtt                                                                                                                    >  1:287607/1‑151 (MQ=255)                                        ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATg                                                                                                               <  2:50725/150‑1 (MQ=255)                                        ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTAAGGGTCTTTTCGCCCTTGAAt                                                                                                                >  1:737241/1‑149 (MQ=255)                                         tgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCCTTTAAt                                                                                                                >  2:1061776/1‑148 (MQ=255)                                               tgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTg                                                                                                        <  2:1003289/150‑1 (MQ=255)                                                gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcg                                                                                                      <  2:330968/151‑1 (MQ=255)                                                gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAAGGTCTTTTCGCCCTTTAATGTCCGtt                                                                                                         >  2:1012842/1‑148 (MQ=255)                                                 tAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAGGATCTTTTCGCCATTTAATGTCAGTTgcgg                                                                                                     >  2:907284/1‑151 (MQ=255)                                                  aGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcgg                                                                                                     <  2:1103448/150‑1 (MQ=255)                                                        tCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTCAAAGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGGaa                                                                                                >  2:305620/1‑149 (MQ=255)                                                            ttCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGTTTAAGGATTTTTGCGCCATTTAATGTCAGGTGCGGCGGAA‑CGGCg                                                                                          >  2:1100878/1‑150 (MQ=255)                                                             tCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGGTGCGGCGTAA‑CGGCGAt                                                                                        >  2:559999/1‑151 (MQ=255)                                                              cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGAt                                                                                        <  1:363457/150‑1 (MQ=255)                                                              cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCCCCCCTCCCGGTTAAGGAGCTTTTCGTCATTGTATGTCACTTGCGGGGGAA‑CCGCGGtt                                                                                       >  1:947145/1‑151 (MQ=255)                                                                  aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGAttttt                                                                                    <  1:208174/150‑1 (MQ=255)                                                                  aaaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGAttttt                                                                                    <  2:351902/150‑1 (MQ=255)                                                                                  aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGTTTAAGGAACTTTTCGCCATTGAATGTAAGGTGCGGCGGAACCGCCGATTTTTT‑CTGTaaa                                                                           >  1:399612/1‑142 (MQ=255)                                                                                        gTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTtggtg                                                             <  2:286437/151‑1 (MQ=255)                                                                                           aTCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTc                                                          <  2:583796/151‑1 (MQ=255)                                                                                               ggAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGt                                                      <  1:1012693/151‑1 (MQ=255)                                                                                                gAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTTAAAGGTCTTTTTGCCATTGAATGTCAGGTGCGGCGGAA‑CGGCGGTTTTTTTCTGTAACGCTGCCCGGTGGGGGTCAgg                                                       >  1:939717/1‑149 (MQ=255)                                                                                                  aaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGGAA‑CGGCGGTTTTTTTCTTTTACGCTGCCAGGTGGTGGGCAGGttt                                                    >  2:56717/1‑150 (MQ=255)                                                                                                   aGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTAc                                                  <  2:538168/151‑1 (MQ=255)                                                                                                   aGCGCGGCGTGAATTGCGCGTAGTGACGCTTGAAGAATTTCATCGGCGGGGGGGCGGCCTACCAGTTAAAGGAAGTTTGTGCCCTTTAAATTCAGTTGGGGGGTAA‑GGGGGTTTTTTTCTTTGAAAGCGGCCCGTTTGGGGGCGGGttt                                                    >  2:22973/1‑149 (MQ=255)                                                                                                      gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTAcgcg                                               <  2:1016319/151‑1 (MQ=255)                                                                                                          cgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAGGTCAAAGATCTTTGCGCCATTGAATGTCAGGTGCGGCGGAA‑CGGCGATTTTTTTCCTTAACGCTGCCAGGTGGGGGTCAGGGTTACgcgcca                                            >  2:123428/1‑148 (MQ=255)                                                                                                            tGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTCCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGGAA‑CGGCGATTTTTTTCTGTAACGCTGCCAGTTGGGGGGCAGGGTTTCGCGCcaa                                           >  2:397072/1‑149 (MQ=255)                                                                                                             gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCCTTTAATGTCAGTTTCGGCGGAA‑AGGCGGTTTTTTTCTGTAACGCTGCCAGGTGGTGGGCAGGGTTACGCGCaaa                                           >  1:447413/1‑148 (MQ=255)                                                                                                               aTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCgg                                      <  2:1021952/151‑1 (MQ=255)                                                                                                                 tGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTTCGCCCTTGAATGTCAGTTGCGGCGGAA‑CGGCGGTTTTTTTCTGTAACGCTGCCCGTTGGTGGGCAGGTTTACGCGCCAACGCgggt                                    >  2:609903/1‑149 (MQ=255)                                                                                                                   cgcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTg                                  <  2:2908/151‑1 (MQ=255)                                                                                                                      gTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTCAGGAACTTTTCGCCCTTGAATGTCAGTTGTGGCGTAA‑AGGCGATTTTTTTCTGTAAAGCTGCCCGGTGGTGGGCCGGGTTTCGCGCAAAAGCGGgttgtt                                >  1:753453/1‑150 (MQ=255)                                                                                                                              gCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTTCTGTAACGCTGCCAGTTGGTGGGCAGGGTTACGCGCAAAAGCGGGTTGTTCATCAATGc                       >  2:551618/1‑151 (MQ=255)                                                                                                                                  gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCaa                    <  1:416557/150‑1 (MQ=255)                                                                                                                                  gCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCaa                    <  2:95931/150‑1 (MQ=255)                                                                                                                                     gAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGa                  <  2:615428/149‑1 (MQ=255)                                                                                                                                       aTTTCATCGGCGGGGGGTCGCCCATCCAGGTCAAGGGTCTTTTCGCCATTGAATGTCAGTTGCGGGGGAA‑CGGCGGTTTTTTTCTTTAAAGCTGCCAGGTGGGGGGCAGGTTTACGCGCAAACGCGGGTTGTTCATCAATGCCAAGa                  >  2:402129/1‑147 (MQ=255)                                                                                                                                         ttCATCGGCGGGGGGTCGCCCATCCCGCTTAAAGGTCTTTGCGCCATTGTATGTCAGTTGCGGCGTAA‑CGGCGGTTTTTTCCTGTAACGCTGCCAGTTGGGGGTCAGGGTTTCGCGCAAAAGCGGTTTGTTCATCAATGCCAAGACgca              >  1:657228/1‑147 (MQ=255)                                                                                                                                                 cggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAc    >  1:64271/1‑151 (MQ=255)                                                                                                                                                   gggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACgg  <  2:64616/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAA‑CGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGG  >  NC_000913/3133153‑3133449

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