Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,048,828	G→T	6.9%	intergenic (+67/+72)	PP_1821 → / ← PP_1822	glutathione S‑transferase family protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,048,828	0	G	T	6.9%	173.0 / 3.5	58	intergenic (+67/+72)	PP_1821/PP_1822	glutathione S‑transferase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base T (2/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

GACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCA  >  NC_002947/2048681‑2048959                                                                                                                                                     |                                                                                                                                    gACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGgcgc                                                                                                                                    >  5:176944/1‑149 (MQ=255)      ggCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGccc                                                                                                                               <  3:12868/149‑1 (MQ=255)        cAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGt                                                                                                                             <  5:321728/149‑1 (MQ=255)         aGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGCGGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGt                                                                                                                             >  1:288224/1‑149 (MQ=255)            cTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGcc                                                                                                                          <  2:177216/148‑1 (MQ=255)             tagcgggcggTAAGGGGGTGGCTGGAGCGGGTCAACAGCAACCCGC‑GGCATGTGCGGATGGTTGGCTGAAATCCTACCTGTTTTCCTGGCGCTTTCTCAATCCTTGTGGGAGCTGGAGTGCCAGCGAACACGGGCGCAGCCCGTGc                                                                                                                           <  4:232967/141‑1 (MQ=255)             tACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAACAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGc                                                                                                                           >  3:232967/1‑146 (MQ=255)               cccGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGCCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAac                                                                                                                       >  6:10491/1‑148 (MQ=255)                ccGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAacca                                                                                                                     >  6:230296/1‑149 (MQ=255)                ccGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCCGGCGTGCCCGCGAACACCGGCGCAGCCCGTGCCAacca                                                                                                                     >  2:183698/1‑149 (MQ=255)                           ggCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCt                                                                                                          <  4:372135/149‑1 (MQ=255)                            gCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTg                                                                                                         <  7:137830/149‑1 (MQ=255)                             ccTGGCTGGAGCGGGTCAGCATGCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGTCGCAGCCCGTGCCAACCACCGCGGTGCCTGc                                                                                                        <  6:331673/149‑1 (MQ=255)                             ccTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGc                                                                                                        >  8:350615/1‑149 (MQ=255)                               tggctggAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCTGCGCCCGTg                                                                                                                            >  8:336566/1‑127 (MQ=255)                               tggctggAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCTGCGCCCGTg                                                                                                                            >  8:274684/1‑127 (MQ=255)                               tggctggAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCTGCGCCCGTg                                                                                                                            <  7:336566/127‑1 (MQ=255)                               tggctggAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCTGCGCCCGTg                                                                                                                            <  7:274684/127‑1 (MQ=255)                                ggctggAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTc                                                                                                     <  2:297213/149‑1 (MQ=255)                                   tggAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTcgcg                                                                                                  >  8:274998/1‑149 (MQ=255)                                     gAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCggg                                                                                                >  8:169931/1‑149 (MQ=255)                                       gCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCt                                                                                              >  5:336143/1‑149 (MQ=255)                                          gatgtaggaggccactgcaatcatcgctttcctggatgagcggtttgccgccatctgactgacgctgcagatcCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTc                                                                                             >  4:267652‑M1/74‑148 (MQ=255)                                                  gcCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTccc                                                                                   <  4:305533/149‑1 (MQ=255)                                                  agtccaaggcaatcatcgcttacctggatgagcggtttgccgccatctgactgacgctgcagatcCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTc                                                                                     >  1:315255‑M1/66‑148 (MQ=255)                                                     aCCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCCGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCcac                                                                                 >  5:107008/1‑148 (MQ=255)                                                     aCCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCcac                                                                                 >  4:178250/1‑148 (MQ=255)                                                       ccGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAg                                                                              <  1:220852/149‑1 (MQ=255)                                                       ccGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAg                                                                              >  3:59421/1‑149 (MQ=255)                                                       ccGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAg                                                                              <  3:178250/149‑1 (MQ=255)                                                        cgc‑gGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAgg                                                                             >  5:88668/1‑149 (MQ=255)                                                                  tgCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCa                                                                     <  5:230296/148‑1 (MQ=255)                                                                   gCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGc                                                                   >  3:83631/1‑149 (MQ=255)                                                                                 gAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCgtcagt                                                             >  1:165445/1‑141 (MQ=255)                                                                                 gAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCgtcagt                                                             <  2:165445/141‑1 (MQ=255)                                                                                  aaaaCCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATggc                                                     <  7:350615/148‑1 (MQ=255)                                                                                        gACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAAc                                              <  7:169931/149‑1 (MQ=255)                                                                                         aCCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAAcc                                             <  4:295855/149‑1 (MQ=255)                                                                                          cgccatctgactgacgctgcagatcCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAgg                                                                             >  2:40255‑M1/26‑116 (MQ=255)                                                                                            tGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCt                                          >  2:212789/1‑149 (MQ=255)                                                                                              ttttCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCa                                        >  3:297019/1‑149 (MQ=255)                                                                                                   cAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAg                                   <  2:288224/149‑1 (MQ=255)                                                                                                    agccgcTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAgg                                  <  3:199077/146‑1 (MQ=255)                                                                                                        gcTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAg                              <  6:176944/149‑1 (MQ=255)                                                                                                                aaTCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGAtt                      >  2:321867/1‑149 (MQ=255)                                                                                                                       GTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTg                >  4:83631‑M1/1‑148 (MQ=255)                                                                                                                        TGGGTGGGGGCGGGCCCGCTAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGa              >  1:183698‑M1/1‑149 (MQ=255)                                                                                                                           GAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACt            <  5:209418‑M1/148‑1 (MQ=255)                                                                                                                           GAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTc           >  5:10491‑M1/1‑149 (MQ=255)                                                                                                                            AGCGGGCGGGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCa          >  6:336143‑M1/1‑149 (MQ=255)                                                                                                                             GCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCAt         >  1:212789‑M1/1‑149 (MQ=255)                                                                                                                                ggCGTGCCCGCGAACACGGCCGCAGCCCGTGACAACCACCGCGGTGCCTCCTTCCCCGGCTCGCCCCCTCCCACAAGGGATCTGCAGCCTCAGACACATGGCGCCAAACCCCTCATCCAGCTAAGCGATGATACCCTTGGACTCAAAc       >  8:304719/1‑148 (MQ=255)                                                                                                                                GGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACa      <  4:156545‑M1/149‑1 (MQ=255)                                                                                                                                   GTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTGAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATc    <  1:189477‑M1/148‑1 (MQ=255)                                                                                                                                    tGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCa  <  1:321867/149‑1 (MQ=255)                                                                                                                                       ccGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCa                                    >  6:323437/1‑112 (MQ=255)                                                                                                                                       ccGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCa                                    <  5:323437/112‑1 (MQ=255)                                                                                                                                                 gggCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCAt         <  7:148142/129‑1 (MQ=255)                                                                                                                                                 gggCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCAt         >  8:148142/1‑129 (MQ=255)                                                                                                                                                     |                                                                                                                                    GACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCA  >  NC_002947/2048681‑2048959

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