Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	947,763	C→A	5.6%	P179T (CCT→ACT)	ycaM →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	947,763	0	C	A	5.6%	218.3 / 3.5	71	P179T (CCT→ACT)	ycaM	putative transporter
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base A (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

GCGTTACAGGGGTTAACGCTGGTGCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGAGAAGATTTCTCCTTACGTTAATCAAA  >  NC_000913/947616‑947910                                                                                                                                                    |                                                                                                                                                    gCGTTACAGGGGTTAACGCTGGTGCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCgcgc                                                                                                                                                     <  8:248875/148‑1 (MQ=255)    ttACAGGGGTTAACGCTGGTGCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGTTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCATGTATGTGGCGATGGCGGTAACCGCGCCTg                                                                                                                                                  >  1:245943/1‑148 (MQ=255)        aGGGGTTAACGCTGGTGCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAAtt                                                                                                                                             <  5:16464/149‑1 (MQ=255)         ggggTTAACGCTGGTGCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTa                                                                                                                                            >  8:153048/1‑149 (MQ=255)                aCGCTGGTGCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGt                                                                                                                                     <  3:37959/149‑1 (MQ=255)                 cGCTGGTTCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTg                                                                                                                                    <  4:93863/149‑1 (MQ=255)                      gTGCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCaa                                                                                                                                               <  5:13608/133‑1 (MQ=255)                      gTGCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCaa                                                                                                                                               >  6:13608/1‑133 (MQ=255)                      gTGCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCat                                                                                                                                 <  1:368073/147‑1 (MQ=255)                       tGCTGTTTATCTTCTTTATCTGGGTTGCTTCACGCGGTATGAAATCTCAGAAAATCGACGGAACTGTGGCACGCATTCCAATGTTTGTTATCTAACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCAAAtt                                                                                                                              >  6:179206/1‑149 (MQ=255)                               aTCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACca                                                                                                                       <  4:59239/148‑1 (MQ=255)                               aTCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACca                                                                                                                       >  3:229727/1‑148 (MQ=255)                                 cttcttTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGATGTTCATATTGCTACcaca                                                                                                                     >  1:298354/1‑148 (MQ=255)                                  ttcgtTATGTGGGTTGCTTCACGCGGAATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATAATGCGAACACaac                                                                                                                   >  6:293549/1‑148 (MQ=255)                                        aTGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATAGTCGGTTCTGTGGCAGGGATTGCTATGGTTGTTATGTGACTCCTGTATGTGTCGATGAAGGTAACCGCGCCTGCATTTACTGAAGTTCATATTGAGACCACAACCATTAc                                                                                                             >  2:271665/1‑149 (MQ=255)                                         tgtcGGTTGCTTCACGCGGTATGAAATCGCTGAAACTCGTCGGTTCTGTGGCAGGGCTTGCTATGTTTGTTATGGCACTCCTGTATGTGGCGATGTCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGTCCTCTAACAATAcc                                                                                                            <  8:167778/149‑1 (MQ=255)                                          gtgGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTAcc                                                                                                            <  4:229727/148‑1 (MQ=255)                                             ggTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTgg                                                                                                         >  3:272149/1‑148 (MQ=255)                                                tGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGCCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGaaa                                                                                                     >  7:145106/1‑149 (MQ=255)                                                   ttCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCGATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGGATATTGCGACCACAAACATTACCTGGGAAACGt                                                                                                  >  4:72449/1‑149 (MQ=255)                                                        gcgGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCAtt                                                                                             >  8:176680/1‑149 (MQ=255)                                                         cgGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTc                                                                                            >  6:228124/1‑149 (MQ=255)                                                          ggTATTAAATCGCTGAAAATCGTCGGTTCTGTCGCAGGGATTGATATGTTTGTTATGTCACTCCTGTATGTGTCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTcc                                                                                           <  5:244681/149‑1 (MQ=255)                                                            tATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGATACGTTCATTCCTc                                                                                         <  8:103792/149‑1 (MQ=255)                                                            tATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTc                                                                                         <  1:225398/149‑1 (MQ=255)                                                            tATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGATATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTc                                                                                         <  6:227667/149‑1 (MQ=255)                                                              tGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCat                                                                                       <  7:40017/149‑1 (MQ=255)                                                                 aaTCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGACTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCAt                                                                                              >  4:306286/1‑139 (MQ=255)                                                                 aaTCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGACTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCAt                                                                                              >  4:301134/1‑139 (MQ=255)                                                                 aaTCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGACTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCAt                                                                                              <  3:301134/139‑1 (MQ=255)                                                                 aaTCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGACTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCAt                                                                                              <  3:306286/139‑1 (MQ=255)                                                                     gCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGCCCACAAACATTACCTGGGCAACGTTCATGCCTCATATCGACt                                                                                >  4:251076/1‑149 (MQ=255)                                                                         aaaaTCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACt                                                                                                                                          >  1:191224/1‑87 (MQ=255)                                                                         aaaaTCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACt                                                                                                                                          <  2:191224/87‑1 (MQ=255)                                                                          aaaTCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTAcc                                                                           <  3:92217/149‑1 (MQ=255)                                                                           aaTCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCt                                                                          <  7:193314/149‑1 (MQ=255)                                                                                tcgGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACAt                                                                      >  8:62051/1‑148 (MQ=255)                                                                                tcgGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACAt                                                                      >  8:309195/1‑148 (MQ=255)                                                                                   gTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTc                                                                                         <  7:1584/126‑1 (MQ=255)                                                                                   gTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTc                                                                                         >  8:1584/1‑126 (MQ=255)                                                                                   gTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTAcc                                                                  <  8:32960/149‑1 (MQ=255)                                                                                         tgGCAGGGATTGCTATGTTTGTTATGGCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTa                                                              <  2:122842/147‑1 (MQ=255)                                                                                               ggATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTACCCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATg                                                      >  7:359100/1‑149 (MQ=255)                                                                                                      tATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTGACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGttt                                               >  3:167399/1‑149 (MQ=255)                                                                                                        tttttgttATGTCACTCCTGTATGTGTCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGAATATTGCGACCACAAACATTACATGTGAAACGTTCATTCCTGATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTc                                             <  5:179206/147‑1 (MQ=255)                                                                                                         gtttgttATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATTTTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTcg                                            >  8:321489/1‑149 (MQ=255)                                                                                                          tttgttATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCata                                                                                      >  5:33817/1‑106 (MQ=255)                                                                                                          tttgttATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCata                                                                                      <  6:33817/106‑1 (MQ=255)                                                                                                               tATGTCACTCCTGTATGTGGCGATGGCTGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTg                                                   <  5:120973/136‑1 (MQ=255)                                                                                                               tATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTg                                                   >  6:120973/1‑136 (MQ=255)                                                                                                                 tGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTgg                                     <  7:176680/148‑1 (MQ=255)                                                                                                                 tGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGGTGGc                                    <  3:41978/149‑1 (MQ=255)                                                                                                                  gTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCg                                   >  5:373283/1‑149 (MQ=255)                                                                                                                  gTCACTCCTGTATCTGGCGATGGCGGTAACCGCGCCTGCCATTACTGAAGTGCATATTGCGACCACAAACAATAGCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACAGTACCACTATTTCAATGCTCGATGTCGCGGTTGGCg                                   >  4:356894/1‑149 (MQ=255)                                                                                                                   tCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTAc                                                                            <  5:176893/107‑1 (MQ=255)                                                                                                                   tCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTAc                                                                            >  6:176893/1‑107 (MQ=255)                                                                                                                   tCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCg                                   <  7:62051/148‑1 (MQ=255)                                                                                                                    cACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCg                                   >  4:270969/1‑147 (MQ=255)                                                                                                                     aCTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAg                                >  2:371265/1‑149 (MQ=255)                                                                                                                      cTCCTGTATGTGGCGATGTCGGTAACAGCGCCTGAAATTACTGAAGTGCATATTGCGACAAAAAACATTACCTGGGAAAGGTGCATTCAGAAATAGGAATTAACAGAAATTCCCAATATTTCAATAATGGTTTTCGAGGTTGGCGGAgg                               >  6:89641/1‑148 (MQ=255)                                                                                                                       tCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAg                                >  5:128299/1‑147 (MQ=255)                                                                                                                          tGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCa                              >  8:222724/1‑146 (MQ=255)                                                                                                                            tATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTgg                                     >  8:269722/1‑137 (MQ=255)                                                                                                                            tATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTgg                                     <  7:269722/137‑1 (MQ=255)                                                                                                                            tATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGagaa                         <  3:270969/149‑1 (MQ=255)                                                                                                                               gtgGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTAc                                                                            <  6:264598/95‑1 (MQ=255)                                                                                                                               gtgGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTAc                                                                            >  5:264598/1‑95 (MQ=255)                                                                                                                                  gCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGAGAAGATTtc                   <  5:228124/149‑1 (MQ=255)                                                                                                                                    gATTGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACAACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGAGAAGATTtctc                 <  2:245943/149‑1 (MQ=255)                                                                                                                                    gATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGAGAAGATTtctc                 <  7:153048/149‑1 (MQ=255)                                                                                                                                      tGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGAGAAGATTTCTCCt               <  4:272149/149‑1 (MQ=255)                                                                                                                                           gTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGAGAAGATTTCTCCTTACGt          <  6:373283/149‑1 (MQ=255)                                                                                                                                                 gcgcCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGAGAAGATTTCTCCTTACGTTAATc     >  8:239472/1‑148 (MQ=255)                                                                                                                                                   gcCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGAGAAGATTTCTCCTTACGTTAATCaaa  <  8:145106/149‑1 (MQ=255)                                                                                                                                                   gcCTGCAATAACTGAAGTGCATATTGCGACCACAAACATTACCTCGGAAACGTTCATACCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGAGAAGATTTCTCCTTACGTTAATCaaa  <  2:87610/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCGTTACAGGGGTTAACGCTGGTGCTGTTTATCTTCTTTATGTGGGTTGCTTCACGCGGTATGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCACTCCTGTATGTGGCGATGGCGGTAACCGCGCCTGCAATTACTGAAGTGCATATTGCGACCACAAACATTACCTGGGAAACGTTCATTCCTCATATCGACTTTACCTACATTACCACTATTTCAATGCTGGTTTTCGCGGTTGGCGGAGCAGAGAAGATTTCTCCTTACGTTAATCAAA  >  NC_000913/947616‑947910

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