Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,346	T→G	6.1%	T13P (ACG→CCG)	trmE ←	GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,346	0	T	G	6.1%	294.5 / 9.5	99	T13P (ACG→CCG)	trmE	GTPase
Reads supporting (aligned to +/- strand):  ref base T (46/47);  new base G (4/2);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e-01

GTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTA  >  NC_002947/6200‑6492                                                                                                                                                   |                                                                                                                                                    gTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTg                                                                                                                                                    <  2:163576/148‑1 (MQ=255)    cAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACATCGCCACGGCCCAGGGCCGGGGCgg                                                                                                                                                >  2:351183/1‑148 (MQ=255)    cAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑a                                                                                                                                               >  2:200774/1‑149 (MQ=255)         ccAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATggcg                                                                                                                                          <  1:58800/149‑1 (MQ=255)           aGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCACAGGGCCGGACAACCTGACAATTCCCACATCGCCACGGCCCAGGGCCGGGGCG‑GTgtggga                                                                                                                                        <  1:351183/149‑6 (MQ=255)                ttCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGc                                                                                                                                                  <  6:276478/135‑1 (MQ=255)                ttCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGc                                                                                                                                                  >  5:276478/1‑135 (MQ=255)                ttCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATgg                                                                                                                                   >  8:249247/1‑149 (MQ=255)                      gtcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCa                                                                                                                              <  6:48230/148‑1 (MQ=255)                         gcgGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGca                                                                                                                          <  4:35957/149‑1 (MQ=255)                               tGGACGGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGGCGGCCTTGGCAGCCAGAGGGGCGGACAACCTGACAATTCCCACTCCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAgtg                                                                                                                     <  5:380106/148‑1 (MQ=255)                               tGGACCGTAATGGGCATGCCGCGGGGTCAGGGTCCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAAATCCCACACCGCCCCGCCCCTGGGCCCTGGCG‑ATGGCCGCGATGGTTTAACGCACAgtg                                                                                                                     >  6:380106/1‑148 (MQ=255)                                 gACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGc                                                                                                                                                  <  1:357334/118‑1 (MQ=255)                                 gACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGc                                                                                                                                                  >  2:357334/1‑118 (MQ=255)                                  aCCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACACCCTGACAATTCCCACACCGCTCCGGCCCTGGGCCGTGGCG‑a                                                                                                                                               >  8:122598/1‑119 (MQ=255)                                   ccGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCAt                                                                                                                >  8:77428/1‑149 (MQ=255)                                   ccGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCAt                                                                                                                >  7:225973/1‑149 (MQ=255)                                     gTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGc                                                                                                              <  5:218125/149‑1 (MQ=255)                                      tAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTGCATGCt                                                                                                             >  6:66759/1‑149 (MQ=255)                                              aTGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGcc                                                                                                     >  7:159700/1‑149 (MQ=255)                                                gCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCctc                                                                                                   >  2:89778/1‑149 (MQ=255)                                                gCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGGGTTAATGCTAGAAAGCct                                                                                                    >  6:265043/1‑148 (MQ=255)                                                 ccGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCctc                                                                                                   >  2:29963/1‑148 (MQ=255)                                                  cgcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaaa                                                                                                        <  8:261215/142‑1 (MQ=255)                                                  cgcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaaa                                                                                                        >  7:261215/1‑142 (MQ=255)                                                    cgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGTATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                                <  1:89778/149‑1 (MQ=255)                                                        gTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGcaca                                                                                                                        >  5:85255/1‑120 (MQ=255)                                                        gTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGcaca                                                                                                                        <  6:85255/120‑1 (MQ=255)                                                        gTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACaa                                                                                           <  1:244139/149‑1 (MQ=255)                                                         tCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATg                                                                                                               <  6:286624/128‑1 (MQ=255)                                                         tCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATg                                                                                                               >  5:286624/1‑128 (MQ=255)                                                           aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATg                                                                                                                                    >  7:279679/1‑105 (MQ=255)                                                           aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATg                                                                                                                                    <  8:279679/105‑1 (MQ=255)                                                           aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCG‑GTgt                                                                                                                                            >  6:366961/1‑96 (MQ=38)                                                           aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCG‑GTgt                                                                                                                                            <  5:366961/97‑2 (MQ=38)                                                              gTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCct                                                                                                    <  6:211806/134‑1 (MQ=255)                                                              gTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCct                                                                                                    >  5:211806/1‑134 (MQ=255)                                                              gTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGa                                                                                     <  1:375241/149‑1 (MQ=255)                                                               tGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGAc                                                                                    <  8:159700/149‑1 (MQ=255)                                                                gCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACa                                                                                   <  8:349519/149‑1 (MQ=255)                                                                       gTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCTGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGc                                                                             >  2:178002/1‑148 (MQ=255)                                                                         gATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaa                                                                           >  4:297194/1‑148 (MQ=255)                                                                             aaCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACg                                                                      >  3:333581/1‑149 (MQ=255)                                                                             aaCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCCATGGTTTCACGTACAGTGTTCATGCTCGAAAGCCTCTACGAAAAAAACGACAGAAACCAAAAACg                                                                      >  6:345657/1‑149 (MQ=255)                                                                              aCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGc                                                                     <  6:16532/149‑1 (MQ=255)                                                                                aGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGccc                                                                   <  8:225973/149‑1 (MQ=255)                                                                                  gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTTGCG‑ATGGCGGCGATGGTTTCACGCACAgtgt                                                                                                                    >  2:247271/1‑98 (MQ=255)                                                                                  gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTTGCG‑ATGGCGGCGATGGTTTCACGCACAgtgt                                                                                                                    <  1:247271/98‑1 (MQ=255)                                                                                  gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCa                                                                 >  4:315827/1‑149 (MQ=255)                                                                                  gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCa                                                                 >  5:37084/1‑149 (MQ=255)                                                                                  gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCa                                                                 >  3:280251/1‑149 (MQ=255)                                                                                  gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTAATGCTCGAAAGCCTCGACGACAAAAACGACAGATAGCAAAAACGCCCCa                                                                 >  1:241299/1‑149 (MQ=255)                                                                                    cagccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCAc                                                                >  3:181726/3‑148 (MQ=255)                                                                                                ccAGAGGGCCGGACAACCTGACAATTCCCACCCCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtttt                                                   <  5:66759/149‑1 (MQ=255)                                                                                                ccAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGccc                                                                   <  1:218239/133‑1 (MQ=255)                                                                                                ccAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGccc                                                                   >  2:218239/1‑133 (MQ=255)                                                                                                ccAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtttt                                                   >  2:329517/1‑149 (MQ=255)                                                                                                  atatggCCGGACAAACTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACATTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATACCAAAAACGCCCCACTAGGGGGCGtttttt                                                 <  4:181726/145‑1 (MQ=255)                                                                                                     gggCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCTCAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTAtt                                              <  1:200774/149‑1 (MQ=255)                                                                                                      ggCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCa                                                                 >  5:83988/1‑129 (MQ=255)                                                                                                      ggCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCa                                                                 <  6:83988/129‑1 (MQ=255)                                                                                                        ccGGACAACCTGACAATTCCCACACCGCCCCGGCCCAGGGCCGTGGCG‑ATGGCGGCGATCGTTTCACGCACACTGTTCATGCTCGAAAGCCTCAACGACAAAAACGACAGATAGAAAAAACGCCCCACTAGGGGGCGTTTTTTAttaa                                            >  5:64372/1‑146 (MQ=255)                                                                                                          ggACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAg                                         <  1:29963/149‑1 (MQ=255)                                                                                                            gcAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGc                                       >  5:346620/2‑149 (MQ=255)                                                                                                              aaCCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTa                                     <  7:249247/149‑1 (MQ=255)                                                                                                                ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCa                                                                                                                              <  7:170584/58‑1 (MQ=255)                                                                                                                ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCa                                                                                                                              >  8:170584/1‑58 (MQ=255)                                                                                                                ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGc                                   >  4:93600/1‑149 (MQ=255)                                                                                                                 cTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCg                                  >  1:42971/1‑149 (MQ=255)                                                                                                                 cTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCg                                  >  4:108720/1‑149 (MQ=255)                                                                                                                  tGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                                    >  5:71269/1‑114 (MQ=255)                                                                                                                  tGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                                    <  6:71269/114‑1 (MQ=255)                                                                                                                  tGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCg                                  >  5:227810/1‑148 (MQ=255)                                                                                                                  tGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGt                                 >  7:225539/1‑149 (MQ=255)                                                                                                                    aCAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGt                                 <  5:265043/147‑1 (MQ=255)                                                                                                                        ttCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGc                           <  8:225539/149‑1 (MQ=255)                                                                                                                           ccACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc                         <  4:280251/148‑1 (MQ=255)                                                                                                                           ccACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCg                        >  8:71309/1‑149 (MQ=255)                                                                                                                            cacaCCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGt                       >  1:230607/1‑149 (MQ=255)                                                                                                                            cacaCCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGt                       <  6:227810/149‑1 (MQ=255)                                                                                                                             ccatcGCCACGGCCCAGGGCCGGGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGt                       >  1:139302/5‑148 (MQ=255)                                                                                                                               acCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc                      >  2:62228/1‑147 (MQ=255)                                                                                                                                tcGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTa                             <  1:278050/138‑1 (MQ=255)                                                                                                                                tcGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTa                             >  2:278050/2‑139 (MQ=255)                                                                                                                                ccGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTa                                                              >  6:282752/1‑106 (MQ=255)                                                                                                                                ccGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTa                                                              <  5:282752/106‑1 (MQ=255)                                                                                                                                ccGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGt                                 <  1:112214/135‑1 (MQ=255)                                                                                                                                ccGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGt                                 >  2:112214/1‑135 (MQ=255)                                                                                                                                ccGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag                   <  1:329517/149‑1 (MQ=255)                                                                                                                                ccGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag                   <  3:93600/149‑1 (MQ=255)                                                                                                                                   ccccGGCCCTGGGCTGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcag                <  1:178002/149‑1 (MQ=255)                                                                                                                                   ccacGGCCCAGGGCCGGGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcag                >  2:120321/4‑149 (MQ=255)                                                                                                                                    cccGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagc               <  6:37084/149‑1 (MQ=255)                                                                                                                                       ggCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagcagc            <  7:77428/149‑1 (MQ=255)                                                                                                                                         cccTGGGCCGTTGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCGTCTACGACAAAAACGACAGATAGCAAATTCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtt          <  5:304874/149‑1 (MQ=255)                                                                                                                                              ggCGGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtttttt      <  6:64372/148‑1 (MQ=255)                                                                                                                                              ggCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttttt       <  3:108720/147‑1 (MQ=255)                                                                                                                                               gCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTgg    <  6:346620/149‑1 (MQ=255)                                                                                                                                                 cGTGGCG‑ATGGCGGCGATGGTTTCACGTACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTATCCCGTCAAGCAGCAGCTTTTTTGGTa  <  5:345657/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCG‑ATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTA  >  NC_002947/6200‑6492

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