Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,996,238	C→T	7.3%	R265H (CGT→CAT)	PP_5252 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,996,238	0	C	T	7.3%	242.1 / 4.8	82	R265H (CGT→CAT)	PP_5252	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (42/34);  new base T (3/3);  total (45/37)
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.18e-01

TTCTGGGTGAACAGGTTGTAGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGATCTCGCCACCTG  >  NC_002947/5996094‑5996385                                                                                                                                                   |                                                                                                                                                    ttCTGCGTGAACAGGTTGTAGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCATGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGACGATGGCGCTGGTCAGGCCCA‑GCCGGATGAGTTCACGCAt                                                                                                                                                 >  4:217710/1‑149 (MQ=255)   gtGGGTGAACTGGTTGTAGGCGATGCGTAGGGTCAATTGAGCGTTGTCGGCCAGCTCCTGGATCACGGAGTAGCTTATCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCAt                                                                                                                                                 <  1:279581/147‑1 (MQ=255)   cTGGGTGAACAGGTTGTAGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATg                                                                                                                                                >  2:79041/1‑148 (MQ=255)     gggTGAACAGGTTGTAGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAAc                                                                                                                                             >  1:261388/1‑149 (MQ=255)         gAACAGGTTGTAGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGc                                                                                                                                         >  6:55774/1‑149 (MQ=255)          aaCAGGTTGTAGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCGCA‑GCCGGTTGAGTTCATGCGTGAACtcaac                                                                                                                                        <  1:28400/149‑5 (MQ=255)               gTTGTAGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTATGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCAA‑GCCGGTTGAGTTAACGCATGAACTGGCGAGTg                                                                                                                                    >  8:214661/1‑148 (MQ=255)                  gTAGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGt                                                                                                                                >  1:153402/1‑149 (MQ=255)                   tAGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGt                                                                                                                                <  1:79041/148‑1 (MQ=255)                    aGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGtt                                                                                                                               >  4:191506/1‑148 (MQ=255)                      gCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGAc                                                                                                                            >  8:292494/1‑149 (MQ=255)                      gCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGAc                                                                                                                            >  6:192293/1‑149 (MQ=255)                            cGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATa                                                                                                                      <  2:153402/149‑1 (MQ=255)                                cGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTc                                                                                                                   >  1:237521/1‑148 (MQ=255)                                cGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGGCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTc                                                                                                                   >  4:318679/1‑148 (MQ=255)                                 ggTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTc                                                                                                                   >  5:11400/1‑147 (MQ=255)                                  gTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAg                                                                                                                >  6:47793/1‑149 (MQ=255)                                  gTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAg                                                                                                                >  4:177195/1‑149 (MQ=255)                                   tCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGc                                                                                                               <  2:237521/149‑1 (MQ=255)                                     aaTTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCgg                                                                                                             >  6:179619/1‑149 (MQ=255)                                          aTCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTGTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGc                                                                                                         <  5:192293/148‑1 (MQ=255)                                          aTCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCt                                                                                                        >  1:278403/1‑149 (MQ=255)                                             gTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAg                                                                                                     <  3:196864/149‑1 (MQ=255)                                               tGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAgg                                                                                                    <  7:166945/148‑1 (MQ=255)                                               tGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGAACTGATACTCCAGCGGCAGCTTAgg                                                                                                    >  8:166945/1‑148 (MQ=255)                                                gTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGc                                                                                                   <  2:278403/148‑1 (MQ=255)                                                      caccTCCTGGCGCTCGCAGTAG‑CCGTCGGCGAAGTGCTGGTCCCCGCCGCCGGCGGCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTgg                                                                                            <  4:251296/146‑1 (MQ=255)                                                      cAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTgg                                                                                            >  8:191833/1‑149 (MQ=255)                                                       aGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGc                                                                                           <  5:353667/149‑1 (MQ=255)                                                         cTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGCTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCGA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCa                                                                                         <  3:217710/149‑1 (MQ=255)                                                            cTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCg                                                                                      >  7:219498/1‑149 (MQ=255)                                                               gATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGc                                                                                                   >  5:113822/1‑133 (MQ=255)                                                               gATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGACCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGc                                                                                                   <  6:113822/133‑1 (MQ=255)                                                                aTCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGc                                                                                  <  2:261388/149‑1 (MQ=255)                                                                 tCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGcccc                                                                                                >  3:20426/1‑134 (MQ=255)                                                                 tCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGcccc                                                                                                <  4:20426/134‑1 (MQ=255)                                                                  cACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCg                                                                                 >  2:329514/1‑148 (MQ=255)                                                                     ggAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTaga                                                                             <  6:11400/149‑1 (MQ=255)                                                                     ggAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTaga                                                                             <  5:55774/149‑1 (MQ=255)                                                                     ggAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTaga                                                                             <  7:292494/149‑1 (MQ=255)                                                                      gAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTagag                                                                            <  8:219498/149‑1 (MQ=255)                                                                               gtcgGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCACGCCGGTTGAGTTCACGCATGAACTGGCGa                                                                                                                                       <  1:201398/83‑1 (MQ=255)                                                                               gtcgGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGa                                                                                                                                       >  2:201398/1‑82 (MQ=255)                                                                               gtcgGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATcg                                                                    >  1:71158/1‑149 (MQ=255)                                                                                tcgGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATcgc                                                                   >  6:53435/1‑149 (MQ=255)                                                                                        gTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCAGGCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAggg                                                            >  1:71539/1‑149 (MQ=255)                                                                                        gTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTg                                                                                    >  7:10150/1‑124 (MQ=255)                                                                                        gTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTg                                                                                    <  8:10150/124‑1 (MQ=255)                                                                                        gTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAgggc                                                           >  4:19142/1‑149 (MQ=255)                                                                                        gTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGGTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGCATCATCGCGTTAgggc                                                           >  4:18592/1‑149 (MQ=255)                                                                                         ttCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAgggcg                                                          >  6:14733/1‑149 (MQ=255)                                                                                          tCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTag                                                                              <  2:183862/128‑1 (MQ=255)                                                                                          tCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTAACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTag                                                                              >  1:183862/1‑128 (MQ=255)                                                                                           cTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAgggcgg                                                         <  5:47793/148‑1 (MQ=255)                                                                                              gTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAgggcgggc                                                       >  2:138901/1‑147 (MQ=255)                                                                                              gTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGa                                                     >  1:181092/1‑149 (MQ=255)                                                                                               tAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGc                                                                                                               >  1:327459/1‑90 (MQ=255)                                                                                               tAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGc                                                                                                               <  2:327459/90‑1 (MQ=255)                                                                                                 gccgccgccGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCa                                                  <  7:191833/149‑1 (MQ=255)                                                                                                  ccgccgccgGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAg                                                 <  2:181092/149‑1 (MQ=255)                                                                                                    gccgccgGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCa                                               <  1:138901/149‑1 (MQ=255)                                                                                                           gCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTaga                                                                             >  1:174470/1‑112 (MQ=255)                                                                                                           gCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTaga                                                                             <  2:174470/112‑1 (MQ=255)                                                                                                                gATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCAATCAGCATGCCGGTgggg                                    <  5:179619/148‑1 (MQ=255)                                                                                                                       cTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGaa                            <  2:71158/149‑1 (MQ=255)                                                                                                                        tGGTCAGGCCCA‑GCCGGTTGAGTTAACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTGAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAAc                           >  2:267288/1‑149 (MQ=255)                                                                                                                          gTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAAc                           <  3:177195/147‑1 (MQ=255)                                                                                                                          gTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAAc                           <  3:19142/147‑1 (MQ=255)                                                                                                                               gCCCT‑GCCGGTTGAGTTCACGCATGAAGTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTcgc                    <  3:191506/149‑1 (MQ=255)                                                                                                                                  cccgCCGGTTGAGGTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCgcgc                  <  8:100951/146‑1 (MQ=255)                                                                                                                                      ccGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGa              >  4:315809/1‑149 (MQ=255)                                                                                                                                        tcgccaGTTCATGCGTGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCa                                               >  2:160235/6‑114 (MQ=255)                                                                                                                                        tcgccaGTTCATGCGTGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCa                                               <  1:160235/109‑1 (MQ=255)                                                                                                                                        ggTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGAtc            >  3:337332/1‑149 (MQ=255)                                                                                                                                         gTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGAtct           >  6:65056/1‑149 (MQ=255)                                                                                                                                          tggttttCATGCTTCAAGTGGCTTGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGAtctc          <  5:57809/144‑1 (MQ=255)                                                                                                                                          gccaGTTCATGCGTGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGAtctc          >  4:93946/4‑149 (MQ=255)                                                                                                                                           ccaGTTCATGCGTGAACTGGCGAGTGGAGTTGACCTGATACTACAGCGGCAGCTTAGGCCCCTTGGCAAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGATCTCg         >  3:40183/3‑149 (MQ=255)                                                                                                                                             aGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGATCTCGc        <  5:53435/148‑1 (MQ=255)                                                                                                                                              gTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCACCGGCAGCTTAGGCCCCTTGGACAGCGTGGCTTAGACGATCATCGCGTTAGGGCGGGCGAGCAGCATGCGGGTGGGGTTGCCGTACTGTTCGCGCTGGATCTCGCCa      >  5:243497/1‑149 (MQ=255)                                                                                                                                                tCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGATCTCGCCAcc    >  8:314689/1‑149 (MQ=255)                                                                                                                                                  aCGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGATCTCGCCACCTg  >  2:286832/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TTCTGGGTGAACAGGTTGTAGGCGATGCGTACGGTCAATTGATCGTTGTCGGCCAGCTCCTGGATCACGGAGTAG‑TCGTCGGGGAAGTTCTGGTAGCCGCCGCCGGCGTCGATGGCGCTGGTCAGGCCCA‑GCCGGTTGAGTTCACGCATGAACTGGCGAGTGGAGTTGACCTGATACTCCAGCGGCAGCTTAGGCCCCTTGGCCAGCGTGGCGTAGAGGATCATCGCGTTAGGGCGGGCGATCAGCATGCCGGTGGGGTTGCCGAACTTGTCGCGCTGGATCTCGCCACCTG  >  NC_002947/5996094‑5996385

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