Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	160,513	A→C	5.6%	F65V (TTC→GTC)	PP_0150 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	160,513	0	A	C	5.6%	226.3 / 2.3	70	F65V (TTC→GTC)	PP_0150	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/33);  new base C (2/2);  total (35/35)
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.85e-01

AGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGC  >  NC_002947/160366‑160660                                                                                                                                                      |                                                                                                                                                    agcCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAg                                                                                                                                                    <  3:18267/149‑1 (MQ=255)  gcCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCAGCGTTGCTGGCGACGAAGc                                                                                                                                                   >  4:190061/1‑149 (MQ=255)     aGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCTCCGTTGCTGGCGACGAAGCTTg                                                                                                                                                <  2:224406/149‑1 (MQ=255)      ggCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGc                                                                                                                                               <  5:63361/149‑1 (MQ=255)      ggCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGc                                                                                                                                               <  3:293033/149‑1 (MQ=255)      ggCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGc                                                                                                                                               <  5:79177/149‑1 (MQ=255)             acgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcg                                                                                                                                        >  6:333438/1‑149 (MQ=255)                 gcTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                    >  8:79338/1‑149 (MQ=255)                     cacaGTGCGAGTATGGCTTCAGCCAATACCTTGTCGCTGGCGGCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGACCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTg                                                                                                                                   <  1:253989/148‑1 (MQ=255)                     cacaGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTgg                                                                                                                                  <  3:40061/147‑1 (MQ=255)                         gTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTcg                                                                                                                            >  5:181449/1‑149 (MQ=255)                         gTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTcg                                                                                                                            >  4:145804/1‑149 (MQ=255)                          tGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCa                                                                                                                                              >  5:57466/1‑130 (MQ=255)                          tGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCa                                                                                                                                              <  6:57466/130‑1 (MQ=255)                          tGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTcgc                                                                                                                           <  5:64059/149‑1 (MQ=255)                              aGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGcac                                                                                                                       <  6:135844/149‑1 (MQ=255)                                    gCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGAGCTTCGCGCACCACttt                                                                                                                 >  6:201238/1‑149 (MQ=255)                                     cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGGTGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                >  7:146187/1‑149 (MQ=255)                                      ttCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                >  3:289719/1‑148 (MQ=255)                                            cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                         <  2:256316/149‑1 (MQ=255)                                            cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                         <  6:181449/149‑1 (MQ=255)                                             aaTACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGGTGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgcg                                                                                                        <  8:146187/149‑1 (MQ=255)                                             aaTACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgcg                                                                                                        <  3:145804/149‑1 (MQ=255)                                                   tGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATc                                                                                                    >  6:77954/1‑147 (MQ=255)                                                     gTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAg                                                                                                 <  5:333438/148‑1 (MQ=255)                                                          tGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg                                                                                           >  6:161854/1‑149 (MQ=255)                                                               ccgcGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTc                                                                                      <  5:161854/149‑1 (MQ=255)                                                               ccgcGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTc                                                                                      <  6:108197/149‑1 (MQ=255)                                                                cgcgcgCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgc                                                                                                                                    <  1:16269/102‑2 (MQ=255)                                                                cgcgcgCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgc                                                                                                                                    >  2:16269/1‑101 (MQ=255)                                                                 gcgcgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCt                                                                                     >  5:165062/1‑148 (MQ=255)                                                                   gcgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAg                                                                                  <  7:25765/149‑1 (MQ=255)                                                                     gcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtc                                                                                                                                      >  5:19278/1‑95 (MQ=255)                                                                     gcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtc                                                                                                                                      <  6:19278/95‑1 (MQ=255)                                                                     gcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATAGTGGTGGTGGTGTCTGAGGt                                                                                >  8:256140/1‑149 (MQ=255)                                                                      cTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTg                                                                               <  4:289719/149‑1 (MQ=255)                                                                          cGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAg                                                                            <  4:134901/148‑1 (MQ=255)                                                                          cGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGt                                                                           <  3:190061/149‑1 (MQ=255)                                                                                tcgctcgcGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggt                                                                                               >  6:87567/1‑123 (MQ=255)                                                                                tcgctcgcGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggt                                                                                               <  5:87567/123‑1 (MQ=255)                                                                                  gctcgcGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGaa                                                                   >  7:78041/1‑149 (MQ=255)                                                                                      gcgcAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTg                                                                                                                                   <  7:86616/81‑1 (MQ=255)                                                                                      gcgcAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTg                                                                                                                                   >  8:86616/1‑81 (MQ=255)                                                                                            ggTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGa                                                                                   <  7:99333/123‑1 (MQ=255)                                                                                            ggTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGa                                                                                   >  8:99333/1‑123 (MQ=255)                                                                                                ttGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacg                                                     >  2:294077/1‑149 (MQ=255)                                                                                                  gAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgca                                                   >  3:71954/1‑149 (MQ=255)                                                                                                   aaGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGGGAAGTCGATCGAACGGCCAAAGgcacgcac                                                  >  8:126955/1‑149 (MQ=255)                                                                                                     ggCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACga                                                >  6:81456/1‑149 (MQ=255)                                                                                                       ccg‑‑ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgatg                                              <  6:165062/149‑1 (MQ=255)                                                                                                        cg‑‑ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACga                                                >  6:208034/1‑146 (MQ=255)                                                                                                        cg‑‑ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACga                                                <  5:208034/146‑1 (MQ=255)                                                                                                         g‑‑ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgatgat                                            <  8:78041/149‑1 (MQ=255)                                                                                                            ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGAGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGatgatg                                           <  5:201238/149‑1 (MQ=255)                                                                                                             cTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGt                                          >  4:303045/1‑149 (MQ=255)                                                                                                                     gTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggtg                                  >  7:202922/1‑149 (MQ=255)                                                                                                                      tACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggtgg                                 >  4:306187/1‑149 (MQ=255)                                                                                                                         gccgcGGATGTCAACGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCCACAATTTGGAGTCGAGGATAGAGGTGGTTGTTTGTGAGGTGTATTAGAAAAAACAGCCAAAAGCACGAACCAAGATGTTGCAGGTGGcca                              >  8:108052/1‑149 (MQ=255)                                                                                                                             cgGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAg                          <  3:306187/149‑1 (MQ=255)                                                                                                                              ggATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAgg                         <  4:71954/149‑1 (MQ=255)                                                                                                                                    cACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGGTCGAGGTGGTGGGGTGGGAGGTGAAGTCGATCGCACGGCCAAAGGCACGCACGATGATGTTTCTGGTGGCCCCCAgagatcgg                   >  2:319951/1‑142 (MQ=255)                                                                                                                                    cACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGccacca                           <  1:319951/141‑1 (MQ=255)                                                                                                                                     aCCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCtgtg                  >  4:298975/1‑149 (MQ=255)                                                                                                                                     aCCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCtgtg                  >  3:26333/1‑149 (MQ=255)                                                                                                                                          tGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTgg             <  8:202922/149‑1 (MQ=255)                                                                                                                                             ccacgacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGc           >  5:333474/4‑148 (MQ=255)                                                                                                                                                 gacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGt      >  6:241502/1‑149 (MQ=255)                                                                                                                                                  acgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTAGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTc     >  8:268551/1‑149 (MQ=255)                                                                                                                                                   cgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCg    <  5:241502/149‑1 (MQ=255)                                                                                                                                                    gaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCgg   >  1:22148/1‑149 (MQ=255)                                                                                                                                                     aaGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGc  <  5:81456/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCG‑‑CCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGC  >  NC_002947/160366‑160660

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