Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,364,420	C→G	5.3%	G72R (GGT→CGT)	folP ←	dihydropteroate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,364,420	0	C	G	5.3%	284.2 / 3.2	89	G72R (GGT→CGT)	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.77e-01

ATCAGCCCGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGG  >  NC_002947/5364272‑5364568                                                                                                                                                     |                                                                                                                                                       aTCAGCCCGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCAcc                                                                                                                                                        >  4:89753/1‑149 (MQ=255)  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                                                                                                                                               <  3:67171/139‑1 (MQ=255)                    gcgcgGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGt                                                                                                                                               >  4:67171/1‑139 (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                                                                                              >  2:380470/1‑149 (MQ=255)                                                         gACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCg                                                                                              >  1:177567/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                                                                                                                                            >  6:182208/1‑88 (MQ=255)                                                                          cGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGa                                                                                                                                            <  5:182208/88‑1 (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:115853/149‑1 (MQ=255)                                                                                         tGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGctga                                                              <  2:160824/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                                            >  3:7868/1‑149 (MQ=255)                                                                                                           ccACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCg                                            <  4:7868/149‑1 (MQ=255)                                                                                                           ccACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCg                                            <  6:26948/149‑1 (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                            <  4:193955/149‑1 (MQ=255)                                                                                                                           ctcGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCggggg                            <  6:8026/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                   <  4:170797/149‑1 (MQ=255)                                                                                                                                    cGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTga                   <  1:185536/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)                                                                                                                                                     |                                                                                                                                                       ATCAGCCCGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATC‑AACGTTGC‑GCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGG  >  NC_002947/5364272‑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: ‑