Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,402,301:1	+A	6.9%	coding (66/330 nt)	PP_2103 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,402,301	1	.	A	6.9%	280.4 / 16.0	72	coding (66/330 nt)	PP_2103	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (40/27);  new base A (3/2);  total (43/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

CACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAG  >  NC_002947/2402165‑2402449                                                                                                                                          |                                                                                                                                                     cACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCagcagc                                                                                                                                          <  4:253403/149‑1 (MQ=255)      aTGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCtt                                                                                                                                     >  3:309230/1‑149 (MQ=255)         ccccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgc                                                                                                                                   <  4:309230/148‑1 (MQ=255)         ccccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgc                                                                                                                                   <  2:34348/148‑1 (MQ=255)         ccccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgc                                                                                                                                   <  6:252458/148‑1 (MQ=255)          cccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcgc                                                                                                                                 <  1:202610/149‑1 (MQ=255)          cccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcg                                                                                                                                  >  5:357694/1‑148 (MQ=255)          cccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcg                                                                                                                                  >  7:57779/1‑148 (MQ=255)           ccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcg                                                                                                                                  >  8:308439/1‑147 (MQ=255)           ccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCa                                                                                                                                >  3:40271/1‑149 (MQ=255)                 tGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGCAAGGCGCTGCTGCGCC‑TGCGCatggaaa                                                                                                                          <  3:342047/149‑7 (MQ=255)                 tGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgc                                                                                                                          >  7:42/1‑149 (MQ=255)                     gcATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCTCCTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCgg                                                                                                                      <  1:120511/149‑1 (MQ=255)                     gcATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGATCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCgg                                                                                                                      >  1:198950/1‑149 (MQ=255)                       aTGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGCAAGGCGCTGCTGCGCCttcgcatg                                                                                                                               <  3:158202/139‑7 (MQ=255)                       aTGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGCAAGGCGCTGCTGCGCC‑TGCGCatg                                                                                                                              >  4:158202/1‑137 (MQ=255)                           ttGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGAGCGGACTCATGACGCATTTCCTGGCGGTGCATTTCCATGAGC‑AGGCGCAGCAGCGCCTTGAGCAGTTCGCGCGGGTTAcg                                                                                                                >  5:127180/1‑149 (MQ=255)                              aTACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCgt                                                                                                             >  3:323444/1‑149 (MQ=255)                                 cccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGtt                                                                                                          >  4:296928/1‑149 (MQ=255)                                  ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGttt                                                                                                         >  4:335629/1‑149 (MQ=255)                                  ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGttt                                                                                                         >  2:364186/1‑149 (MQ=255)                                  ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGttt                                                                                                         <  3:335629/149‑1 (MQ=255)                                  ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAACCAACTGACCGGACTCATGAAGGATTTCCTCGCGGTGCATTTCCATGCGC‑AGGTGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGttt                                                                                                         >  5:318612/1‑149 (MQ=255)                                         ttCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTgg                                                                                                       >  1:135059/1‑144 (MQ=255)                                         ttCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTgg                                                                                                       <  2:135059/144‑1 (MQ=255)                                            aTGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCt                                                                                                <  8:42/148‑1 (MQ=255)                                            aTGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTg                                                                                               >  6:340583/1‑149 (MQ=255)                                             tGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTgg                                                                                              <  7:308439/149‑1 (MQ=255)                                               tAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGc                                                                                             <  1:364186/148‑1 (MQ=255)                                               tAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCa                                                                                            <  3:296928/149‑1 (MQ=255)                                                    aGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCt                                                                                        >  3:163478/1‑148 (MQ=255)                                                           cATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCAGGGCaa                                                                                >  1:163453/1‑149 (MQ=255)                                                           cATACCGCGCACGCGGCGCAGCGGATCAAGCAACTGCCCGGACTCATGACGGATTACCTGGCCGTGCATTTCCATGCGC‑AGGCGCACCAGCCCCTTGCGCAGTTCCCGCGGGTTACGCGTGTTTgg                                                                                                       >  5:33954/1‑126 (MQ=255)                                                           cATACCGCCGACGGGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGGATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTgg                                                                                                       <  6:33954/126‑1 (MQ=255)                                                                cgcgCACGCGGCGCAGCGGTTCAAGCAACTGACCGGCCTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGgc                                                                           >  6:22993/1‑149 (MQ=255)                                                                        cggcgCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTAc                                                                                                                 <  7:321576/103‑1 (MQ=255)                                                                        cggcgCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTAc                                                                                                                 >  8:321576/1‑103 (MQ=255)                                                                              aGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTg                                                             >  7:291590/1‑149 (MQ=255)                                                                               gCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTgg                                                            >  2:75009/1‑149 (MQ=255)                                                                               gCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTgg                                                            >  5:70551/1‑149 (MQ=255)                                                                                  gTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCaa                                                         <  8:57779/149‑1 (MQ=255)                                                                                   ttCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAg                                                        <  2:76795/149‑1 (MQ=255)                                                                                      aaGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtc                                                     >  3:256491/1‑149 (MQ=255)                                                                                       aGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTg                                                             <  5:364308/140‑1 (MQ=255)                                                                                       aGCAACTGACCGGACTCATGACGGATTTCATGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGACTTGCGCAGTTCGCGCGCGTTACGCGTGTTTGTCAGTTCTGGCAGGCTCATGGAAACAGGCGCTCCCAGTCCTTg                                                             >  6:364308/1‑140 (MQ=255)                                                                                           aCTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTcc                                                                      >  2:321490/1‑127 (MQ=255)                                                                                           aCTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTcc                                                                      <  1:321490/127‑1 (MQ=255)                                                                                               aCCGGACTCATGACGCATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTAGCGCGGGTTACCCGAGTTTGGCAGTTCTGGCAGGCTCATGGCAAAACCCGCTCCCTGTCATTGGCAAGTTCGTCGAggct                                            >  6:319916/1‑147 (MQ=255)                                                                                                  ggACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGgc                                         >  4:240612/1‑149 (MQ=255)                                                                                                  ggACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGgc                                         >  7:2783/1‑149 (MQ=255)                                                                                                      tCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGcc                                      >  3:137432/1‑148 (MQ=255)                                                                                                        aTGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGaa                                   <  6:70551/149‑1 (MQ=255)                                                                                                         tGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAg                                  >  8:178523/1‑149 (MQ=255)                                                                                                         tGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAg                                  >  2:194828/1‑149 (MQ=255)                                                                                                           aCGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAggg                                <  8:291590/149‑1 (MQ=255)                                                                                                           aCGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAggg                                >  6:130256/1‑149 (MQ=255)                                                                                                           aCGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGGCGAGGGTGGCGCCGAAgg                                 <  5:22993/148‑1 (MQ=255)                                                                                                            cGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAggg                                <  4:40271/148‑1 (MQ=255)                                                                                                              gATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCga                             <  1:194828/149‑1 (MQ=255)                                                                                                              gATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCga                             <  3:240612/149‑1 (MQ=255)                                                                                                               aTTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgat                            <  4:256491/149‑1 (MQ=255)                                                                                                                tttCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGTAAGTTCGTCGAGGGTGGCGCCGAAGGGCgat                            >  2:242574/1‑148 (MQ=255)                                                                                                                tttCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgat                            >  5:344702/1‑148 (MQ=255)                                                                                                                    cTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATtc                       >  5:251376/1‑149 (MQ=255)                                                                                                                       gCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATtcgtc                    >  4:140927/1‑149 (MQ=255)                                                                                                                       gCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATtcg                      <  5:340583/147‑1 (MQ=255)                                                                                                                         ggTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATtcgtcg                   >  4:239925/1‑148 (MQ=255)                                                                                                                                  ccATGCCC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACAccgccg         >  3:313292/1‑149 (MQ=255)                                                                                                                                   cATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACCCcgccgc        >  2:225919/1‑149 (MQ=255)                                                                                                                                    aTGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCgccgcc       <  6:127180/149‑1 (MQ=255)                                                                                                                                     tGCGCAAGGCGCTGCTGCGCC‑TGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCt      >  7:350771/1‑149 (MQ=255)                                                                                                                                     tGCGCAAGGCGCTGCTGCGCC‑TGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCt      >  2:340170/1‑149 (MQ=255)                                                                                                                                       cgc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGTAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTc    <  1:242574/149‑1 (MQ=255)                                                                                                                                       cgc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTc    <  5:130256/149‑1 (MQ=255)                                                                                                                                        gc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCtt     >  5:141724/1‑147 (MQ=255)                                                                                                                                         c‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAg  <  4:313292/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                     CACCGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAG  >  NC_002947/2402165‑2402449

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