Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,048,829	C→G	6.9%	intergenic (+68/+71)	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,829	0	C	G	6.9%	185.4 / 4.4	58	intergenic (+68/+71)	PP_1821/PP_1822	glutathione S‑transferase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/25);  new base G (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 = 9.38e-01

GACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTC  >  NC_002947/2048681‑2048977                                                                                                                                                      |                                                                                                                                                     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                                                                                                                                              <  7:274684/127‑1 (MQ=255)                               tggctggAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCTGCGCCCGTg                                                                                                                                              <  7:336566/127‑1 (MQ=255)                               tggctggAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCTGCGCCCGTg                                                                                                                                              >  8:274684/1‑127 (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                                                                                                >  3:59421/1‑149 (MQ=255)                                                       ccGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAg                                                                                                <  3:178250/149‑1 (MQ=255)                                                       ccGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAg                                                                                                <  1:220852/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)                                                                                                                                                      cAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGcttc  <  6:107008/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GACCTGGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTC  >  NC_002947/2048681‑2048977

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