Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,901,656	G→A	6.9%	A374T (GCC→ACC)	PP_4309 →	NCS1 family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,901,656	0	G	A	6.9%	215.2 / 4.1	72	A374T (GCC→ACC)	PP_4309	NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base A (3/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

GCCATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCAT  >  NC_002947/4901514‑4901801                                                                                                                                               |                                                                                                                                                   gCCATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCAccc                                                                                                                                              >  1:190059/1‑149 (MQ=255)        cTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGTACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGt                                                                                                                                       <  3:358620/149‑1 (MQ=255)         tGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTg                                                                                                                                      >  5:294197/1‑149 (MQ=255)           gTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTg                                                                                                                                      <  6:294197/147‑1 (MQ=255)               tCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATc                                                                                                                                >  6:355261/1‑149 (MQ=255)                 ggCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATc                                                                                                                                <  5:214801/147‑1 (MQ=255)                 ggCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCACCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATc                                                                                                                                <  5:4377/147‑1 (MQ=255)                           ccTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAAc                                                                                                                    <  8:258069/149‑1 (MQ=255)                               ccTGATCGTCACCGTAGCGGTGAACATCCAGGCCAACTTCGTCGCCCCGGCGTTCGTCCTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTg                                                                                                                 <  8:347515/148‑1 (MQ=255)                                cTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTa                                                                                                               <  6:205832/149‑1 (MQ=255)                                cTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTa                                                                                                               >  5:205832/1‑149 (MQ=255)                                   aTCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAAc                                                                                                                    <  5:250034/141‑1 (MQ=255)                                   aTCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAAc                                                                                                                    >  6:250034/1‑141 (MQ=255)                                   aTCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTac                                                                                                              <  2:190059/147‑1 (MQ=255)                                     cGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCt                                                                                                                                    <  4:14053/123‑1 (MQ=255)                                     cGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCt                                                                                                                                    >  3:14053/1‑123 (MQ=255)                                       tCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTg                                                                                                                                      <  3:354981/119‑1 (MQ=255)                                       tCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTg                                                                                                                                      >  4:354981/1‑119 (MQ=255)                                       tCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGcc                                                                                                                           >  6:1581/1‑130 (MQ=255)                                       tCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGcc                                                                                                                           <  5:1581/130‑1 (MQ=255)                                       tCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAAc                                                                                                                    <  6:249658/137‑1 (MQ=255)                                       tCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAAc                                                                                                                    >  5:249658/1‑137 (MQ=255)                                         aCCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGcc                                                                                                       <  4:112042/148‑1 (MQ=255)                                         aCCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGcc                                                                                                       <  6:168621/148‑1 (MQ=255)                                             tAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTac                                                                                                              <  7:131931/137‑1 (MQ=255)                                             tAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTac                                                                                                              >  8:131931/1‑137 (MQ=255)                                                        atcatGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑a                                                                                      >  4:180774/1‑149 (MQ=255)                                                        atcatGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑a                                                                                      >  2:313150/1‑149 (MQ=255)                                                        atcatGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑a                                                                                      <  5:355261/149‑1 (MQ=255)                                                          catGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACt                                                                                    >  1:146797/1‑149 (MQ=255)                                                            tGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTc                                                                                  <  2:146797/149‑1 (MQ=255)                                                             ggCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTcc                                                                                 >  8:356768/1‑149 (MQ=255)                                                             ggCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTc                                                                                  <  8:123562/148‑1 (MQ=255)                                                             ggCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTc                                                                                  >  7:123562/1‑148 (MQ=255)                                                                 aaCTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑Ac                                                                                     >  6:206825/1‑141 (MQ=255)                                                                 aaCTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑Ac                                                                                     <  5:206825/141‑1 (MQ=255)                                                                  aCTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGTTGATCGTGT‑ACTTCCTGTCg                                                                            <  3:180774/149‑1 (MQ=255)                                                                        tcgCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTg                                                                      >  4:107368/1‑149 (MQ=255)                                                                                 cTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTg                                                             >  8:181592/1‑149 (MQ=255)                                                                                    tCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGc                                                            <  1:51890/147‑1 (MQ=255)                                                                                       tACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGc                                                       >  4:84686/1‑149 (MQ=255)                                                                                           cAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCACCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGc                                                   >  3:71367/1‑149 (MQ=255)                                                                                            aGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGc                                                   <  1:191937/148‑1 (MQ=255)                                                                                            aGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCt                                                  >  1:110693/1‑149 (MQ=255)                                                                                            aGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCt                                                  <  7:164664/149‑1 (MQ=255)                                                                                            aGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGACTGGGCGCCCTGCTAGGCCCGCt                                                  >  8:164664/1‑149 (MQ=255)                                                                                             gCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTg                                                 >  8:260124/1‑149 (MQ=255)                                                                                                    ggCGCCGCGCCACCTGAACTTCAGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTCT‑ACATCCTGTCGGGCCTGGGCGCCCTGATAGGCCCGCTGTACGGct                                          >  2:84176/1‑148 (MQ=255)                                                                                                       tccgcgccACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACTGCCCGCTGGTGATCGTGTCACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCTCTGTACGGCGTg                                        <  1:84176/148‑1 (MQ=255)                                                                                                               aCCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTcc                               <  1:313150/149‑1 (MQ=255)                                                                                                                    aaCTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGacta                          >  6:174645/1‑149 (MQ=255)                                                                                                                    aaCTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGact                           <  1:20437/148‑1 (MQ=255)                                                                                                                    aaCTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGact                           <  1:26293/148‑1 (MQ=255)                                                                                                                    aaCTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGac                            <  3:84686/147‑1 (MQ=255)                                                                                                                      cTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGactact                        <  7:181592/149‑1 (MQ=255)                                                                                                                      cTTCCACCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGactact                        >  2:264645/1‑149 (MQ=255)                                                                                                                       ttCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTa                                               >  8:78264/1‑125 (MQ=255)                                                                                                                       ttCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTa                                               <  7:78264/125‑1 (MQ=255)                                                                                                                       ttCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTg                       >  6:280715/1‑149 (MQ=255)                                                                                                                                      cTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAgg        <  2:110693/149‑1 (MQ=255)                                                                                                                                         aTCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGc       <  7:356768/147‑1 (MQ=255)                                                                                                                                         aTCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTg     >  5:267933/1‑149 (MQ=255)                                                                                                                                         aTCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTgg                      >  1:351792/1‑132 (MQ=255)                                                                                                                                         aTCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTgg                      <  2:351792/132‑1 (MQ=255)                                                                                                                                         aTCAGCACCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTg     <  4:71367/149‑1 (MQ=255)                                                                                                                                          tCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTc                                                                                  <  1:366579/71‑1 (MQ=255)                                                                                                                                          tCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTc                                                                                  >  2:366579/1‑71 (MQ=255)                                                                                                                                          tCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTg     >  6:361970/1‑148 (MQ=255)                                                                                                                                          tCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTg     >  1:291438/1‑148 (MQ=255)                                                                                                                                          tCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGc    <  7:260124/149‑1 (MQ=255)                                                                                                                                          tCAGCACCGCCAGGGTGGCGCTGATTCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTg     >  8:262693/1‑148 (MQ=255)                                                                                                                                            aGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCTGACTACTGGTTGCTGCGCAAAGGCTGCAt  >  7:323049/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                   GCCATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCTTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGT‑ACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCAT  >  NC_002947/4901514‑4901801

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