Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,632,085	G→T	5.2%	E327* (GAG→TAG) ‡	PP_5590 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,632,085	0	G	T	5.2%	252.3 / 2.9	76	E327* (GAG→TAG) ‡	PP_5590	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base T (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGG  >  NC_002947/3631939‑3632229                                                                                                                                                   |                                                                                                                                                 cGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAg                                                                                                                                                >  1:1205/1‑149 (MQ=255) cGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAg                                                                                                                                                <  5:385541/149‑1 (MQ=255)        ggCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGa                                                                                                                                                 <  5:190249/141‑1 (MQ=255)        ggCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGa                                                                                                                                                 >  6:190249/1‑141 (MQ=255)          ccACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGttgctg                                                                                                                                       >  2:196178/1‑149 (MQ=255)           cACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTgg                                                                                                                                      <  1:196178/149‑1 (MQ=255)             ccgccgGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTGGcg                                                                                                                                    <  6:314056/149‑2 (MQ=255)                     gTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCa                                                                                                                            >  2:248246/1‑149 (MQ=255)                        aTGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCt                                                                                                                         <  6:224031/149‑1 (MQ=255)                              ttCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAAccc                                                                                                                   >  3:11757/1‑149 (MQ=255)                              ttCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAAccc                                                                                                                   >  7:219191/1‑149 (MQ=255)                                  aTGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGctgctg                                                                                                                                       <  6:404885/125‑1 (MQ=255)                                  aTGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGctgctg                                                                                                                                       >  5:404885/1‑125 (MQ=255)                                  aTGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGAc                                                                                                               <  7:236310/149‑1 (MQ=255)                                  aTGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGAc                                                                                                               >  2:69268/1‑149 (MQ=255)                                          tACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATc                                                                                                         <  5:268766/147‑1 (MQ=255)                                             gtgtTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTg                                                                                                    >  4:154709/1‑149 (MQ=255)                                             gtgtTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTg                                                                                                    >  5:347759/1‑149 (MQ=255)                                               gtTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTgg                                                                                                   >  2:289799/1‑148 (MQ=255)                                                  cAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAg                                                                                                <  1:248246/148‑1 (MQ=255)                                                  cAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGc                                                                                               >  8:264741/1‑149 (MQ=255)                                                   aaCGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCt                                                                                                     <  7:236398/142‑1 (MQ=255)                                                   aaCGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCt                                                                                                     >  8:236398/1‑142 (MQ=255)                                                   aaCGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCt                                                                                              >  7:220061/1‑149 (MQ=255)                                                    aCGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTg                                                                                             <  8:65529/149‑1 (MQ=255)                                                     cGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTg                                                                                             >  1:279480/1‑148 (MQ=255)                                                             tCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtg                                                                                    >  3:316802/1‑149 (MQ=255)                                                                  cGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTgtgggt                                                                               <  1:272370/149‑1 (MQ=255)                                                                  cGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTgtgggt                                                                               >  2:272370/1‑149 (MQ=255)                                                                   gCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTgtgggtg                                                                              <  6:347759/149‑1 (MQ=255)                                                                        gAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGa                                                                         <  1:289799/149‑1 (MQ=255)                                                                         aGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAg                                                                        <  5:252927/149‑1 (MQ=255)                                                                          gTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAgg                                                                       >  8:134381/1‑149 (MQ=255)                                                                               aGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACg                                                                  <  1:22209/149‑1 (MQ=255)                                                                                gCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGc                                                                 >  5:298113/1‑149 (MQ=255)                                                                                  tggtggGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTTGTGGAGGCGACGCt                                                                >  2:186980/1‑148 (MQ=255)                                                                                    gtggGCAAGAAGGAGATGTACATACCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTgaa                                                             <  4:11757/149‑1 (MQ=255)                                                                                      gggCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGATGCGCTGCTGGCCCCCCATCACCTGAACCCCGAc                                                                                                               <  7:260022/97‑1 (MQ=255)                                                                                         cAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATg                                                        >  3:221871/1‑149 (MQ=255)                                                                                           agaagGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGc                                                      <  8:219191/149‑1 (MQ=255)                                                                                               ggAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTTGTGGAGGCGACGCTGAAGAATGGCCAgc                                                  <  1:186980/149‑1 (MQ=255)                                                                                               ggAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAgc                                                  <  2:1205/149‑1 (MQ=255)                                                                                                       aCATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATg                                                        <  8:118236/135‑1 (MQ=255)                                                                                                       aCATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATg                                                        >  7:118236/1‑135 (MQ=255)                                                                                                            cccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGcc                                     <  7:264741/149‑1 (MQ=255)                                                                                                             ccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCCCCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCcg                                    >  4:144116/1‑149 (MQ=255)                                                                                                             ccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCcg                                    >  4:229971/1‑149 (MQ=255)                                                                                                              cTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgt                                                                                       <  8:268047/97‑1 (MQ=255)                                                                                                              cTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgt                                                                                       >  7:268047/1‑97 (MQ=255)                                                                                                                 cAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGt                                 <  4:221871/148‑1 (MQ=255)                                                                                                                 cAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGt                                 <  2:279480/148‑1 (MQ=255)                                                                                                                   aCGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGc                              <  8:220061/149‑1 (MQ=255)                                                                                                                    cGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgt                                                                                       >  4:76575/1‑91 (MQ=255)                                                                                                                     gCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCAc                            >  2:394913/1‑149 (MQ=255)                                                                                                                     gCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCAc                            <  6:298113/149‑1 (MQ=255)                                                                                                                       caacaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGGGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACct                          >  8:175735/1‑149 (MQ=255)                                                                                                                         acaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAg                                                                        >  7:156050/1‑101 (MQ=255)                                                                                                                         acaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAg                                                                        <  8:156050/101‑1 (MQ=255)                                                                                                                          caaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACctac                        >  5:255009/1‑148 (MQ=255)                                                                                                                           aaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACctacta                      >  1:190553/1‑149 (MQ=255)                                                                                                                            aGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCAtgtg                                          <  3:305244/128‑1 (MQ=255)                                                                                                                            aGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTCGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCAtgtg                                          >  4:305244/1‑128 (MQ=255)                                                                                                                              gtgtTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGtggtgg                                                                          >  3:120763/1‑94 (MQ=255)                                                                                                                              gtgtTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGtggtgg                                                                          <  4:120763/94‑1 (MQ=255)                                                                                                                              gtgtTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTgaa                                                             >  6:235264/1‑107 (MQ=255)                                                                                                                              gtgtTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTgaa                                                             <  5:235264/107‑1 (MQ=255)                                                                                                                                 ttCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTcga                >  4:381273/1‑149 (MQ=255)                                                                                                                                  tCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTcga                >  8:58971/1‑148 (MQ=255)                                                                                                                                   cACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTcgac               >  1:354279/1‑148 (MQ=255)                                                                                                                                     ccGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGaa            <  6:255009/149‑1 (MQ=255)                                                                                                                                     ccGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGaa            <  1:394913/149‑1 (MQ=255)                                                                                                                                      cGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAg                                                                        >  7:364880/1‑88 (MQ=255)                                                                                                                                      cGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAg                                                                        <  8:364880/88‑1 (MQ=255)                                                                                                                                        ccAGCAGCGCCTCGGCGTTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGa          >  7:314206/1‑148 (MQ=255)                                                                                                                                          aGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGacac       <  3:229971/149‑1 (MQ=255)                                                                                                                                               cGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTggg  <  2:354279/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGG  >  NC_002947/3631939‑3632229

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