Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,143,479	A→T	5.2%	T87S (ACC→TCC)	PP_1899 →	ExbD/TolR family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,143,479	0	A	T	5.2%	232.0 / 2.3	77	T87S (ACC→TCC)	PP_1899	ExbD/TolR family protein
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base T (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 = 9.33e-01

CCGCGAAACGCAACTTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGA  >  NC_002947/2143331‑2143623                                                                                                                                                     |                                                                                                                                                  ccGCGAAACGCAACTTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCa                                                                                                                                                   >  2:156854/1‑149 (MQ=255)   gcgAAACGCAACTTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCAcc                                                                                                                                                 <  8:106216/149‑1 (MQ=255)   gcgAAACGCAACTTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCAcc                                                                                                                                                 >  4:80277/1‑149 (MQ=255)        aCGCAACTTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAg                                                                                                                                            >  7:26224/1‑149 (MQ=255)                 cgcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGcc                                                                                                                                     <  1:191507/147‑1 (MQ=255)                  gcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCAt                                                                                                                                   <  5:109143/148‑1 (MQ=255)                  gcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATa                                                                                                                                  <  1:60436/149‑1 (MQ=255)                       gAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgc                                                                                                                              <  5:26125/148‑1 (MQ=255)                        aGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                            <  1:156854/149‑1 (MQ=255)                          cTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                            >  2:216625/1‑147 (MQ=255)                          cTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGa                                                                                                                          <  3:197357/149‑1 (MQ=255)                           tGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCTCGCTCAGGGTGGcc                                                                                                                                     >  8:147726/1‑137 (MQ=255)                           tGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCTCGCTCAGGGTGGcc                                                                                                                                     <  7:147726/137‑1 (MQ=255)                            gCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGt                                                                                                                        >  3:104836/1‑149 (MQ=255)                             cccGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTc                                                                                                                       >  8:356946/1‑149 (MQ=255)                              ccGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCg                                                                                                                      >  6:93633/1‑149 (MQ=255)                                gAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCggg                                                                                                                    <  2:90227/149‑1 (MQ=255)                                  aGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGggcg                                                                                                                  <  4:242554/149‑1 (MQ=255)                                    ccaccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGggcggc                                                                                                                <  7:6450/149‑1 (MQ=255)                                       ccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Gac                                                                                                            <  6:381585/149‑1 (MQ=255)                                       ccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Ga                                                                                                             >  6:328486/1‑148 (MQ=255)                                       ccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Ga                                                                                                             >  6:334207/1‑148 (MQ=255)                                        caGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCggg                                                                                                                    <  6:161740/141‑1 (MQ=255)                                        caGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCggg                                                                                                                    >  5:161740/1‑141 (MQ=255)                                          gcgcCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Gacaa                                                                                                          >  5:333221/1‑148 (MQ=255)                                            gcCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGGGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Gacaac                                                                                                         >  1:336324/1‑147 (MQ=255)                                               gAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGtt                                                                                                    <  8:26224/149‑1 (MQ=255)                                                 ctgtgcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GAcaacaa                                                                                                       >  4:206454/5‑144 (MQ=255)                                                    ggcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGc                                                                                               >  2:315860/1‑149 (MQ=255)                                                    ggcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGc                                                                                               <  4:104836/149‑1 (MQ=255)                                                    ggcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGc                                                                                               >  4:130395/1‑149 (MQ=255)                                                     gcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCt                                                                                              >  6:216614/1‑149 (MQ=255)                                                      cgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTg                                                                                             <  5:93633/149‑1 (MQ=255)                                                      cgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTg                                                                                             <  1:315860/149‑1 (MQ=255)                                                      cgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTg                                                                                             <  3:80277/149‑1 (MQ=255)                                                       gccggccgACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCGTGAGCGAGGCCATTGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTgg                                                                                            <  2:178944/149‑1 (MQ=255)                                                            ccgACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATc                                                                                       <  1:216625/149‑1 (MQ=255)                                                             cgACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCa                                                                                      <  6:398347/149‑1 (MQ=255)                                                                 ccGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCAGGCCATCAgcgc                                                                                  <  7:356946/149‑1 (MQ=255)                                                                    ggCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGa                                                                               <  2:336324/149‑1 (MQ=255)                                                                        aGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGc                                                                           <  3:130395/149‑1 (MQ=255)                                                                        aGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGc                                                                           <  5:216614/149‑1 (MQ=255)                                                                           tggtggAAATCACCTTCAGCGCCGCAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAg                                                                        <  7:76099/149‑1 (MQ=255)                                                                               ggAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGAccc                                                                    >  1:357301/1‑149 (MQ=255)                                                                                gAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                            <  1:195909/93‑1 (MQ=255)                                                                                gAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                            >  2:195909/1‑93 (MQ=255)                                                                                  aaTCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGc                                                                 <  6:333221/149‑1 (MQ=255)                                                                                       ccATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAgcg                                                                                   >  8:219072/1‑126 (MQ=255)                                                                                       ccATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAgcg                                                                                   <  7:219072/126‑1 (MQ=255)                                                                                       ccATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAg                                                            >  2:213117/1‑149 (MQ=255)                                                                                        cATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGc                                                                                                                                      <  6:228272/75‑1 (MQ=255)                                                                                        cATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGc                                                                                                                                      >  5:228272/1‑75 (MQ=255)                                                                                        cATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAgg                                                           >  2:334457/1‑149 (MQ=255)                                                                                         aTCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCAGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAgg                                                           <  1:101671/149‑1 (MQ=255)                                                                                         aTCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGc                                                          >  4:186232/1‑149 (MQ=255)                                                                                                ccGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGt                                                    >  2:315051/1‑148 (MQ=255)                                                                                                ccGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTc                                                   <  1:213117/149‑1 (MQ=255)                                                                                                        gtgtATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAAt                                           >  6:162473/1‑149 (MQ=255)                                                                                                                 gTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCAGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGAtgctgc                                  >  2:129030/1‑149 (MQ=255)                                                                                                                   gAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCg                                >  6:187935/1‑149 (MQ=255)                                                                                                                   gAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCg                                >  6:194844/1‑149 (MQ=255)                                                                                                                       aacCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGc                              >  7:162558/1‑147 (MQ=255)                                                                                                                             ctgctgCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCgg                      >  5:15706/1‑149 (MQ=255)                                                                                                                                ctgcCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGttt                   >  7:316022/1‑149 (MQ=255)                                                                                                                                   cccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                 <  5:194844/148‑1 (MQ=255)                                                                                                                                   cccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                 <  6:15706/148‑1 (MQ=255)                                                                                                                                    ccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGc                                                                                           >  4:97459/1‑73 (MQ=255)                                                                                                                                    ccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGc                                                                                           <  3:97459/73‑1 (MQ=255)                                                                                                                                    ccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGc               >  7:34785/1‑149 (MQ=255)                                                                                                                                     cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGa                                        <  7:351045/123‑1 (MQ=255)                                                                                                                                     cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGa                                        >  8:351045/1‑123 (MQ=255)                                                                                                                                     cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGc               >  5:93459/1‑148 (MQ=255)                                                                                                                                      aaaaGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCa             <  6:93459/149‑1 (MQ=255)                                                                                                                                            tctatGGCCTCGCTCAGGGTGGCCATAGAGCGCGAGTCGGGCGGC‑GGCAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                 >  5:114946/5‑139 (MQ=255)                                                                                                                                            tctatGGCCTCGCTCAGGGTGGCCATAGAGCGCGAGTCGGGCGGC‑GGCAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                 <  6:114946/135‑1 (MQ=255)                                                                                                                                             aCCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTgcgc      <  1:334457/149‑1 (MQ=255)                                                                                                                                             aCCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTgcgc      <  5:187935/149‑1 (MQ=255)                                                                                                                                                  gCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGa  <  1:315051/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCGCGAAACGCAACTTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGA  >  NC_002947/2143331‑2143623

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑