Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,364,421	C→G	5.2%	P71P (CCG→CCC)	folP ←	dihydropteroate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,364,421	0	C	G	5.2%	282.1 / 3.1	89	P71P (CCG→CCC)	folP	dihydropteroate synthase
Reads supporting (aligned to +/- strand):  ref base C (46/38);  new base G (2/3);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TCAGCCCGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGG  >  NC_002947/5364273‑5364568                                                                                                                                                     |                                                                                                                                                      tCAGCCCGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGCCCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCcccc                                                                                                                                                       >  3:193955/1‑149 (MQ=255) tCAGCCCGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGAAGGCCTCGACCATTGGCGCCACACCCACCAGCTCCTCGGTCACAGACACTGCACGCGCAccc                                                                                                                                                       >  6:372718/1‑149 (MQ=255)     cccGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCCCCCGGgc                                                                                                                                                   >  8:185566/1‑149 (MQ=255)       cGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGgcgc                                                                                                                                                 <  2:352806/149‑1 (MQ=255)         gCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCCCGGGTGCGCGt                                                                                                                                               <  2:254177/149‑1 (MQ=255)         gCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGt                                                                                                                                               <  4:122084/149‑1 (MQ=255)           ccccAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTg                                                                                                                                              >  6:65906/1‑148 (MQ=255)               aGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTc                                                                                                                                         <  1:163457/149‑1 (MQ=255)                 gcgcgcGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGc                                                                                                                                       <  8:90573/149‑1 (MQ=255)                 gcgcgcGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGc                                                                                                                                       >  5:8026/1‑149 (MQ=255)                   gcgcgGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGt                                                                                                                                               >  4:67171/1‑139 (MQ=255)                   gcgcgGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGt                                                                                                                                               <  3:67171/139‑1 (MQ=255)                      cgGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCg                                                                                                                                  <  2:59924/149‑1 (MQ=255)                       ggCTTCTCGCATGAGAGCAGGCGTGGAGGTATCGACCGAAATGACCAGGTCGTGGCGACTGTTGATGGCCTGGACAATTGGCGCCACACGCTCCTGCTCCTCTGTCACGGAGTCTGCACGCGCACCGGGGCGAGTCCACGCGcccgggt                                                                                                                                 <  1:253305/149‑6 (MQ=255)                            cACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATc                                                                                                                             <  7:263007/148‑1 (MQ=255)                             aCGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCg                                                                                                                            <  2:333735/148‑1 (MQ=255)                                 aTGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCCCACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCTATC‑ag                                                                                                                      >  4:271885/1‑148 (MQ=255)                                  tGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AAc                                                                                                                     <  4:187518/149‑1 (MQ=255)                                      agcagGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCAAGCTCCTCGGTAACCGAGACTGCACGCGCACCCGGGCGCGTGGCAGCGCCACCGATAACGAAC‑acaattg                                                                                                                 >  8:53240/1‑143 (MQ=255)                                         agGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGTGCGTGGACTCGCCACCGATATCGATC‑AACGTTgc‑gc                                                                                                             >  8:99943/1‑149 (MQ=255)                                         agGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTgc‑gc                                                                                                             >  7:323944/1‑149 (MQ=255)                                          ggCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑Gcc                                                                                                            <  4:11756/149‑1 (MQ=255)                                               ggAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCcggcgg                                                                                                       >  6:193672/1‑149 (MQ=255)                                                   gTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAAc                                                                                                   >  5:26948/1‑149 (MQ=255)                                                        gACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCg                                                                                              >  1:177567/1‑149 (MQ=255)                                                        gACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCg                                                                                              >  2:380470/1‑149 (MQ=255)                                                         aCCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATTGATC‑AACGTTGC‑GCCGGCGGCAACCATCGc                                                                                             >  3:291634/1‑149 (MQ=255)                                                         aCCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGc                                                                                             >  2:185536/1‑149 (MQ=255)                                                           cGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCtt                                                                                           >  2:9609/1‑149 (MQ=255)                                                           cGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCtt                                                                                           >  5:198046/1‑149 (MQ=255)                                                             aaaTGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTTGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCg                                                                                         >  2:335104/1‑149 (MQ=255)                                                                 gACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCCCCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGt                                                                                     >  3:358114/1‑149 (MQ=255)                                                                  aCCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCCCGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAAACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCg                                                                                      <  1:44886/148‑1 (MQ=255)                                                                  aCCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGt                                                                                     <  3:345232/148‑1 (MQ=255)                                                                    cACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGgc                                                                                  >  1:183178/1‑149 (MQ=255)                                                                        tCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAg                                                                              >  6:150636/1‑149 (MQ=255)                                                                         cGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGa                                                                                                                                            <  5:182208/88‑1 (MQ=255)                                                                         cGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGa                                                                                                                                            >  6:182208/1‑88 (MQ=255)                                                                           aGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGc                                                                           <  5:372718/149‑1 (MQ=255)                                                                             gCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCg                                                                                      >  4:98354/1‑136 (MQ=255)                                                                             gCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCg                                                                                      <  3:98354/136‑1 (MQ=255)                                                                              cGACGGTTTATGGGCTCGACGATTGGCTCCACACGCTCTAGCTCGTAGGTCACAGCGACTGCACGCGGACGCGGGTGCGTGGAGTCGCCACCGATTTCTTTC‑AACGTTGC‑GCCGGAGGCAACCATCGCTTCGGCGTGGGGCAGGGGCtc                                                                        <  6:65096/149‑1 (MQ=255)                                                                               gACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCtcg                                                                       >  3:298233/1‑149 (MQ=255)                                                                               gACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCtcg                                                                       >  7:194214/1‑149 (MQ=255)                                                                                  tGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCtcgtcg                                                                    >  7:181062/1‑149 (MQ=255)                                                                                  tGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTAACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACAATCGCTACGGCGTGGCGCAGGGCAtcgtcg                                                                    >  1:32533/1‑149 (MQ=255)                                                                                    ttGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTcgcg                                                                  >  5:368744/1‑149 (MQ=255)                                                                                    ttGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTcgcg                                                                  >  2:168292/1‑149 (MQ=255)                                                                                     tGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCgcgc                                                                 >  1:139153/1‑149 (MQ=255)                                                                                        tGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGctga                                                              <  2:160824/149‑1 (MQ=255)                                                                                        tGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGctga                                                              <  2:115853/149‑1 (MQ=255)                                                                                          gCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGctgact                                                            >  4:331810/1‑149 (MQ=255)                                                                                           ccTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGctgactg                                                           <  4:298233/149‑1 (MQ=255)                                                                                            cTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTg                                                                                    >  4:333711/1‑123 (MQ=255)                                                                                            cTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTg                                                                                    <  3:333711/123‑1 (MQ=255)                                                                                            cTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGctgactg                                                           >  7:146394/1‑148 (MQ=255)                                                                                            cTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGctgactg                                                           >  5:253634/1‑148 (MQ=255)                                                                                               gACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGaa                                                         <  4:172898/147‑1 (MQ=255)                                                                                               gACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGc                                                       <  2:177567/149‑1 (MQ=255)                                                                                                         gcCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATc                                             >  3:167512/1‑149 (MQ=255)                                                                                                         gcCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATc                                             >  6:99422/1‑149 (MQ=255)                                                                                                          ccACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCg                                            <  6:26948/149‑1 (MQ=255)                                                                                                          ccACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCg                                            <  4:7868/149‑1 (MQ=255)                                                                                                          ccACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCg                                            >  3:7868/1‑149 (MQ=255)                                                                                                           cacaCGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCg                                            >  1:92911/1‑148 (MQ=255)                                                                                                            acacGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGga                                          <  3:89753/149‑1 (MQ=255)                                                                                                             cacGCTCCAGCTCCTCGGTCACCGAGACTGCACGGGCAGCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGga                                          <  3:197302/148‑1 (MQ=255)                                                                                                              acGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGgaga                                        <  1:119738/149‑1 (MQ=255)                                                                                                              acGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGgaga                                        <  4:348767/149‑1 (MQ=255)                                                                                                              acGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGgag                                         >  5:210656/1‑148 (MQ=255)                                                                                                              acGCTCCAGCTCCTCGGTCACCGAGACAGCACGCCCACCCGGCCGCGTGGACACCCCACCCATATCGATC‑AACCTTTC‑CACGGCGGCAACCAACCCATCGGCGTGGCGCAGGGCCTCGTCGCTCTTACATAAGCGACCGCCATCCga                                          >  4:197302/1‑147 (MQ=255)                                                                                                                       ctcctcGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCgg                               >  3:170797/1‑149 (MQ=255)                                                                                                                       ctcctcGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCgg                               >  6:337704/1‑149 (MQ=255)                                                                                                                          ctcGGTCACCGAGACTGCACGCGCACCCGGGCGCGTTGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCgggg                             <  1:335104/148‑1 (MQ=255)                                                                                                                          ctcGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCggggg                            <  6:8026/149‑1 (MQ=255)                                                                                                                          ctcGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCggggg                            <  4:193955/149‑1 (MQ=255)                                                                                                                          ctcGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCggg                              <  3:181150/147‑1 (MQ=255)                                                                                                                             ggTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGt                           <  5:99422/147‑1 (MQ=255)                                                                                                                               tCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGata                       >  7:99992/1‑149 (MQ=255)                                                                                                                                 aCCGCGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATCTCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTGCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATAtt                     <  8:16048/149‑1 (MQ=255)                                                                                                                                 aCCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATAtt                     >  3:168945/1‑149 (MQ=255)                                                                                                                                   cGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTga                   <  1:185536/149‑1 (MQ=255)                                                                                                                                   cGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTga                   <  4:170797/149‑1 (MQ=255)                                                                                                                                      gagtccACGCGCCCGGGTGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTgaga                 >  6:325967/6‑148 (MQ=255)                                                                                                                                      gaCTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGaa                <  2:139153/149‑1 (MQ=255)                                                                                                                                         tccACGCGCCCGGGTGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGaa                                                         >  1:26468/3‑105 (MQ=255)                                                                                                                                         tccACGCGCCCGGGTGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGaa                                                         <  2:26468/103‑1 (MQ=255)                                                                                                                                          gCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATAcc            >  4:12926/1‑149 (MQ=255)                                                                                                                                            aCGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCa          <  6:368744/149‑1 (MQ=255)                                                                                                                                              gcgcACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGCAGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCAt         >  1:83384/1‑149 (MQ=255)                                                                                                                                              gcgcACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTCGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGCCCGCCATCGGAGCAGGAATCGGGGGTGATATTGAGAATACCCATg        >  1:181472/1‑149 (MQ=255)                                                                                                                                                   cccGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATg   <  7:53240/149‑1 (MQ=255)                                                                                                                                                    ccGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTga                   <  1:191265/132‑1 (MQ=255)                                                                                                                                                    ccGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTga                   >  2:191265/1‑132 (MQ=255)                                                                                                                                                    ccGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATgg  <  1:312664/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TCAGCCCGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGG  >  NC_002947/5364273‑5364568

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