Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,768,071	G→A	6.6%	G29D (GGC→GAC)	PP_5600 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,768,071	0	G	A	6.6%	179.1 / 2.1	60	G29D (GGC→GAC)	PP_5600	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base A (2/2);  total (32/28)
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.79e-01

AAGATGTAGAATCCCTCGCCAACCTCCCAACCTGATCATTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCT  >  NC_002947/3767923‑3768215                                                                                                                                                     |                                                                                                                                                 aagaTGTAGAATCCCTCGCCAACCTCCCAACCTGATCATTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCgccatggg                                                                                                                                                  >  2:385062/1‑149 (MQ=255)   gaTGTAGAATCCCTCGCCAACCTCCCAACCTGATCATTCAGTTTTCACGCCGAGTCTCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCgccatgggc                                                                                                                                                 >  5:323268/1‑148 (MQ=255)              ccTCGCCAACCTCCCAACCTGATCATTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAAt                                                                                                                                      <  6:215500/148‑1 (MQ=255)              ccTCGCCAACCTCCCAACCTGATCATTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAAt                                                                                                                                      <  4:166621/148‑1 (MQ=255)              acacGCCAACCTCCCAAACTGATGATTCACTTTTCACGCCGAGACCCTCCATGAAGCTTTCACACAACACACGATCACATCTCCCCTCAGTGCTCTATTTAAGCATCCTGTTCGGCTCGTTTCAGTGAGCCATGGGCCATGGGCAAAtc                                                                                                                                     >  2:59148/4‑148 (MQ=255)                    cAACCTCCCACCCTGATCATTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGTGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATTTGCCCATGGcccat                                                                                                                               <  2:104505/149‑4 (MQ=255)                    cAACCTCCCAACCTGATCTTTCAGTTTTCACGGCGTGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTTTTTGAGCATCCTGTTCGGCGCGGTGCTGGGCGCCAGGGGCCATGGGCAAATGGCCggg                                                                                                                               <  1:59148/149‑1 (MQ=255)                    cAACCTCCCAACCTGATCATTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCggg                                                                                                                               <  8:220904/149‑1 (MQ=255)                          cccAACCTGATCATTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTgcgc                                                                                                                         >  7:42027/1‑149 (MQ=255)                          cccAACCTGATCATTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTgcgc                                                                                                                         >  3:188135/1‑149 (MQ=255)                            cAACCTGATCATTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCct                                                                                                                       <  6:193165/149‑1 (MQ=255)                                      aTTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTAt                                                                                                              <  1:173968/148‑1 (MQ=255)                                      aTTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATg                                                                                                             <  7:358047/149‑1 (MQ=255)                                       ttCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGACCATGGGCAAATGGCCGGGTTGCGCctggc                                                                                                                    >  5:207783/1‑141 (MQ=255)                                       ttCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGACCATGGGCAAATGGCCGGGTTGCGCctggc                                                                                                                    >  5:202634/1‑141 (MQ=255)                                       ttCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGACCATGGGCAAATGGCCGGGTTGCGCctggc                                                                                                                    <  6:202634/141‑1 (MQ=255)                                       ttCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGACCATGGGCAAATGGCCGGGTTGCGCctggc                                                                                                                    <  6:207783/141‑1 (MQ=255)                                        tCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGga                                                                                                           >  6:46050/1‑149 (MQ=255)                                                  cGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATgccgcc                                                                                                 >  3:68331/1‑149 (MQ=255)                                                   gCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATgcc                                                                                                    <  8:299689/145‑1 (MQ=255)                                                   gCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATgcc                                                                                                    >  7:299689/1‑145 (MQ=255)                                                           ccGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACg                                                                                        >  6:92517/1‑149 (MQ=255)                                                           ccGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCAGTCCGCCATGGGCCATGGGCAAATGGCCGGGTAGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACg                                                                                        >  2:42767/1‑149 (MQ=255)                                                                aTGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCctggc                                                                                                                    <  1:189442/116‑1 (MQ=255)                                                                aTGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCctggc                                                                                                                    >  2:189442/1‑116 (MQ=255)                                                                         gCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCgg                                                                          >  2:240534/1‑149 (MQ=255)                                                                           ccaccaCACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCa                                                                        >  3:108203/1‑149 (MQ=255)                                                                            caccacACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCAt                                                                       >  8:345888/1‑149 (MQ=255)                                                                             accacACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATg                                                                      >  5:257269/1‑149 (MQ=255)                                                                             accacACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGGTGCGCCTGGCTGGTATGGATGATGCCGCCTTGTGCACGGCCGACAGCAGCGGCATg                                                                      <  4:51407/149‑1 (MQ=255)                                                                                 cacacGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTcc                                                                  <  4:68331/149‑1 (MQ=255)                                                                                      gATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGga                                                                                                           >  2:195723/1‑103 (MQ=255)                                                                                      gATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGga                                                                                                           <  1:195723/103‑1 (MQ=255)                                                                                      gATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGt                                                             <  2:34969/149‑1 (MQ=255)                                                                                       aTCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGt                                                             >  1:73454/1‑148 (MQ=255)                                                                                        tCACATGTCGCCTGGGTGCTGTCTTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGa                                                           <  8:305354/149‑1 (MQ=255)                                                                                        tCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGt                                                             >  3:188311/1‑147 (MQ=255)                                                                                              gTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGa                                                     >  2:227088/1‑149 (MQ=255)                                                                                                 gCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGc                                                                                      <  1:283740/113‑1 (MQ=255)                                                                                                 gCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGc                                                                                      >  2:283740/1‑113 (MQ=255)                                                                                                  ccTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACtt                                                  >  8:312159/1‑148 (MQ=255)                                                                                                  ccTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTc                                                 <  1:240534/149‑1 (MQ=255)                                                                                                   cTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCa                                                <  8:53560/149‑1 (MQ=255)                                                                                                      ggTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAg                                               <  8:42027/147‑1 (MQ=255)                                                                                                      ggTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAg                                               <  7:345888/147‑1 (MQ=255)                                                                                                       gTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGa                                            >  7:57628/1‑149 (MQ=255)                                                                                                        tGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGAt                                           >  5:186343/1‑149 (MQ=255)                                                                                                         gCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATc                                          >  8:227620/1‑149 (MQ=255)                                                                                                           tGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATcc                                         >  8:21299/1‑148 (MQ=255)                                                                                                            gTATTACAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGt                                                                                             >  7:316637/1‑95 (MQ=255)                                                                                                            gTATTACAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGt                                                                                             <  8:316637/95‑1 (MQ=255)                                                                                                              aTTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGccc                                     >  1:316992/1‑149 (MQ=255)                                                                                                                  cAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTc                                                 >  4:317557/1‑133 (MQ=255)                                                                                                                  cAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTc                                                 <  3:317557/133‑1 (MQ=255)                                                                                                                  cAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGt                                 >  7:325043/1‑149 (MQ=255)                                                                                                                            ttCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCcggcg                        <  1:385062/148‑1 (MQ=255)                                                                                                                                    gTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGc                <  6:257269/148‑1 (MQ=255)                                                                                                                                     ttGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGt              >  1:256503/1‑149 (MQ=255)                                                                                                                                        cTGTGCGCCATGGGGCATGGGCAAATCGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCACCTGGTGCACGGCCGACAGCAGCGGCATGGACCAAGGAGAGCAAGACTTCAGCGATCCGACCGTGTCTGCCGGCGGTGACTGCCTGAt           >  1:110558/1‑149 (MQ=255)                                                                                                                                          gtgCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCg         >  3:186691/1‑149 (MQ=255)                                                                                                                                           tgCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATc          <  7:21299/147‑1 (MQ=255)                                                                                                                                             cgcCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCa      <  5:92517/149‑1 (MQ=255)                                                                                                                                              gccatgggccatggCCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATCGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAg     >  8:265440/1‑149 (MQ=255)                                                                                                                                                 atgggccatgggcAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCGTGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCt  >  3:105696/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AAGATGTAGAATCCCTCGCCAACCTCCCAACCTGATCATTCAGTTTTCACGCCGAGTCCCGCCATGAAGCTGGCCCACCACACACGATCACATGTCGCCTGGGTGCTGTATTTCAGCATCCTGTTCGGCTCGTTGCTGTGCGCCATGGGCCATGGGCAAATGGCCGGGTTGCGCCTGGCTGGTATGGATGATGCCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCT  >  NC_002947/3767923‑3768215

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