Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,148,412	A→G	5.2%	H117R (CAT→CGT) ‡	msuE →	FMN reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,148,412	0	A	G	5.2%	234.8 / 2.3	77	H117R (CAT→CGT) ‡	msuE	FMN reductase
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base G (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

GCACCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTG  >  NC_002947/3148265‑3148558                                                                                                                                                    |                                                                                                                                                    gCACCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCAt                                                                                                                                                    >  3:376886/1‑149 (MQ=255)   aCCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCAttt                                                                                                                                                  <  3:283244/149‑1 (MQ=255)   aCCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCAt                                                                                                                                                    <  2:363374/147‑1 (MQ=255)   aCCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCAt                                                                                                                                                    >  1:363374/1‑147 (MQ=255)     ccATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGc                                                                                                                                                <  3:295192/149‑1 (MQ=255)           tCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCa                                                                                                                                          >  7:159025/1‑149 (MQ=255)           tCTGCTGCAGGCCCTGGGGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAAATGCGGTTGCAGGTCa                                                                                                                                          <  1:183727/149‑1 (MQ=255)              gctgcAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCg                                                                                                                                      <  2:299823/149‑1 (MQ=255)               ctgcAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑Tc                                                                                                                                       <  5:157320/147‑1 (MQ=255)                 gcAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAg                                                                                                                                    >  7:65856/1‑148 (MQ=255)                  cAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGGCCTGCAACCGCATTTGCAGGCCA‑TCGAg                                                                                                                                    <  6:248669/147‑1 (MQ=255)                  cAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGc                                                                                                                                   <  6:355322/148‑1 (MQ=255)                   aGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCaa                                                                                                                                 <  8:159025/149‑1 (MQ=255)                       ccTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGcc                                                                                                                              >  7:180067/1‑148 (MQ=255)                           gtggAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGAcc                                                                                                                          <  4:246852/148‑1 (MQ=255)                           gtggAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGAcc                                                                                                                          <  4:241999/148‑1 (MQ=255)                           gtggAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGAc                                                                                                                           <  4:77152/147‑1 (MQ=255)                                gCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTAc                                                                                                                    >  2:276058/1‑149 (MQ=255)                                 cGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGa                                                                                                                                     >  6:226850/1‑131 (MQ=255)                                 cGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGa                                                                                                                                     <  5:226850/131‑1 (MQ=255)                                  ggCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTg                                                                                                                  <  1:365880/149‑1 (MQ=255)                                    cTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGt                                                                                                                <  3:103239/149‑1 (MQ=255)                                         ggccgTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCa                                                                                                           <  2:375161/149‑1 (MQ=255)                                              tgtTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCg                                                                                                      <  6:3302/149‑1 (MQ=255)                                              tgtTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCGGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCg                                                                                                      <  5:124968/149‑1 (MQ=255)                                                   ccGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGgtgtg                                                                                                 >  1:326444/1‑148 (MQ=255)                                                           ggAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGc                                                                                                                                                >  4:19976/1‑95 (MQ=255)                                                            gagggtcACTGTGTGCGCATCGCCGAACAGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTc                                                                                        <  8:175845/143‑1 (MQ=255)                                                            gAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTc                                                                                        >  3:129178/1‑149 (MQ=255)                                                            gAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTc                                                                                        <  6:196075/149‑1 (MQ=255)                                                              ggTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTcc                                                                                       >  5:7259/1‑148 (MQ=255)                                                              ggTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTcc                                                                                       >  5:255642/1‑148 (MQ=255)                                                              ggTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTcc                                                                                       >  5:1899/1‑148 (MQ=255)                                                               gTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTa                                                                                     <  2:382150/149‑1 (MQ=255)                                                               gTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTa                                                                                     >  4:102885/1‑149 (MQ=255)                                                                   aCTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACAcc                                                                                 <  7:169823/149‑1 (MQ=255)                                                                     tGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCgg                                                                               >  6:8832/1‑149 (MQ=255)                                                                        gTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGttt                                                                            <  6:2562/149‑1 (MQ=255)                                                                          gcgcATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGt                                                                          >  3:165471/1‑149 (MQ=255)                                                                              aTCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTc                                                                        >  4:197035/1‑147 (MQ=255)                                                                              aTCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACACGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCaa                                                                      <  1:276058/149‑1 (MQ=255)                                                                               tCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAg                                                                     >  6:287445/1‑149 (MQ=255)                                                                                 gCCGCACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCa                                                                   >  8:107592/1‑149 (MQ=255)                                                                                  ccGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCAc                                                                  >  3:140189/1‑149 (MQ=255)                                                                                   cGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCAcc                                                                 >  3:292942/1‑149 (MQ=255)                                                                                    gAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCt                                                                >  2:386328/1‑149 (MQ=255)                                                                                        tGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTc                                                            >  6:82927/1‑149 (MQ=255)                                                                                              ccAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTg                                                      >  5:44736/1‑149 (MQ=255)                                                                                              ccAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTg                                                      >  1:157146/1‑149 (MQ=255)                                                                                                    tGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGAt                                                >  2:104710/1‑149 (MQ=255)                                                                                                    tGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGAt                                                >  4:252207/1‑149 (MQ=255)                                                                                                      gCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCa                                              <  4:376886/149‑1 (MQ=255)                                                                                                       ccGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCAc                                             <  4:129178/149‑1 (MQ=255)                                                                                                              gCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTTTCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCg                                       <  2:326444/148‑1 (MQ=255)                                                                                                              gCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGc                                      >  8:124887/1‑149 (MQ=255)                                                                                                                  gAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTc                                    >  8:268690/1‑147 (MQ=255)                                                                                                                    gcgcGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAgg                                <  5:82927/149‑1 (MQ=255)                                                                                                                      gcgATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGc                              <  4:165471/149‑1 (MQ=255)                                                                                                                        gATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGt                            >  6:3433/1‑149 (MQ=255)                                                                                                                          tAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGc                          <  5:8832/149‑1 (MQ=255)                                                                                                                            gcgcCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGccc                        <  5:287445/149‑1 (MQ=255)                                                                                                                               ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                     <  6:1899/149‑1 (MQ=255)                                                                                                                               ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                     <  6:255642/149‑1 (MQ=255)                                                                                                                               ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                     >  4:4790/1‑149 (MQ=255)                                                                                                                               ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                     <  3:197035/149‑1 (MQ=255)                                                                                                                               ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                     <  6:7259/149‑1 (MQ=255)                                                                                                                                  tcgatggccTGCAAATGCGGTTGCAGGTCAGCCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtggt                   >  3:232351/8‑149 (MQ=38)                                                                                                                                     cTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCtt               >  1:357128/1‑149 (MQ=255)                                                                                                                                      tggccTGCAAATGCGGTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGAc                                                         <  1:29342/103‑1 (MQ=255)                                                                                                                                      tggccTGCAAATGCGGTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGAc                                                         >  2:29342/4‑106 (MQ=255)                                                                                                                                          cTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTgccgc          >  5:290610/1‑149 (MQ=255)                                                                                                                                           tGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTc                                                                        >  7:289136/1‑86 (MQ=255)                                                                                                                                           tGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTc                                                                        <  8:289136/86‑1 (MQ=255)                                                                                                                                             cAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCAc       <  8:180067/149‑1 (MQ=255)                                                                                                                                              aaCCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCAcc      <  7:268690/149‑1 (MQ=255)                                                                                                                                                  gCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCggtg  <  2:227520/149‑1 (MQ=255)                                                                                                                                                  gCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCggtg  <  6:290610/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCACCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTG  >  NC_002947/3148265‑3148558

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