Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,048,830	A→C	8.1%	intergenic (+69/+70)	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,830	0	A	C	8.1%	161.7 / 4.2	59	intergenic (+69/+70)	PP_1821/PP_1822	glutathione S‑transferase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base C (3/2);  total (31/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.95e-01

GGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCT  >  NC_002947/2048686‑2048978                                                                                                                                                  |                                                                                                                                                     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                                                                                <  2:165445/141‑1 (MQ=255)                                                                            gAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCgtcagt                                                                                >  1:165445/1‑141 (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)                                                                                                                                                  aGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGcttct  <  6:88668/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GGCAGACTACCCGGCGGTAAGGGCCTGGCTGGAGCGGGTCAGCAGCCACCCGC‑GGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACACGGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCT  >  NC_002947/2048686‑2048978

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