Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	474,155	T→A	6.4%	Q277L (CAG→CTG)	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	474,155	0	T	A	6.4%	181.4 / 3.0	62	Q277L (CAG→CTG)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base A (2/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

CGATTCCAGTGCGCGCCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTA  >  NC_002947/474007‑474291                                                                                                                                                     |                                                                                                                                          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)                                                                                                                                                     |                                                                                                                                          CGATTCCAGTGCGCGCCAGGCGGCCTTGCGCCCGGCCTGGTCGCCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATC‑TCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTA  >  NC_002947/474007‑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: ‑