Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,994,755	Δ1 bp	8.2%	coding (1482/1599 nt)	PP_5251 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,994,755	0	C	.	8.2%	184.3 / 11.0	61	coding (1482/1599 nt)	PP_5251	transporter
Reads supporting (aligned to +/- strand):  ref base C (23/31);  new base . (3/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

CTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTG  >  NC_002947/5994608‑5994901                                                                                                                                                    |                                                                                                                                                    cTTCCAGCAGGCGGCGCGGGAGCAGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTCTGGCTGGAGTTCTTCTTCGAGTAATCGTGGGTGGAGCATCTGCGACATCacaacc                                                                                                                                                    >  8:73280/1‑149 (MQ=255)    ccAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCg                                                                                                                                                   <  7:247437/147‑1 (MQ=255)     cagcagGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCggg                                                                                                                                                 <  7:239615/148‑1 (MQ=255)        cagGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGAc                                                                                                                                             <  7:173136/149‑1 (MQ=255)          ggcggcGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGAccc                                                                                                                                           <  6:58350/149‑1 (MQ=255)            cggcgCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑gt                                                                                                                                        <  1:18674/149‑1 (MQ=255)             ggcgCGGGAGCTGGAAGCGATGCGAAAGCGCAATGGAGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATGTGAGACATCACCACCGGGTGACCC‑gtg                                                                                                                                       <  7:73280/149‑1 (MQ=255)             ggcgCGGGAGATGGAATCGATGCGCAAGGGCAATGGAGCATTGTCGTGGGGGTTGATTCGGGACAGTGCTGATCCGGCGTTGTGGCTTGAGTTCTTCTTCGAGTAATCGTGGCAGGAGCATCTGCGACATCACCACCGGGTGAGCC‑gtg                                                                                                                                       <  3:42611/149‑1 (MQ=255)                gcgGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTGCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑gtggtg                                                                                                                                    <  1:63515/149‑1 (MQ=255)                          gAAGTGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTg                                                                                                                               >  2:61419/1‑144 (MQ=255)                               gATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAgg                                                                                                                     >  7:1065/1‑149 (MQ=255)                                   cgcaagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGgcgc                                                                                                                  <  8:7073/148‑1 (MQ=255)                                    gcaagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCgg                                                                                                                >  7:143938/1‑149 (MQ=255)                                    gcaagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCgg                                                                                                                <  8:1065/149‑1 (MQ=255)                                        gcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCa                                                                                                            <  4:51321/149‑1 (MQ=255)                                                gcgcATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtgctg                                                                                                    <  5:264825/149‑1 (MQ=255)                                                 cgcATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGa                                                                                                                       >  7:178343/1‑129 (MQ=255)                                                 cgcATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGa                                                                                                                       <  8:178343/129‑1 (MQ=255)                                                    aTTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGa                                                                                                <  8:143938/149‑1 (MQ=255)                                                      tGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAgg                                                                                                          >  1:291797/1‑137 (MQ=255)                                                      tGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAgg                                                                                                          <  2:291797/137‑1 (MQ=255)                                                          gTGGGGGGTGATGCGGGTCAGTGCTGATCCGGCGTTGTGGCTGCTGTTCTTCTTCGAGGACTCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGATGCGCGGGTCAGGATGCTGCAGACCGAt                                                                                           <  5:333570/148‑1 (MQ=255)                                                           tGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCA‑CGGGTCACCCGGTGGTGatg                                                                                                                                 <  3:200605/109‑3 (MQ=255)                                                           tGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCA‑CGGGTCACCCGGTGGTGatg                                                                                                                                 >  4:200605/1‑107 (MQ=255)                                                            gggggTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCgg                                                                                                                                                  <  6:321340/92‑1 (MQ=255)                                                            gggggTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCgg                                                                                                                                                  >  5:321340/1‑92 (MQ=255)                                                               ggTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTGTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATAACCAACGGGTGACCA‑GTGGAGATATGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGtg                                                                                     >  4:29311/1‑148 (MQ=255)                                                               ggTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGtt                                                                                     >  8:52634/1‑149 (MQ=255)                                                                gTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGtt                                                                                     >  5:183022/1‑148 (MQ=255)                                                                 ttGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACAGATGGGttt                                                                                     >  8:112456/1‑147 (MQ=255)                                                                    aTGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGt                                                                                >  1:148451/1‑149 (MQ=255)                                                                     tGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACgggtc                                                                                                                                                >  4:23440/1‑83 (MQ=255)                                                                     tGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACgggtc                                                                                                                                                <  3:23440/85‑3 (MQ=255)                                                                                    aTCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtgctg                                                                                                    <  2:59218/113‑1 (MQ=255)                                                                                    aTCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtgctg                                                                                                    >  1:59218/1‑113 (MQ=255)                                                                                     tCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTAt                                                               >  7:293767/1‑149 (MQ=255)                                                                                      ccGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGc                                                                              >  8:191154/1‑133 (MQ=255)                                                                                      ccGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGc                                                                              <  7:191154/133‑1 (MQ=255)                                                                                       cGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGTCATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTAt                                                               <  3:13062/147‑1 (MQ=255)                                                                                       cGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTAtt                                                              <  7:52634/148‑1 (MQ=255)                                                                                       cGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAAGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTg                                                             <  6:67568/149‑1 (MQ=255)                                                                                          cGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCt                                                          <  2:148451/149‑1 (MQ=255)                                                                                                  cTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCTATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGa                                                  <  1:61419/149‑1 (MQ=255)                                                                                                      agatcttcttcGATGAATCGTGGCTGGAGCATCTGCGACATCACCACAGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGTATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCAt                                              <  8:236790/146‑1 (MQ=255)                                                                                                           ttcttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCAt                                                                   >  5:221445/1‑123 (MQ=255)                                                                                                           ttcttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCAt                                                                   <  6:221445/123‑1 (MQ=255)                                                                                                             cttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAgg                                          <  8:106444/146‑1 (MQ=255)                                                                                                              ttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGcaca                                      >  3:314506/1‑149 (MQ=255)                                                                                                              ttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGAGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGcaca                                      >  1:258916/1‑149 (MQ=255)                                                                                                                 gAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCa                                   >  6:169916/1‑149 (MQ=255)                                                                                                                         gTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGcc                           <  8:293767/149‑1 (MQ=255)                                                                                                                            tcgggcgCATCGGCGAAATCACCAACGGGTGACCC‑GTGGTGAACTGAACAGCGCCGCTCGGGGAAGGATTCTGCCGACCGATGGTTGTGAGGGTCGGATTCCGAATCTATTTGCTGGTGGGGATCATAAGGCACACCATTGATGCCGGc                        <  3:29311/143‑1 (MQ=255)                                                                                                                                gAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGTGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCGGGTGGGGAGCATAAGGCCCACCATTGATGCCGGcccgt                    >  7:75769/1‑147 (MQ=255)                                                                                                                                gAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGAt                              <  4:137293/139‑1 (MQ=255)                                                                                                                                gAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGAt                              >  3:137293/1‑139 (MQ=255)                                                                                                                                 aGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTa                   <  2:258916/149‑1 (MQ=255)                                                                                                                                    aTCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACgg                >  8:185656/1‑149 (MQ=255)                                                                                                                                      cTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGcc              >  2:149138/1‑149 (MQ=255)                                                                                                                                      cTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGGGCGGATGCGGCATCTAGTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGcc              >  4:268318/1‑149 (MQ=255)                                                                                                                                            gttCACCA‑CGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAgg                                          <  4:285338/113‑1 (MQ=255)                                                                                                                                            gttCACCA‑CGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAgg                                          >  3:285338/3‑115 (MQ=255)                                                                                                                                            gttCACCA‑CGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTcg        >  4:197340/3‑149 (MQ=255)                                                                                                                                                accaccGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGTGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCggg    >  3:179139/1‑149 (MQ=255)                                                                                                                                                   accGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGtt  >  4:27741/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTG  >  NC_002947/5994608‑5994901

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑