Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,901,659	A→G	9.7%	T375A (ACC→GCC)	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,659	0	A	G	9.7%	162.5 / 7.5	62	T375A (ACC→GCC)	PP_4309	NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base G (4/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

GCCATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACG  >  NC_002947/4901514‑4901806                                                                                                                                                   |                                                                                                                                                    gCCATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCAccc                                                                                                                                                  <  4:268349/149‑1 (MQ=255) gCCATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCAccc                                                                                                                                                  >  5:293291/1‑149 (MQ=255)   cATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTg                                                                                                                                                >  5:99428/1‑149 (MQ=255)            tGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTg                                                                                                                                       <  7:13352/149‑1 (MQ=255)                 ggCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATc                                                                                                                                    <  3:263808/147‑1 (MQ=255)                 ggCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCt                                                                                                                                  >  4:41302/1‑149 (MQ=255)                     gCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGccc                                                                                                                              >  2:297264/1‑149 (MQ=255)                      ccTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGcc                                                                                                                               <  6:99428/147‑1 (MQ=255)                        tGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTg                                                                                                                            >  6:158025/1‑148 (MQ=255)                          gCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGaa                                                                                                                         <  3:148432/149‑1 (MQ=255)                          gCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TGTCGTACTCAGCAACCTCGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGaa                                                                                                                         >  1:221115/1‑149 (MQ=255)                           ccTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAAc                                                                                                                        >  7:87176/1‑149 (MQ=255)                               ccTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGt                                                                                                                    <  8:87176/149‑1 (MQ=255)                                 tGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGCATTTCGTACTCAGCAACCTGGCGCCGCGCCACCTCAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGt                                                                                                                    >  1:76914/1‑148 (MQ=255)                                   ctCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTac                                                                                                                  <  6:225076/146‑1 (MQ=255)                                   aTCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTac                                                                                                                  <  3:181706/147‑1 (MQ=255)                                   aTCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TGTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTacaa                                                                                                                <  2:221115/149‑1 (MQ=255)                                    tCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCAGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTacaa                                                                                                                >  1:172747/1‑148 (MQ=255)                                       tCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGc                                                                                                            >  7:146461/1‑149 (MQ=255)                                           cGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGc                                                                                                        <  3:41302/149‑1 (MQ=255)                                                     aaCATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCgt                                                                                              <  2:76914/149‑1 (MQ=255)                                                     aaCATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCT‑GCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCgtg                                                                                             <  5:89852/149‑1 (MQ=255)                                                              gCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTcc                                                                                      <  1:310248/148‑1 (MQ=255)                                                                cAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACATGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACACCAGCCCGCGGGTGATCCTGTACTTCCTGt                                                                                   >  8:206513/1‑149 (MQ=255)                                                                 caCTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTa                                                                                                                   <  4:71957/115‑1 (MQ=255)                                                                      cgtcgCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGc                                                                              <  7:127027/148‑1 (MQ=255)                                                                           ccccGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACACCAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGgc                                                                        >  3:287745/1‑149 (MQ=255)                                                                            cccGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCGGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGgcg                                                                       <  6:293291/149‑1 (MQ=255)                                                                               ggC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGccc                                                                    <  6:172832/149‑1 (MQ=255)                                                                                gC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCt                                                                   >  1:203979/1‑149 (MQ=255)                                                                                 c‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTg                                                                  >  7:212239/1‑149 (MQ=255)                                                                                        tACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGc                                                            >  8:70675/1‑149 (MQ=255)                                                                                           tCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCg                                                         >  4:310404/1‑149 (MQ=255)                                                                                           tCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCg                                                         >  7:8560/1‑149 (MQ=255)                                                                                             aGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACATCAGCCCGCTGGTGAt                                                                                                 <  5:238568/107‑1 (MQ=255)                                                                                             aGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCt                                                       >  6:52915/1‑149 (MQ=255)                                                                                             aGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTCGCGCTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGAt                                                                                                 >  6:238568/1‑107 (MQ=255)                                                                                                  ccTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACg                                                  <  5:158025/149‑1 (MQ=255)                                                                                                    tGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGc                                                >  2:278143/1‑149 (MQ=255)                                                                                                      gcgccgcgccACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGt                                              <  2:172747/149‑1 (MQ=255)                                                                                                                  cTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCgg                                  <  8:146461/149‑2 (MQ=255)                                                                                                                  cTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGa                                  >  5:205839/1‑149 (MQ=255)                                                                                                                   tGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGac                                 >  2:205325/1‑149 (MQ=255)                                                                                                                       cTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCt                                                       <  6:226642/123‑1 (MQ=255)                                                                                                                       cTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCt                                                       >  5:226642/1‑123 (MQ=255)                                                                                                                          ccgccgTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTgg                           <  2:272082/148‑1 (MQ=255)                                                                                                                                tGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTtgctgc                    <  5:27782/149‑1 (MQ=255)                                                                                                                                 gCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTgcg                   >  3:136045/1‑149 (MQ=255)                                                                                                                                     ggCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCaaa               <  6:205839/149‑1 (MQ=255)                                                                                                                                      gCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAg              <  2:203979/149‑1 (MQ=255)                                                                                                                                      caggATCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAg              >  7:72799/4‑149 (MQ=255)                                                                                                                                      caggATCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAg              >  7:22524/4‑149 (MQ=255)                                                                                                                                       cTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAgg             <  4:136045/149‑1 (MQ=255)                                                                                                                                         gATCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCt           >  1:264212/1‑149 (MQ=255)                                                                                                                                          aTCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGTCt           <  3:310404/148‑1 (MQ=255)                                                                                                                                          aTCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGc            <  7:58275/147‑1 (MQ=255)                                                                                                                                          aTCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTtgctgc                    <  4:181913/139‑1 (MQ=255)                                                                                                                                          aTCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTtgctgc                    >  3:181913/1‑139 (MQ=255)                                                                                                                                           tCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGc         <  2:264212/149‑1 (MQ=255)                                                                                                                                               cgcCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCt           <  4:267226/143‑1 (MQ=255)                                                                                                                                               cgcCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCt           >  3:267226/1‑143 (MQ=255)                                                                                                                                                  cACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACg  >  6:83816/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GCCATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACG  >  NC_002947/4901514‑4901806

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