Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	474,152	T→A	6.4%	Q278L (CAA→CTA) ‡	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	474,152	0	T	A	6.4%	188.0 / 2.8	62	Q278L (CAA→CTA) ‡	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base A (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

TGTCGATTCCAGTGCGCGCCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTA  >  NC_002947/474004‑474291                                                                                                                                                     |                                                                                                                                             tGTCGATTCCAGTGCGCGCCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCtt                                                                                                                                              >  5:60136/1‑149 (MQ=255)   tCGATTCCAGTGCGCGCCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGc                                                                                                                                            >  1:66986/1‑149 (MQ=255)    cGATTCCAGTGCGCGCCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCt                                                                                                                                           >  4:127471/1‑149 (MQ=255)       ttCCAGTGCGCGCCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTggg                                                                                                                                        <  3:79573/149‑1 (MQ=255)            gTGCGCGCCAGGCGGCCGTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGGGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCCTTGCGCAGGCCTTGCTGGGCCAgg                                                                                                                                   <  2:164916/149‑1 (MQ=255)                gcgcCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCaa                                                                                                                               <  2:168120/149‑1 (MQ=255)                 cgcCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAAt                                                                                                                              <  6:329477/149‑1 (MQ=255)                  gcCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATg                                                                                                                             >  7:70177/1‑149 (MQ=255)                  gcCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCCCTGGCCCAGCAAGGCCTGCGCAATg                                                                                                                             >  7:139093/1‑149 (MQ=255)                  gcCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCCCTGGCCCAGCAAGGCCTGCGCAATg                                                                                                                             <  4:154800/149‑1 (MQ=255)                    cAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGAc                                                                                                                           >  6:263249/1‑149 (MQ=255)                    cAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCCCTGGCCCAGCAAGGCCTGCGCAATgcc                                                                                                                           <  8:139093/149‑3 (MQ=38)                         ggcCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCa                                                                                                                      >  4:108536/1‑149 (MQ=255)                           ccTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATg                                                                                                                    <  6:60136/149‑1 (MQ=255)                             ttGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATg                                                                                                                    >  4:192580/1‑147 (MQ=255)                                 gcCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGccc                                                                                                              >  4:122304/1‑149 (MQ=255)                                   ccGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTc                                                                                                            >  1:143174/1‑149 (MQ=255)                                    cGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTcg                                                                                                           <  2:356663/149‑1 (MQ=255)                                    cGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTcg                                                                                                           >  8:242413/1‑149 (MQ=255)                                     ggCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTcga                                                                                                          >  8:211876/1‑149 (MQ=255)                                      gCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCCGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTcgac                                                                                                         >  1:113035/1‑149 (MQ=255)                                           gTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTcgacgacga                                                                                                    <  3:192580/149‑1 (MQ=255)                                            tCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCTGTGACAAGGGTGGCCACGGCATTGCGCAGTCCTTGCTGGGCCAGGGAAATTACGTCCATGTAGCCCTCGACGACgat                                                                                                   >  8:135419/1‑149 (MQ=255)                                              gCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACgatg                                                                                                  <  5:124698/148‑1 (MQ=255)                                                     aaaCAGAACAGCACGCTGGGCACCACGCGGAACAGGAGCTTGAGGTGCTCTTCGCTGGTGGAGGTTCCAAGGGTGTCCACGGAATTGCGCAGTCCTTGCTGGGCAAGGGCAATGAAGTCCATGTAGCCCTCGACGACGATGATC‑tcgtc                                                                                         >  8:298197/1‑149 (MQ=255)                                                        cagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGA‑GATCATCGTCGTCGAg                                                                                      <  2:66986/149‑1 (MQ=255)                                                         agaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAgg                                                                                     <  3:127471/149‑1 (MQ=255)                                                         agaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAgg                                                                                     >  8:224638/1‑149 (MQ=255)                                                              agCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCCATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCg                                                                                >  6:266616/1‑149 (MQ=255)                                                              agCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCg                                                                                >  8:36071/1‑149 (MQ=255)                                                               gCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCgg                                                                               <  3:191675/149‑1 (MQ=255)                                                                       ggCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGc                                                                         <  2:113035/147‑1 (MQ=255)                                                                           ccacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTCGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgc                                                                   <  4:275483/149‑1 (MQ=255)                                                                            cacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTgcgcg                                                                  >  6:186194/1‑149 (MQ=255)                                                                              cgcgGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGcc                                                                <  7:211876/149‑1 (MQ=255)                                                                                cgGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCt                                                               <  5:145021/148‑1 (MQ=255)                                                                                cgGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCt                                                               >  6:145021/1‑148 (MQ=255)                                                                                   aaCAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTa                                                           <  5:186194/149‑1 (MQ=255)                                                                                          gcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCg                                                    >  8:21478/1‑149 (MQ=255)                                                                                              gAGGTGCTCTTCGCTTGTGGCGGTGCCAAGGGTGTCCAAGGCATTGCGCAGGCCTTGCTGGGCCATGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTAc                                                 <  7:135419/148‑1 (MQ=255)                                                                                              gAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTa                                                  <  5:263249/147‑1 (MQ=255)                                                                                                ggTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACTGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTcga                                                                                                          <  6:309850/90‑1 (MQ=255)                                                                                                ggTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTcga                                                                                                          >  5:309850/1‑90 (MQ=255)                                                                                                    ctctTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTa                                                           <  6:16021/132‑1 (MQ=255)                                                                                                    ctctTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTa                                                           >  5:16021/1‑132 (MQ=255)                                                                                                    ctctTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCt                                          <  7:21478/149‑1 (MQ=255)                                                                                                       ttCGCTGGTGGCGGTGCCAAGGGTGGCCTCGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGcc                                       >  1:204746/1‑149 (MQ=255)                                                                                                          gCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCt                                    >  3:80234/1‑149 (MQ=255)                                                                                                                    gTGCCAAGGGTGGCCACGGCATTGAGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCTTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCGCCTGACCCTTGTGGATAAc                          >  3:186222/1‑149 (MQ=255)                                                                                                                      gCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCg                        >  3:245268/1‑149 (MQ=255)                                                                                                                        cAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                      <  8:70177/149‑1 (MQ=255)                                                                                                                        cAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                      >  8:65968/1‑149 (MQ=255)                                                                                                                        cAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                      <  2:204746/149‑1 (MQ=255)                                                                                                                              tgtccACGGCATTGCGCAGGCCTTGCTTGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTca                <  4:186222/146‑2 (MQ=255)                                                                                                                                gCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCgg              <  7:36071/149‑1 (MQ=255)                                                                                                                                gCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCgg              <  5:266616/149‑1 (MQ=255)                                                                                                                                      tcATTGCCCTGGCCCAGCAAGGCCTGAGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGt         >  6:305520/2‑148 (MQ=255)                                                                                                                                      gCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGtt        <  7:242413/149‑1 (MQ=255)                                                                                                                                       cATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTc       <  4:245268/149‑1 (MQ=255)                                                                                                                                            cgcAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGt   <  3:108536/148‑1 (MQ=255)                                                                                                                                            cgcAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTa  >  2:71326/1‑149 (MQ=255)                                                                                                                                                 gCCTTGCTGGGCCAGGGCAATGACGTCCATTTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAAcc                         >  2:331464/1‑121 (MQ=255)                                                                                                                                                 gCCTTGCTGGGCCAGGGCAATGACGTCCATTTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAAcc                         <  1:331464/121‑1 (MQ=255)                                                                                                                                                    ttGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCg                        <  8:167015/119‑1 (MQ=255)                                                                                                                                                    ttGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCg                        >  7:167015/1‑119 (MQ=255)                                                                                                                                                     |                                                                                                                                             TGTCGATTCCAGTGCGCGCCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTA  >  NC_002947/474004‑474291

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