Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,901,656	G→A	10.1%	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	10.1%	160.0 / 6.6	59	A374T (GCC→ACC)	PP_4309	NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (4/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

GCCATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGC  >  NC_002947/4901514‑4901799                                                                                                                                                |                                                                                                                                                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:22524/4‑149 (MQ=255)                                                                                                                                      caggATCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAg       >  7:72799/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    >  3:267226/1‑143 (MQ=255)                                                                                                                                               cgcCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCt    <  4:267226/143‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                GCCATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGC  >  NC_002947/4901514‑4901799

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