Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,148,413	T→G	5.5%	H117Q (CAT→CAG) ‡	msuE →	FMN reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,148,413	0	T	G	5.5%	226.9 / 4.5	73	H117Q (CAT→CAG) ‡	msuE	FMN reductase
Reads supporting (aligned to +/- strand):  ref base T (32/37);  new base G (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

GCACCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTA  >  NC_002947/3148265‑3148561                                                                                                                                                     |                                                                                                                                                      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:241999/148‑1 (MQ=255)                           gtggAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGAcc                                                                                                                             <  4:246852/148‑1 (MQ=255)                           gtggAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGAc                                                                                                                              <  4:77152/147‑1 (MQ=255)                                gCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTAc                                                                                                                       >  2:276058/1‑149 (MQ=255)                                 cGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGa                                                                                                                                        <  5:226850/131‑1 (MQ=255)                                 cGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGa                                                                                                                                        >  6:226850/1‑131 (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:1899/1‑148 (MQ=255)                                                              ggTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTcc                                                                                          >  5:255642/1‑148 (MQ=255)                                                               gTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTa                                                                                        >  4:102885/1‑149 (MQ=255)                                                               gTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTa                                                                                        <  2:382150/149‑1 (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                                                         >  1:157146/1‑149 (MQ=255)                                                                                              ccAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTg                                                         >  5:44736/1‑149 (MQ=255)                                                                                                    tGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGAt                                                   >  4:252207/1‑149 (MQ=255)                                                                                                    tGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGAt                                                   >  2:104710/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:7259/149‑1 (MQ=255)                                                                                                                               ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                        >  4:4790/1‑149 (MQ=255)                                                                                                                               ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                        <  6:255642/149‑1 (MQ=255)                                                                                                                               ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                        <  6:1899/149‑1 (MQ=255)                                                                                                                               ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                        <  3:197035/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                                                                           <  8:289136/86‑1 (MQ=255)                                                                                                                                           tGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTc                                                                           >  7:289136/1‑86 (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)                                                                                                                                                     tttGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCCCCGGTGGTa  <  3:102885/149‑1 (MQ=255)                                                                                                                                                     tttGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTa  <  2:357128/149‑1 (MQ=255)                                                                                                                                                     tttGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTa  <  4:292942/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCACCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTA  >  NC_002947/3148265‑3148561

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