Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,033,102	0	T	C	41.3%	95.4 / 13.4	63	S322P (TCC→CCC)	hyaA	hydrogenase 1, small subunit
Reads supporting (aligned to +/- strand):  ref base T (32/5);  new base C (0/26);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.76e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGAT  >  NC_000913/1032953‑1033245                                                                                                                                                       |                                                                                                                                                ataACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCAtt                                                                                                                                                 <  2:421001/150‑1 (MQ=255)    aCGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCGTTCCAc                                                                                                                                             <  1:607158/151‑1 (MQ=255)     cGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCAcc                                                                                                                                            >  2:284350/1‑151 (MQ=255)         tGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCAccccc                                                                                                                                         >  2:144466/1‑150 (MQ=255)             cctccACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATAc                                                                                                                                    >  2:433237/1‑151 (MQ=255)                   cacGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCgg                                                                                                                              >  1:1032532/1‑151 (MQ=255)                         ggAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAAc                                                                                                                        >  1:1167638/1‑151 (MQ=255)                         gaaaTGAGGGTTTTTCTTCCCCATTCCAGTCGGGTCGCGGCTGCCTGGGCTGGGGGGAAATTGGTTTCTGGGATGGGGTTTCTTTCCACAGCCGGGTGGTGGATATTCCGAAAAGGGG‑TCCTCATCCCACCGCCGATACCGTCGGttt                                                                                                                           <  2:647797/146‑1 (MQ=255)                          gAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAAcc                                                                                                                       >  2:902282/1‑151 (MQ=255)                             tgatgGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACcgcg                                                                                                                    >  1:432873/1‑151 (MQ=255)                               atgGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCCCCGATACCGTCGGTTTAACCGCGCt                                                                                                                  >  2:603662/1‑151 (MQ=255)                                  gCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTgg                                                                                                               >  1:114353/1‑151 (MQ=255)                                           tcCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGc                                                                                                      >  1:292238/1‑151 (MQ=255)                                             ccAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGc                                                                                                      >  2:398613/1‑149 (MQ=255)                                                aTCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCgg                                                                                                 >  2:830510/1‑151 (MQ=255)                                                 tCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCgg                                                                                                 >  1:52045/1‑150 (MQ=255)                                                  cgaGTGTGATCACGTCGGTCTGGGCTTGGGGAAAAATGTTTTCTGGGATCGCGTTTCTTTCTACACCGGCGGGTTCTATTTCCCGAAAGTGGG‑TACGCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAg                                                                                                    <  2:808163/144‑1 (MQ=255)                                                      tCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTg                                                                                           >  1:779490/1‑151 (MQ=255)                                                         ggTCGCGGCTCCGTGGGCTGTGCGGAAAATGTTTTGTGGGCTCGGGGTTCGTTCTACACCCGCGTGGTCGATATCCCGAAAAGGGT‑TACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGgt                                                                                         <  1:198901/150‑1 (MQ=255)                                                         ggTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGgtg                                                                                        >  2:1205553/1‑151 (MQ=255)                                                           tcccgtcTGCCGGGGCTGTGCGAAAAATGTTTTCTGGGATCGCGGTTCGTTCTACGGCCGCGTGGTCGATATCCCGCAAAGGGG‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTgg                                                                                          <  2:637525/141‑1 (MQ=255)                                                           tCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGCTCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGtgtg                                                                                      >  2:140472/1‑151 (MQ=255)                                                            cACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGc                                                                                     >  2:357897/1‑151 (MQ=255)                                                                    ccggggCGGTGGGAAAAGTGGTTTCGGGGATCGGGGTTCGTTCTACGGCCGGGGGGTCGATATCCCGCAAGGGGG‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCAc                                                                                   <  2:288997/142‑1 (MQ=255)                                                                       gggcggggcgAAAATTGTTTTCGGGGATCGGGTTTCTTTCTCCAGCCGGGTGGTCGATTTTCCGCAAGGGGG‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCg                                                                            <  1:407821/149‑1 (MQ=255)                                                                       gggCTGGGCGAAAAATGTTTTCGGGGACCGCGTTTCTTTCTACAGCCGCGTGGTCAATATCCCGCAAATGGG‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCg                                                                            <  2:662108/149‑1 (MQ=255)                                                                         gCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCTGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGcc                                                                          >  1:167032/1‑149 (MQ=255)                                                                            gggggaaaaaGGGTTTCTGGGATCGCGGTTCTTTCTACAGCGGCGGGGTCGATATCCCGCAAAGGGG‑TCCTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAgc                                                                       <  1:921141/143‑1 (MQ=255)                                                                                gaaaaaTGTTTTCGGGGATCGCGGTTTTTTCTACAGCCGCGTGGTCGATATTCCGCAAAGGGG‑TACTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGtt                                                                 <  1:872211/149‑1 (MQ=255)                                                                                   aaaTGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGAc                                                              >  1:1126927/1‑151 (MQ=255)                                                                                       gtttggtGGGACCGCGGTTCTTTCTACAGCCGGGTGGTCAATATCCCGCAAGGGGG‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCa                                                            <  2:1213800/143‑1 (MQ=255)                                                                                              gggATCGCGTTTCTTTCTACGGCCGGGTGGTCGATATCCCACAAGGGGG‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGAc                                                    <  2:48692/150‑1 (MQ=255)                                                                                               ggATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACg                                                   >  2:787206/1‑150 (MQ=255)                                                                                               ggATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACg                                                   >  2:341735/1‑150 (MQ=255)                                                                                                gATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCACCGCCGTTGACCAGCGCAGACggc                                                 >  1:137165/1‑149 (MQ=255)                                                                                                  tCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCAt                                               >  2:546443/1‑151 (MQ=255)                                                                                                   cgcgGTTCTTTCTACAGCCGCGTGGTCGATATCCCGAAAATGGT‑TACTCTTTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCAt                                               <  1:237833/150‑1 (MQ=255)                                                                                                    ggggttgtttcTACCGCCGTGTGGTGGATTTTCCGAAAATGGT‑CACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATa                                              <  1:951329/142‑1 (MQ=255)                                                                                                          cgttcTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGc                                        >  1:414023/1‑150 (MQ=255)                                                                                                            ttctccgcccGGGTGTTCGATATCCCGCAAGGGGT‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCa                                       <  2:62501/141‑1 (MQ=255)                                                                                                                 cacccggttGGTCGATTTCCCGAAAGTGGT‑TACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACa                                <  2:1167452/143‑1 (MQ=255)                                                                                                                    ccGCGTGGTCTATATCCCGCAATTGGT‑TCCTCATCCCCCCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAg                               <  2:912850/149‑1 (MQ=255)                                                                                                                    ccGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCAACGCCGATACCGTCGGTTTAACCGAGCTTGGTGTGGTGGCAGCGGCTGTTGGTGTGCCCGCATTCCGCAGCGCCGTTGAGTAGCGCAGGAGTCAAAACCAGCAACCTACAGa                              >  2:868479/1‑150 (MQ=255)                                                                                                                      gcgTGGTCGATTTCCCGAAAGGGGG‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAg                               <  2:614642/147‑1 (MQ=255)                                                                                                                         tGGTCGATATCCCGCAAGGGGT‑TACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAAcc                          <  1:10006/149‑1 (MQ=255)                                                                                                                          ggTCGATTTCCCGAAAATGGG‑TACTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGa                        <  1:323284/150‑1 (MQ=255)                                                                                                                          ggTCAATATCCCGCAATTGGT‑TACTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGaa                       <  1:214463/151‑1 (MQ=255)                                                                                                                            tcaatatCCCGCAAATGGG‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAAc                      <  1:1129948/147‑1 (MQ=255)                                                                                                                             cGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACAt                    >  2:78325/1‑151 (MQ=255)                                                                                                                                tatTCAGCAAAGGGT‑TCCTCATCCCCCCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACa                     <  2:964933/147‑1 (MQ=255)                                                                                                                                   tCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGcc               >  1:433029/1‑150 (MQ=255)                                                                                                                                   cccacAAAGGGT‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCa              <  2:1107304/147‑1 (MQ=255)                                                                                                                                    cccaaaaggggttccTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCa              <  2:139875/137‑1 (MQ=255)                                                                                                                                    ccGAAATTGGT‑TACTCTTTCCGCCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCa                  <  2:615930/146‑1 (MQ=255)                                                                                                                                     cacAAATGGGTTCCTCATCCCCCCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCAAACATCAGcc               <  2:504207/147‑1 (MQ=255)                                                                                                                                     cGCAAATGGG‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAgg            <  1:541218/151‑1 (MQ=255)                                                                                                                                       cAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGc           >  1:231313/1‑150 (MQ=255)                                                                                                                                        aaaTGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGACACCGAACATCAGCCAGGCa          >  2:85651/1‑150 (MQ=255)                                                                                                                                           tGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATg       >  1:766074/1‑150 (MQ=255)                                                                                                                                            ggg‑TCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAg     <  1:158035/151‑1 (MQ=255)                                                                                                                                            ggg‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATg       >  1:356968/1‑149 (MQ=255)                                                                                                                                             ggttccTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAg     <  1:193323/146‑1 (MQ=255)                                                                                                                                                tACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGAt  >  2:607074/1‑151 (MQ=255)                                                                                                                                                tACTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGa   <  1:916714/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                ATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGG‑TACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGAT  >  NC_000913/1032953‑1033245

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