Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,148,408	C→A	5.3%	P116T (CCG→ACG)	msuE →	FMN reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,148,408	0	C	A	5.3%	244.6 / 3.4	75	P116T (CCG→ACG)	msuE	FMN reductase
Reads supporting (aligned to +/- strand):  ref base C (34/37);  new base A (2/2);  total (36/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 = 5.16e-01

ATCACGCACCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACC  >  NC_002947/3148260‑3148554                                                                                                                                                     |                                                                                                                                                    aTCACGCATCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAAc                                                                                                                                                     <  5:100439/149‑1 (MQ=255)      gCACCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCAt                                                                                                                                                >  3:376886/1‑149 (MQ=255)        aCCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCAttt                                                                                                                                              <  3:283244/149‑1 (MQ=255)        aCCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCAt                                                                                                                                                >  1:363374/1‑147 (MQ=255)        aCCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCAt                                                                                                                                                <  2:363374/147‑1 (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                                                                                                                                 <  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                                                                                    <  6:196075/149‑1 (MQ=255)                                                                 gAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTc                                                                                    >  3:129178/1‑149 (MQ=255)                                                                   ggTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTcc                                                                                   >  5:255642/1‑148 (MQ=255)                                                                   ggTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTcc                                                                                   >  5:1899/1‑148 (MQ=255)                                                                   ggTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTcc                                                                                   >  5:7259/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                                                  >  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:7259/149‑1 (MQ=255)                                                                                                                                    ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                 <  6:255642/149‑1 (MQ=255)                                                                                                                                    ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                 <  6:1899/149‑1 (MQ=255)                                                                                                                                    ccTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtg                 >  4:4790/1‑149 (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                                                     >  2:29342/4‑106 (MQ=255)                                                                                                                                           tggccTGCAAATGCGGTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGAc                                                     <  1:29342/103‑1 (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)                                                                                                                                                     |                                                                                                                                                    ATCACGCACCCATGGTCTGCTGCAGGCCCTGGTGGAGCGGCTGCCGGCCGTGTTGCCGAAGCTGGAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCA‑TCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACC  >  NC_002947/3148260‑3148554

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