Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,920,273	C→T	8.0%	R498Q (CGG→CAG)	acs ←	acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,920,273	0	C	T	8.0%	175.0 / 4.3	62	R498Q (CGG→CAG)	acs	acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base T (2/3);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.96e-01

GCGCATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTCCTTGCT  >  NC_002947/3920125‑3920418                                                                                                                                                     |                                                                                                                                                   gcgcATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGcc                                                                                                                                                    <  4:94025/149‑1 (MQ=255)   gcATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATG‑CCCGGGTGCAGCCGCa                                                                                                                                                 <  2:193430/149‑1 (MQ=255)    cATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCa                                                                                                                                                 <  2:177012/149‑1 (MQ=255)         ggAATTTCTGGACCTTGCCCGTCACCGTCATCGGGAACTAGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTTCCCTCGCCCACTCACGCAACTCTACCTGTCTGACAGCATGCCCCGGGTGCGGCTGCAcccgg                                                                                                                                            <  6:55958/149‑3 (MQ=255)         ggAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCTGGGTGCGGCTGCAcccgg                                                                                                                                            >  8:302807/1‑147 (MQ=255)            aTTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCTGGGTGCGGCTGCACCCggggc                                                                                                                                         <  7:302807/149‑3 (MQ=255)            aTTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGcc                                                                                                                                         <  6:315749/149‑1 (MQ=255)                  gCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCACGGGTGCAGCCGCACCCAgcc                                                                                                                                          >  1:166161/1‑140 (MQ=255)                  gCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCACGGGTGCAGCCGCACCCAgcc                                                                                                                                          <  2:166161/142‑3 (MQ=255)                     ccTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATctct                                                                                                                                <  7:78154/149‑1 (MQ=255)                      cTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCtt                                                                                                                               >  7:313867/1‑149 (MQ=255)                         gCCCGTCACCGTCATCGGGATCTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCa                                                                                                                             <  4:292101/148‑1 (MQ=255)                                 ccgtccTCGGGAACTCGTCGACAAATCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTtg                                                                                                                    <  6:34343/149‑1 (MQ=255)                                     caTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTtgctg                                                                                                                 >  2:73210/1‑148 (MQ=255)                                      aTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCCCTGACAGCATGCCCCGGGTGCAGCCGCCCCCCGGCCACGATCTCTTCACCGTACTtgctg                                                                                                                 >  4:97701/1‑147 (MQ=255)                                         gggAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAggg                                                                                                            >  5:343197/1‑149 (MQ=255)                                                   cgACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATc                                                                                                  <  8:250387/149‑1 (MQ=255)                                                       aaaaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCTTACTTGCTGCAGGGCACGCCAATCACCt                                                                                              >  6:288631/1‑149 (MQ=255)                                                       aaaaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCt                                                                                              <  5:295455/149‑1 (MQ=255)                                                       aaaaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCt                                                                                              >  6:5234/1‑149 (MQ=255)                                                       aaaaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCt                                                                                              >  2:67425/1‑149 (MQ=255)                                                       aaaaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCt                                                                                              >  2:197640/1‑149 (MQ=255)                                                            ggAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGcac                                                                                          >  8:238838/1‑148 (MQ=255)                                                                       agCACCTCGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCAccgc                                                                             >  6:331159/2‑149 (MQ=255)                                                                          aCCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCgg                                                                          <  1:197640/149‑1 (MQ=255)                                                                            cTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCgggt                                                                        >  3:130820/1‑149 (MQ=255)                                                                               gAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCgggtgg                                                                      <  5:5234/148‑1 (MQ=255)                                                                               gAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCgggtgg                                                                      <  5:331159/148‑1 (MQ=255)                                                                               gAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCgggtgg                                                                      <  8:313867/148‑1 (MQ=255)                                                                                    gcgcAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGggtgggtg                                                                   >  1:159030/1‑146 (MQ=255)                                                                                    gcgcAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGggtgggtg                                                                   >  1:146953/1‑146 (MQ=255)                                                                                      gcAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTgaagaa                                                              >  2:200591/1‑149 (MQ=255)                                                                                      gcAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTgaag                                                                >  6:190649/1‑147 (MQ=255)                                                                                      gcAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTgaag                                                                >  3:283610/1‑147 (MQ=255)                                                                                         aTCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTgaagaagaa                                                           <  5:84153/149‑1 (MQ=255)                                                                                          tCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAAc                                                          <  2:159030/149‑1 (MQ=255)                                                                                             ttGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAActc                                                        >  4:196829/1‑148 (MQ=255)                                                                                               gCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAActcct                                                      <  6:343197/148‑1 (MQ=255)                                                                                                         aCTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTcgcgc                                           <  4:130820/149‑1 (MQ=255)                                                                                                          cTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCAccgc                                                                             >  3:117078/1‑114 (MQ=255)                                                                                                          cTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCAccgc                                                                             <  4:117078/114‑1 (MQ=255)                                                                                                             aCGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCgg                                         <  1:73210/147‑1 (MQ=255)                                                                                                               gCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTa                                      >  8:271873/1‑148 (MQ=255)                                                                                                               gCAACTCTACCTCACTGACAGCATGCCCCGAGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTa                                      >  8:126675/1‑148 (MQ=255)                                                                                                                    tctACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCCACCGCCGGGTGGGTGAAGAAGAACtcctcc                                                    <  1:203356/130‑1 (MQ=255)                                                                                                                    tctACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACtcctcc                                                    >  2:203356/1‑129 (MQ=255)                                                                                                                          tCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcgcc                          >  8:191602/1‑149 (MQ=255)                                                                                                                              tGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTgaagaagaa                                                           <  6:54239/112‑1 (MQ=255)                                                                                                                                aCAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcgccgccgc                     >  5:287137/1‑148 (MQ=255)                                                                                                                                  aGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGaa                  <  1:67425/149‑1 (MQ=255)                                                                                                                                  aGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCACGCGGGTAAATGTTCTCGCCGCCGCGaa                  >  5:29407/1‑149 (MQ=255)                                                                                                                                   gCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAAt                 >  3:273549/1‑149 (MQ=255)                                                                                                                                      tGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGAt              <  3:196829/149‑1 (MQ=255)                                                                                                                                       gCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTtc                               >  2:198752/1‑131 (MQ=255)                                                                                                                                       gCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTtc                               <  1:198752/131‑1 (MQ=255)                                                                                                                                        ccccGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATc             <  7:238838/148‑1 (MQ=255)                                                                                                                                          ctgggTGCGGCTGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATc             >  1:282992/3‑146 (MQ=255)                                                                                                                                          ctgggTGCGGCTGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATc             <  2:282992/144‑1 (MQ=255)                                                                                                                                          ccGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATg          >  5:81240/1‑149 (MQ=255)                                                                                                                                          ccGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATg          >  5:325078/1‑149 (MQ=255)                                                                                                                                            gggTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTc        <  1:200591/149‑1 (MQ=255)                                                                                                                                                  aGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTCCTtgct  <  5:190649/149‑1 (MQ=255)                                                                                                                                                  aGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTCCTtgct  >  8:43686/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCGCATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGGGCACGCCAATCACCTGCACATCGG‑CCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGTCCTTGCT  >  NC_002947/3920125‑3920418

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