Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,585,086	0	A	C	41.8%	127.2 / 15.8	67	Y580D (TAC→GAC)	ggt	gamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (37/2);  new base C (0/28);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 7.42e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGG  >  NC_000913/3584937‑3585233                                                                                                                                                       |                                                                                                                                                    ttattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAc                                                                                                                                                    >  2:962075/1‑151 (MQ=255)  tattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcc                                                                                                                                                   >  1:1113661/1‑151 (MQ=255)           aGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt                                                                                                                                           >  2:6185/1‑150 (MQ=255)                       aaTTCTTTCCCGGGGACCGAAAAAACCGCTGGAAAAATA‑CCAAAGGAATACAAAAAAAAAAAAGGCTCCCTTGGCCTTCCCCTAAAAAAATAGCCCTCTCCCCACGAAGAGGGCCGCTACCTTTATTCCCCCGCCGTTAAATCATCCACCg                                                                                                                              <  2:747069/151‑1 (MQ=255)                             ttcTCGGGGACCGGAAAAACCGCTGAAAAAATATCTACAGGATTAC‑AAAAAAAAAAAGCCTCCCTTCGGCTTGCCCTGAAAAAATAGCCCTCTCCCCACGAAGAGGCCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAgc                                                                                                                           <  1:745212/148‑1 (MQ=255)                                tctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGt                                                                                                                     >  2:279476/1‑151 (MQ=255)                                  tGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                                                   >  1:513287/1‑151 (MQ=255)                                  ggggAACGGAAAACCCGCTGGAAAAATATCTACAGGATTAA‑AAAAAAAAAAAGGCCCCTTTCGGCTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTACCTTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTc                                                                                                                    <  2:884526/149‑1 (MQ=255)                                   gggAACGCAAAACCCGCTGAAAAAATATCTACAGGATTAC‑AAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCCCTCCCCACGAAGGGGCCCGCTACCCTTATTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTc                                                                                                                    <  1:399181/149‑1 (MQ=255)                                     gAACGGCAAAACCGCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGa                                                                                                                 >  2:695967/1‑150 (MQ=255)                                     gAACGGCAAAACCGCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGAt                                                                                                                >  1:65768/1‑151 (MQ=255)                                                  gCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                    >  2:257757/1‑150 (MQ=255)                                                  gCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                   >  1:663302/1‑151 (MQ=255)                                                     ggAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                                >  1:669729/1‑151 (MQ=255)                                                     ggAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                                >  1:441526/1‑151 (MQ=255)                                                     ggAAAAATATCAACAGAATTAC‑AAAAAAAAAAAGGCTCCCTGCGGCTTCCCCTAAAAAAATAGCCCTCTCCCCACGAAGAGGGCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                                 <  2:967150/150‑1 (MQ=255)                                                      gaagaaTATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                               >  1:859481/1‑151 (MQ=255)                                                      gaagaaTATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                               >  1:601795/1‑151 (MQ=255)                                                       aagaaTATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                              >  1:477184/1‑151 (MQ=255)                                                           atatCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                           >  1:853838/1‑150 (MQ=255)                                                               caaCAGAATTAC‑AAAAAAAAAAGGGCTCCCTGCGGCTTCCCCTGAAAAAATAGCCCTCTCCCCACGAAGAGGGCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                          <  2:518325/145‑1 (MQ=255)                                                                 aCAGTATTAC‑AAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                     <  2:194268/150‑1 (MQ=255)                                                                  cAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                   >  1:637720/1‑151 (MQ=255)                                                                     gaattaaaaaaaaaaaaaGGCTCCCTCCGTCTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                   <  1:516270/141‑1 (MQ=255)                                                                      aTTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTTATGCCACGAAGAGTGCCGCTAACCTTAGTACCCCGCCGTTAAGTCATCCACCGAGCGCGGGTCGGATCCGCCGTACAACTCACCGTCCGGCCCAAGCa                                                                                 >  2:202095/1‑149 (MQ=255)                                                                        ttccaaaaaaaaaaaGGCTCCCTCCGGCTTGCCCTGAAAAAAGAGCCCTCTCCCCAGGAAGAGGCCCGCTACCCTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATa                                                                               <  1:396619/145‑1 (MQ=255)                                                                            aaaaaaaaaaaGGCTCCCTTGGGCTCGCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTAACTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                            <  1:14217/148‑1 (MQ=255)                                                                             aaaaaaaaaaGCCTCCCTTCGCCTTCCCTTAAAAAAATACCCCCCTCCCCCCGAAGAGGCCCCCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGc                                                                           <  1:851580/149‑1 (MQ=255)                                                                              aaaaaaaaaGGCTCCCTTCGGCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAAGAGGCCCGCTAACCTAATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                        <  1:754565/148‑1 (MQ=255)                                                                              aaaaaaaaaGGCTACCTTCGGCTGCCCCTAAAAAAAACCCCCTCTCCCCAAGAAGGGGCCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                        <  1:211831/148‑1 (MQ=255)                                                                               aTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAACAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                        >  2:1013670/1‑150 (MQ=255)                                                                                aaaaaagGGCCACCTGCGCCTTCCCCAGAAAAAATGCCCCTCTCCCCAGGAAGAGGGCCGCTACCCTAAGCACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                      <  2:929607/150‑1 (MQ=255)                                                                                    aaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                  >  1:354264/1‑151 (MQ=255)                                                                                      aGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                                >  2:772247/1‑151 (MQ=255)                                                                                       ggCTCCCTTCGGCTTGCCCTGAAAAAATAGCCCTCTCCCCCGGAAGAGGGCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                               <  2:222918/151‑1 (MQ=255)                                                                                         cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                             >  1:672737/1‑151 (MQ=255)                                                                                                 ggCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                     >  2:505048/1‑151 (MQ=255)                                                                                                   cTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                   >  2:468313/1‑151 (MQ=255)                                                                                                   cTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                   >  1:447041/1‑151 (MQ=255)                                                                                                      gCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                 >  1:1061087/1‑150 (MQ=255)                                                                                                      ccccTAAAAAAATACCCCCCTCCCCACGAAGGGGCCCGCAACCCTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCa                                                <  2:652665/150‑1 (MQ=255)                                                                                                       cccTGAAAAAATCCCCCTCTCCCCACGAAGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                 <  1:817954/149‑1 (MQ=255)                                                                                                       cccTAAAAAAATACCCCTTTCCCCCCGAAGGGGCCCGCTACCTTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                      <  2:549291/144‑1 (MQ=255)                                                                                                       cccAGAAAAAATGCCCCTCTCCCCACGAAGAGGCCCTTTACCTTTGGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                 <  1:93068/149‑1 (MQ=255)                                                                                                           aaaaaaaTAGCCTTCTCCCCAGGAGGGGGGCCGTTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGcc                                           <  2:254739/148‑1 (MQ=255)                                                                                                            aCAAAATACCCCTCTCCCCACGAAGAGGCCCGCAAACCTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                            <  2:915933/149‑1 (MQ=255)                                                                                                                           ccccACGAAGGGCCCCTCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTg                           <  1:697133/150‑1 (MQ=255)                                                                                                                              cACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                        >  1:1085124/1‑151 (MQ=255)                                                                                                                                 gAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                      >  1:1085696/1‑150 (MQ=255)                                                                                                                                  aaGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCca                     <  1:642235/150‑1 (MQ=255)                                                                                                                                  aaGAGGCCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcag                    <  1:403671/151‑1 (MQ=255)                                                                                                                                   agagGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                   >  2:443376/1‑151 (MQ=255)                                                                                                                                   agagGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                   >  1:250388/1‑151 (MQ=255)                                                                                                                                    ggggcccGCAACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagca                  <  1:623019/146‑1 (MQ=255)                                                                                                                                     agGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                   >  2:110918/1‑149 (MQ=255)                                                                                                                                      gggCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCagcagc                >  2:758140/1‑151 (MQ=255)                                                                                                                                      ggccccttacccTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                        <  2:1031954/133‑1 (MQ=255)                                                                                                                                      ggcccccTACCCTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagcag                 <  2:713905/144‑1 (MQ=255)                                                                                                                                      ggcccGTTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCagcagc                <  1:974670/148‑1 (MQ=255)                                                                                                                                        gCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCtt              >  1:507255/1‑151 (MQ=255)                                                                                                                                              acccTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCg         <  1:279511/148‑1 (MQ=255)                                                                                                                                              aaCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt        >  1:198507/1‑151 (MQ=255)                                                                                                                                              aaCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt        >  2:1073279/1‑151 (MQ=255)                                                                                                                                               aCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt        >  1:712967/1‑150 (MQ=255)                                                                                                                                                ccTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTa       <  1:38494/150‑1 (MQ=255)                                                                                                                                                ccTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTAt      >  1:315043/1‑151 (MQ=255)                                                                                                                                                    aGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCgg  >  1:1052524/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                    TTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTAC‑GAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGG  >  NC_000913/3584937‑3585233

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