Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,901,663	T→G	9.2%	L376R (CTG→CGG) ‡	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,663	0	T	G	9.2%	180.7 / 7.4	63	L376R (CTG→CGG) ‡	PP_4309	NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base T (26/31);  new base G (4/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGC  >  NC_002947/4901516‑4901809                                                                                                                                                     |                                                                                                                                                   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                       >  3:181913/1‑139 (MQ=255)                                                                                                                                        aTCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTtgctgc                       <  4:181913/139‑1 (MQ=255)                                                                                                                                         tCAGCACCGCCAGGGTGGCGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGc            <  2:264212/149‑1 (MQ=255)                                                                                                                                             cgcCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCt              >  3:267226/1‑143 (MQ=255)                                                                                                                                             cgcCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCt              <  4:267226/143‑1 (MQ=255)                                                                                                                                                cACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACg     >  6:83816/1‑149 (MQ=255)                                                                                                                                                   ccTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGc  <  5:83816/149‑1 (MQ=255)                                                                                                                                                    cTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTg                               <  7:217100/119‑1 (MQ=255)                                                                                                                                                    cTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTg                               >  8:217100/1‑119 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CATTCCTGGTGCTCGGCTGCCTGGCCTTCCTGATCGTCACCGTAGCGGTGAACATCATGGCCAACTTCGTCGCCCCGGC‑TTTCGTACTCAGCAACCTGGCGCCGCGCCACCTGAACTTCCGCCGTGCCGGGCTGATCAGCGCCACCCTGGCGGTGCTGATCCTGCCCTGGAACCTGTACAACAGCCCGCTGGTGATCGTGTACTTCCTGTCGGGCCTGGGCGCCCTGCTAGGCCCGCTGTACGGCGTGATCATGTCCGACTACTGGTTGCTGCGCAAAGGCTGCATCAACGTGC  >  NC_002947/4901516‑4901809

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