Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,033,783	G→C	8.6%	V36L (GTG→CTG)	dfp →	bifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,033,783	0	G	C	8.6%	209.7 / 7.0	67	V36L (GTG→CTG)	dfp	bifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base C (4/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

AAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGCC  >  NC_002947/6033639‑6033917                                                                                                                                                  |                                                                                                                                       aaaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTgcg                                                                                                                                    >  4:270700/1‑149 (MQ=255) aaaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTgc                                                                                                                                     <  4:190931/148‑1 (MQ=255) aaaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTgc                                                                                                                                     <  3:82929/148‑1 (MQ=255)  aaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTgcgc                                                                                                                                   >  4:12860/1‑149 (MQ=255)  aaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGACGCGCACCTgcgc                                                                                                                                   >  2:309896/1‑149 (MQ=255)     gtcgtcGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCCGTCTCCTGGAGCACGGCGCGCAGGTGCGcgt                                                                                                                                 <  1:47920/149‑1 (MQ=255)       cgtcgCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGACGCGCACCTGCGcgccgt                                                                                                                              <  1:309896/149‑1 (MQ=255)             ccGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCgtcgtc                                                                                                                             >  4:58424/1‑144 (MQ=255)             ccGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCgtcgtc                                                                                                                             <  3:58424/144‑1 (MQ=255)             ccGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGAc                                                                                                                        >  5:57272/1‑149 (MQ=255)             ccGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGACGCGCACCTGCGcgccgtgc                                                                                                                            >  6:138448/1‑143 (MQ=255)             ccGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGACGCGCACCTGCGcgccgtgc                                                                                                                            <  5:138448/145‑3 (MQ=255)                 tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGcgtc                                                                                                                                <  4:141464/137‑1 (MQ=255)                 tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGcgtc                                                                                                                                >  3:141464/1‑137 (MQ=255)                     aCTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtg                                                                                                                >  1:299987/1‑149 (MQ=255)                       tGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtggt                                                                                                              >  8:340005/1‑149 (MQ=255)                         gAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACGAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGgtggtg                                                                                                             >  8:117467/1‑148 (MQ=255)                         gAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtggt                                                                                                              <  3:12860/147‑1 (MQ=255)                          aaaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCa                                                                                                           <  5:341138/149‑1 (MQ=255)                           aaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCGTGACCCGGGGTGGTGCa                                                                                                           <  4:343676/148‑1 (MQ=255)                             aaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGaa                                                                                                        >  2:145218/1‑149 (MQ=255)                             aaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGaa                                                                                                        <  3:305524/149‑1 (MQ=255)                                  aaGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATtcat                                                                                                   <  1:166490/149‑1 (MQ=255)                                  aaGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATtcat                                                                                                   <  3:270700/149‑1 (MQ=255)                                        aTGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCAcccc                                                                                             >  5:240827/1‑149 (MQ=255)                                        aTGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CAGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACACGGGGTGGTGCAGAATTCATCAcccc                                                                                             >  5:246061/1‑149 (MQ=255)                                            aGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCt                                                                                          >  2:113085/1‑148 (MQ=255)                                             gCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGa                                                                                        >  2:247614/1‑149 (MQ=255)                                               ggCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGAcc                                                                                      <  4:125224/149‑1 (MQ=255)                                                    tATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGc                                                                                  <  1:145218/148‑1 (MQ=255)                                                      tCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAgg                                                                               <  6:3784/149‑1 (MQ=255)                                                          aaGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCt                                                                           <  2:129743/149‑1 (MQ=255)                                                            gcgcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGt                                                                         <  6:57272/149‑1 (MQ=255)                                                                aTCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                                   >  3:198417/1‑103 (MQ=255)                                                                aTCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                                   <  4:198417/103‑1 (MQ=255)                                                                         ggCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGt                                                            <  6:246061/149‑1 (MQ=255)                                                                         ggCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGt                                                            <  6:240827/149‑1 (MQ=255)                                                                                gtggCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATg                                                     <  7:340005/149‑1 (MQ=255)                                                                                   gcggcATTGCCGCCTACGAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGAc                                                  <  7:117467/149‑1 (MQ=255)                                                                                   gcggcATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCCTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGAc                                                  <  1:113085/149‑1 (MQ=255)                                                                                     ggcAATGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGAc                                                  <  1:176359/147‑1 (MQ=255)                                                                                     ggcAATGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGAc                                                  >  2:176359/1‑147 (MQ=255)                                                                                        aTTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCt                                             <  1:247614/149‑1 (MQ=255)                                                                                        aTTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCt                                             <  3:13383/149‑1 (MQ=255)                                                                                              gccTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGAccc                                       >  3:17539/1‑149 (MQ=255)                                                                                                  aCAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGgcc                                   >  8:128966/1‑149 (MQ=255)                                                                                                   cAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGgccg                                  >  7:118134/1‑149 (MQ=255)                                                                                                     aaaGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGgccgcc                                >  3:31141/1‑149 (MQ=255)                                                                                                     aaaGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGAGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGCCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCCgccgcc                                >  5:206963/1‑149 (MQ=255)                                                                                                      aaGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAAGTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGGACACATGGACCTGCTAGACCCGgccg                                  >  8:263070/1‑146 (MQ=255)                                                                                                              aGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACAt                                                      <  8:181143/118‑1 (MQ=255)                                                                                                              aGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACAt                                                      >  7:181143/1‑118 (MQ=255)                                                                                                                cTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCAt                     <  7:128966/149‑1 (MQ=255)                                                                                                                    ttCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGc                 <  4:17539/149‑1 (MQ=255)                                                                                                                       g‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGc                 <  6:121012/146‑1 (MQ=255)                                                                                                                       g‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGc                 >  5:121012/1‑146 (MQ=255)                                                                                                                       g‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCa               >  3:258372/1‑148 (MQ=255)                                                                                                                       g‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCa               >  1:65922/1‑148 (MQ=255)                                                                                                                         ccGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCata             <  8:118134/149‑1 (MQ=255)                                                                                                                         cactctcCTGGAGCACGGCGCGCAGCTGCGCGTCGTCAAGACCCGGGGTGGTGCAGAATTAAACACCCCGCTGACCCTGCAGGCGCTCTCCCGCCACCCCGAACACATGGACCTGCTCGACCCGGCAGCCGAACCTGCCATGGGCCata             >  6:150123/4‑149 (MQ=255)                                                                                                                                cTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGa         <  3:154783/146‑1 (MQ=255)                                                                                                                                cTGGAGCACGGCGCGCAGGTGCGCGGCGTCATGACCCGGGGTTGTGCAGATTTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCt      <  6:206963/149‑1 (MQ=255)                                                                                                                                   gAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGc   <  2:299987/149‑1 (MQ=255)                                                                                                                                   gAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGc   >  6:174357/1‑149 (MQ=255)                                                                                                                                    atgacgacgCGCACCTGCGCGCCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGc   >  3:112219/4‑148 (MQ=255)                                                                                                                                      cACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGCCCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGcc  >  8:308303/1‑147 (MQ=255)                                                                                                                                      cACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGcc  >  4:219504/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                       AAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGCC  >  NC_002947/6033639‑6033917

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑