Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	262,752	T→A	8.2%	L90Q (CTG→CAG)	PP_0211 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	262,752	0	T	A	8.2%	180.1 / 5.2	61	L90Q (CTG→CAG)	PP_0211	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base A (3/2);  total (32/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 = 9.88e-01

CCGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGA  >  NC_002947/262606‑262900                                                                                                                                                   |                                                                                                                                                      cagTGCGGCGCGCCTGTCCGGAGTGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATTGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTg                                                                                                                                                      <  8:171337/146‑1 (MQ=255)   gTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGt                                                                                                                                                    >  3:111897/1‑148 (MQ=255)            gcCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGa                                                                                                                                         <  2:203271/149‑1 (MQ=255)              cTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAg                                                                                                                                       >  6:80587/1‑149 (MQ=255)               tGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCGAGGGTATGGCGTACAACTGGTGTTCCGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGc                                                                                                                                      >  7:212117/1‑149 (MQ=255)                 tCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGc                                                                                                                                      <  5:231968/147‑1 (MQ=255)                      aGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATccc                                                                                                                               <  3:169823/149‑1 (MQ=255)                       gggATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATc                                                                                                                                 <  6:318916/146‑1 (MQ=255)                       gggATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATc                                                                                                                                 >  5:318916/1‑146 (MQ=255)                               ccTGGTGCGCGACGATGTGCGCGTGGAGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACACTAGGGCTTACTAGCGTGACCTGGACA‑GCAAGTGAAGCTGATCCCTTCAGTAGc                                                                                                                      >  2:269869/1‑149 (MQ=255)                                cTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGcc                                                                                                                     <  4:111897/149‑1 (MQ=255)                                     gcgcGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGC‑TGACCAGGTCacgcaa                                                                                                                                             >  4:245551/1‑115 (MQ=255)                                     gcgcGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGC‑TGACCAGGTCacgcaa                                                                                                                                             <  3:245551/120‑6 (MQ=255)                                         gacgaTGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGt                                                                                                            >  5:269637/1‑149 (MQ=255)                                          acgaTGTGCGCGTGGTGCAAATCGAGTCGCGAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTg                                                                                                           <  8:212117/149‑1 (MQ=255)                                           cgaTGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTg                                                                                                           >  3:89653/1‑148 (MQ=255)                                                tgCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGaa                                                                                                      <  6:33438/148‑1 (MQ=255)                                                    cgTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAg                                                                                                  >  7:117208/1‑148 (MQ=255)                                                     gtggtgCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGga                                                                                                <  6:125550/149‑1 (MQ=255)                                                     gtggtgCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGga                                                                                                >  4:64632/1‑149 (MQ=255)                                                          gCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGAcc                                                                                           >  2:123866/1‑149 (MQ=255)                                                           cAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCg                                                                                          <  6:269637/149‑1 (MQ=255)                                                           cAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCg                                                                                          >  8:185293/1‑149 (MQ=255)                                                            aaaTCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGc                                                                                         <  2:185105/149‑1 (MQ=255)                                                              aTCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGAt                                                                                                                                  >  6:97831/1‑106 (MQ=255)                                                              aTCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGAt                                                                                                                                  <  5:97831/106‑1 (MQ=255)                                                               tCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTa                                                                                      <  7:185293/149‑1 (MQ=255)                                                                  aGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGc                                                                                    >  6:178061/1‑148 (MQ=255)                                                                      gCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAAc                                                                                <  1:123866/148‑1 (MQ=255)                                                                                 gCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGa                                                                     <  2:294525/148‑1 (MQ=255)                                                                                         cTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCt                                                            >  3:110373/1‑149 (MQ=255)                                                                                                    cGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCg                                                 <  8:117208/149‑1 (MQ=255)                                                                                                       tGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGaca                                              >  6:347657/1‑149 (MQ=255)                                                                                                           cATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACAcccc                                          >  5:117885/1‑149 (MQ=255)                                                                                                           cATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTAGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGAcccccc                                          >  1:292505/1‑149 (MQ=255)                                                                                                             tGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGa                                        <  3:64632/149‑1 (MQ=255)                                                                                                              gAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCa                                                                                                                           >  4:120182/1‑65 (MQ=255)                                                                                                              gAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCa                                                                                                                           <  3:120182/65‑1 (MQ=255)                                                                                                               aGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAgt                                      <  5:178061/149‑1 (MQ=255)                                                                                                                gcgcGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAgtg                                     >  5:78367/1‑149 (MQ=255)                                                                                                                   cgGGATCTACCCTTGGGCTTTCTTGCGTGAGCTGGTCA‑GGAAGTGAAGCTGATCCCTTCAGTTGCCGGACGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACGTTCTCTTTGGCGACACCCCGAGTGTgg                                  <  1:269869/149‑1 (MQ=255)                                                                                                                   cgGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTgg                                  <  3:164738/149‑1 (MQ=255)                                                                                                                   cgGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTgg                                  <  4:110373/149‑1 (MQ=255)                                                                                                                      gATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTa                               >  8:121012/1‑149 (MQ=255)                                                                                                                           aCCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCAGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAggg                          >  2:322732/1‑149 (MQ=255)                                                                                                                              cTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGgcg                        >  4:144894/1‑148 (MQ=255)                                                                                                                              cTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCAGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGgcg                        >  2:226701/1‑148 (MQ=255)                                                                                                                                 gggCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCaa                     <  6:78367/148‑1 (MQ=255)                                                                                                                                  ggCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGctc                                                        >  6:158682/1‑112 (MQ=255)                                                                                                                                  ggCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAAc                    <  3:144894/148‑1 (MQ=255)                                                                                                                                   cttcaCTTGC‑TGACCAGGTCACGCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTg                  >  6:326456/5‑149 (MQ=255)                                                                                                                                    ttcaCTTGC‑TGACCAGGTCACGCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTgg                 >  5:286978/4‑149 (MQ=255)                                                                                                                                    ttcaCTTGC‑TGACCAGGTCACGCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTgg                 <  5:326456/146‑1 (MQ=255)                                                                                                                                    cTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGACCTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTgg                 <  2:35236/149‑1 (MQ=255)                                                                                                                                    cTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTgg                 <  5:347657/149‑1 (MQ=255)                                                                                                                                         ttGCGTGACCTGGTCA‑GTAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAAc            >  5:265893/1‑149 (MQ=255)                                                                                                                                            cGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTg         >  7:287136/1‑149 (MQ=255)                                                                                                                                             gTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGc        >  4:57153/1‑149 (MQ=255)                                                                                                                                             gTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGc        >  4:151565/1‑149 (MQ=255)                                                                                                                                                aCCTGGTCA‑GTAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCAc     <  6:265893/149‑1 (MQ=255)                                                                                                                                                  cTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCAcc    >  2:217954/1‑148 (MQ=255)                                                                                                                                                   tGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGa  >  6:302277/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CCGTGCGGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGA  >  NC_002947/262606‑262900

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