Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,041,127	T→A	7.4%	Q17H (CAA→CAT)	PP_5293 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,041,127	0	T	A	7.4%	236.0 / 7.5	81	Q17H (CAA→CAT)	PP_5293	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/39);  new base A (4/2);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

TGAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTT  >  NC_002947/6040984‑6041275                                                                                                                                                 |                                                                                                                                                      tGAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑ATGcc                                                                                                                                                   >  3:387981/1‑148 (MQ=255) tGAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑ATGcc                                                                                                                                                   >  5:38838/1‑148 (MQ=255) tGAAAGTTCAGCGCCTCTGTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCg                                                                                                                                                  >  1:95062/1‑149 (MQ=255)  gAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCg                                                                                                                                                  >  3:360617/1‑148 (MQ=255)  gAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCg                                                                                                                                                  >  5:353707/1‑148 (MQ=255)   aaaGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGcc                                                                                                                                                   <  7:210944/146‑1 (MQ=255)   aaaGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGa                                                                                                                                                >  2:52792/1‑149 (MQ=255)     aGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGAtt                                                                                                                                              >  2:27054/1‑149 (MQ=255)         cAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCg                                                                                                                                           <  3:115219/148‑1 (MQ=255)         cAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATATCgg                                                                                                                                          >  6:208710/1‑149 (MQ=255)            cgcCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CAt                                                                                                                                      >  3:278351/1‑149 (MQ=255)               ctctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTgg                                                                                                                                   >  2:381631/1‑149 (MQ=255)               ctctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑ATGCCGAAATCCGG‑CAATgg                                                                                                                                   <  6:38838/149‑1 (MQ=255)                   ttGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCgg                                                                                                                               >  6:183925/1‑149 (MQ=255)                    tGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑ATGCCGAAATCCGG‑CAATGGGCga                                                                                                                               <  4:387981/148‑2 (MQ=255)                     gCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCggg                                                                                                                              <  4:28972/148‑1 (MQ=255)                      cccTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCgggg                                                                                                                             <  4:355358/148‑1 (MQ=255)                       ccTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGt                                                                                                                            >  4:197645/1‑148 (MQ=255)                       ccTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGt                                                                                                                            >  2:59240/1‑148 (MQ=255)                           ggCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATc                                                                                                                        <  4:5667/148‑1 (MQ=255)                            gCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGa                                                                                                                      <  5:395214/149‑1 (MQ=255)                             cGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGAc                                                                                                                     <  3:321317/149‑1 (MQ=255)                              gTCGAAAAAAGGGGCTTGCTCGCGCTTGAGGTGCCATACCACCAGCGGGCCCACCAGGGGCCCAAGCGGCACTACCACGCCCAACAGCGCGGAGAGGGGGCACAACATCCCCCATTTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGAc                                                                                                                     <  1:52792/149‑1 (MQ=255)                               tCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAg                                                                                                                   <  3:341261/149‑1 (MQ=255)                                   tAATAGGGTCTTGCTCGCGGTTGAGGTGCCATAGCACCTGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACTTCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTc                                                                                                               <  1:87349/149‑1 (MQ=255)                                    aaaaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCg                                                                                                              <  1:27054/149‑1 (MQ=255)                                    aaaaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGACGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTc                                                                                                               >  7:89245/1‑148 (MQ=255)                                     aaaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCg                                                                                                              >  7:251752/1‑148 (MQ=255)                                      aaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGtt                                                                                                                <  4:129437/145‑1 (MQ=255)                                      aaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGtt                                                                                                                >  3:129437/1‑145 (MQ=255)                                      aaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGaa                                                                                                            >  5:311489/1‑149 (MQ=255)                                      aaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGaa                                                                                                            >  3:355438/1‑149 (MQ=255)                                      aaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGaa                                                                                                            <  8:251752/149‑1 (MQ=255)                                      aaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGaa                                                                                                            >  2:58874/1‑149 (MQ=255)                                        gggTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATc                                                                                                          <  7:128152/149‑1 (MQ=255)                                          gTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGc                                                                                                        <  5:226297/149‑1 (MQ=255)                                               gCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGGCCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTc                                                                                                      <  1:58874/147‑1 (MQ=255)                                               gCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCa                                                                                                     <  8:26842/147‑1 (MQ=255)                                                 tCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGcc                                                                                                 >  5:250669/1‑149 (MQ=255)                                                      ccgctggTGCTATGGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCa                                                                                            >  8:335065/6‑149 (MQ=255)                                                           ggTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTc                                                                                                               <  7:79338/125‑1 (MQ=255)                                                           ggTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTc                                                                                                               >  8:79338/1‑125 (MQ=255)                                                           ggTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctcc                                                                                       <  6:28393/149‑1 (MQ=255)                                                               ccATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTgg                                                                                   <  6:353707/149‑1 (MQ=255)                                                                  tAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCg                                                                                <  7:28964/149‑1 (MQ=255)                                                                    gCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGt                                                                              >  3:314420/1‑149 (MQ=255)                                                                     gaCCTGCGGGCCGAGCAGGTGCCCAAGCGGCACTACGAGGCCCAACAGCGCTGAGAGGTGGCATAACATCGCCCA‑TTGGCGGATTTCGG‑CATTGGGCGGGGTGAGCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGtt                                                                             <  3:112842/148‑1 (MQ=255)                                                                     cACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGtt                                                                                                                <  2:85528/114‑1 (MQ=255)                                                                     cACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGtt                                                                                                                >  1:85528/1‑114 (MQ=255)                                                                      aCCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCggg                                                                                                                              >  8:318423/1‑99 (MQ=255)                                                                      aCCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCggg                                                                                                                              <  7:318423/99‑1 (MQ=255)                                                                         aGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctcc                                                                                       <  6:245459/135‑1 (MQ=255)                                                                         aGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATGTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCAcctcc                                                                                       >  5:245459/1‑135 (MQ=255)                                                                          gCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtcagt                                                                        >  8:306926/1‑149 (MQ=255)                                                                           cGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtcagtc                                                                       <  3:197645/149‑1 (MQ=255)                                                                                   gcagGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtg                                                               <  5:183925/149‑1 (MQ=255)                                                                                       gTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtgctgc                                                           >  1:256260/1‑149 (MQ=255)                                                                                           ccAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGctgcctgc                                                       <  7:306926/149‑1 (MQ=255)                                                                                            cAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCctgct                                                      >  3:117688/1‑149 (MQ=255)                                                                                                    aCTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTc                                               >  2:180113/1‑148 (MQ=255)                                                                                                       aCCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCCTGCCCCGCACCGCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGaa                                            >  7:85119/1‑148 (MQ=255)                                                                                                        ccAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGaaa                                           <  5:43105/148‑1 (MQ=255)                                                                                                           ggCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtct                                       <  1:381631/149‑1 (MQ=255)                                                                                                            gCCCAACAGCGCGGCGAGGGGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtctc                                      <  3:305587/149‑1 (MQ=255)                                                                                                            gCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtctc                                      <  6:311489/149‑1 (MQ=255)                                                                                                                   aGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCTACAGGTTCGATTCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCTGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATg                               <  1:59240/149‑1 (MQ=255)                                                                                                                        ggAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCtt                          <  2:331826/149‑1 (MQ=255)                                                                                                                              gTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGaa                                                                                                            <  1:213690/61‑1 (MQ=255)                                                                                                                              gTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGaa                                                                                                            >  2:213690/1‑61 (MQ=255)                                                                                                                              gTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTgcgc                    >  5:323489/1‑149 (MQ=255)                                                                                                                                 gCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGaaa                                           >  8:349522/1‑123 (MQ=255)                                                                                                                                   aGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGGCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAg               >  1:323200/1‑149 (MQ=255)                                                                                                                                   aGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGATAATCTCGTTCATGCCCTTCTGCGCCAAAg               >  7:105175/1‑149 (MQ=255)                                                                                                                                    gAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGc              <  6:323489/149‑1 (MQ=255)                                                                                                                                       cATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAg           >  5:115535/1‑149 (MQ=255)                                                                                                                                       cATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGGTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAg           <  5:264169/149‑1 (MQ=255)                                                                                                                                        cccGCCCA‑ATGCCGAAATCCGG‑CAATGGGCGGGGTGATCGACAGGTTCGAATCGCTCATCCCCCGCACCTCCTTCGCCGGTTTCACTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCCTTCATGCCCTTCTGCGCCAAAGCGAg           >  3:312342/3‑148 (MQ=255)                                                                                                                                        aTCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGc          <  2:323200/149‑1 (MQ=255)                                                                                                                                           gCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATg       >  2:26461/1‑149 (MQ=255)                                                                                                                                           gCCCA‑ATGCCGAAATCCGG‑CAATGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATg       >  1:381698/1‑149 (MQ=255)                                                                                                                                            cccA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATgg      <  5:208710/149‑1 (MQ=255)                                                                                                                                                 ttGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGtt  <  4:117688/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      TGAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCA‑TTGCCGGATTTCGG‑CATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCGTT  >  NC_002947/6040984‑6041275

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