Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	563,005	C→A	5.4%	P30T (CCT→ACT) ‡	rpoA →	DNA‑directed RNA polymerase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	563,005	0	C	A	5.4%	240.3 / 3.1	74	P30T (CCT→ACT) ‡	rpoA	DNA‑directed RNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base A (2/2);  total (35/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 = 8.01e-01

TCAACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTAA  >  NC_002947/562858‑563152                                                                                                                                                     |                                                                                                                                                    tCAACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGcc                                                                                                                                                    >  1:158959/1‑149 (MQ=255)  cAACGAAAACCTGATTGTCGAGCTCTACTCGAAGTAAGGGCAAGAAAATAGGTGCATCAATGCAGATTTCGGTAAATGAGTACCTGACTCCCCTACACAT‑TGATGTGCACGAAGTCACTCCGACCCGCGCCAAAATTACGCACGAGcc                                                                                                                                                    >  6:222323/1‑148 (MQ=255)   aaCGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCct                                                                                                                                                   <  1:275073/148‑1 (MQ=255)    aCGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCctct                                                                                                                                                 >  1:288515/1‑149 (MQ=255)       aaaaCCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGa                                                                                                                                              <  6:203492/149‑1 (MQ=255)         aaCCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGAGGCTCGAGc                                                                                                                                            <  7:86147/149‑1 (MQ=255)            cTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                         <  6:359257/149‑1 (MQ=255)            cTGATTGGCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGATCCTGACTCCCCGCCACAT‑TGATCTGCAGGAAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                         <  7:295977/149‑1 (MQ=255)              gATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGc                                                                                                                                                     <  6:152135/135‑1 (MQ=255)              gATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGc                                                                                                                                                     >  5:152135/1‑135 (MQ=255)               aTTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGAGGCTCGAGCGtaatt                                                                                                                                      <  7:289014/149‑5 (MQ=255)                ttGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGttt                                                                                                                                     <  6:101634/149‑1 (MQ=255)                   tCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAG‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTCGTTTCgg                                                                                                                                  >  6:86700/1‑149 (MQ=255)                     gAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGc                                                                                                                                 <  5:266003/148‑1 (MQ=255)                         tctACTCCAAGTAAGGTCTAGAAAATAGGTGCATCCATGCAGATTTCTGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                         <  1:123147/136‑1 (MQ=255)                         tctACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATAc                                                                                                                            >  3:232876/1‑149 (MQ=255)                         tctACTACAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTTCAGGTAGTCAGTCCGACAAGCGCCAAAATTACGATAGAGCCTCTCGAGCGTg                                                                                                                                         >  2:123147/1‑136 (MQ=255)                             cTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTg                                                                                                                        <  5:387292/149‑1 (MQ=255)                               ccAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTggg                                                                                                                      >  3:243046/1‑149 (MQ=255)                                cAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGtt                                                                                                                                      >  5:202687/1‑132 (MQ=255)                                cAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGtt                                                                                                                                      <  6:202687/132‑1 (MQ=255)                                         gCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCt                                                                                                                         >  7:101821/1‑136 (MQ=255)                                         gCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCt                                                                                                                         <  8:101821/136‑1 (MQ=255)                                           tAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTgcg                                                                                                          >  2:333536/1‑149 (MQ=255)                                                     gTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCa                                                                                                                               >  2:50960/1‑118 (MQ=255)                                                     gTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCa                                                                                                                               <  1:50960/118‑1 (MQ=255)                                                        cATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCa                                                                                                     <  2:142135/141‑1 (MQ=255)                                                        cATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCa                                                                                                     >  1:142135/1‑141 (MQ=255)                                                        cATCCATGCAGATTTCGGTAAATGAGGTCCTGACTCCCCGACACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCTCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgtt                                                                                              <  8:153405/148‑1 (MQ=255)                                                         aTCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACATGTGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCtctc                                                                                                                                                >  1:78877/1‑98 (MQ=255)                                                         aTCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCtctc                                                                                                                                                <  2:78877/97‑1 (MQ=255)                                                         aTCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttgt                                                                                            >  6:42261/1‑149 (MQ=255)                                                              tGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGtcctc                                                                                        >  6:178064/1‑148 (MQ=255)                                                                cAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttg                                                                                             <  8:257405/141‑1 (MQ=255)                                                                cAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCtgttg                                                                                             >  7:257405/1‑141 (MQ=255)                                                                cAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCAt                                                                                     <  2:158959/149‑1 (MQ=255)                                                                 aGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATg                                                                                    <  5:365998/149‑1 (MQ=255)                                                                    tttCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCt                                                                                 <  4:232876/149‑1 (MQ=255)                                                                        ggTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCt                                                                             >  5:35981/1‑149 (MQ=255)                                                                            aaTGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGc                                                                         <  2:288515/149‑1 (MQ=255)                                                                                    ccTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTc                                                                  >  2:378572/1‑148 (MQ=255)                                                                                     cTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTc                                                                  >  2:278407/1‑147 (MQ=255)                                                                                     cTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTc                                                                  >  6:126757/1‑147 (MQ=255)                                                                                     cTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGa                                                                <  6:35981/149‑1 (MQ=255)                                                                                       gACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAgg                                                              >  2:300956/1‑149 (MQ=255)                                                                                        aCTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGc                                                             <  1:300956/149‑1 (MQ=255)                                                                                          tCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGcc                                                            >  7:180426/1‑148 (MQ=255)                                                                                          tCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGcc                                                            >  1:278330/1‑148 (MQ=255)                                                                                            cccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAgg                                                              <  5:205978/144‑1 (MQ=255)                                                                                            cccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAgg                                                              >  6:205978/1‑144 (MQ=255)                                                                                            cccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGcc                                                            >  3:370249/1‑146 (MQ=255)                                                                                            cccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCga                                                          <  2:59201/148‑1 (MQ=255)                                                                                             ccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                         >  5:916/1‑68 (MQ=255)                                                                                             ccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTg                                                                                                                                         <  6:916/68‑1 (MQ=255)                                                                                             ccGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCgaga                                                        <  4:370249/149‑1 (MQ=255)                                                                                              agCCAGAT‑TGATGTGCAGGTTGTCAGTCCGACCGGGTCCAAAATTACGCTCGAGCCTCTCGAGCGTGGATTCGGCCTTACCCTGGGCAACGCGCTGCGCCGCATCCTGGTGTCCTGCATGCCTGGGTGTGCAGTAGTCGAGGCCGAGAt                                                       <  5:222323/148‑1 (MQ=255)                                                                                                ccACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCg                                                     <  1:378572/149‑1 (MQ=255)                                                                                                 cacaT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGa                                                    <  8:180426/149‑1 (MQ=255)                                                                                                    aT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATgg                                                 >  8:222013/1‑149 (MQ=255)                                                                                                     t‑tGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTg                                                                                                            >  6:231095/1‑89 (MQ=255)                                                                                                     t‑tGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTg                                                                                                            <  5:231095/89‑1 (MQ=255)                                                                                                            tgCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACt                                           <  4:243046/148‑1 (MQ=255)                                                                                                               aGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCAt                                       <  5:42261/149‑1 (MQ=255)                                                                                                               aGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCAt                                       <  5:86700/149‑1 (MQ=255)                                                                                                                      cagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCTGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGa                                                                >  7:336499/1‑117 (MQ=255)                                                                                                                      cagtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCTGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGa                                                                <  8:336499/117‑1 (MQ=255)                                                                                                                        gtcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTcc                              >  3:219669/1‑149 (MQ=255)                                                                                                                         tcCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCAACCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCcg                             >  1:121165/1‑149 (MQ=255)                                                                                                                              cccGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATc                        <  1:333536/149‑1 (MQ=255)                                                                                                                               ccGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATc                        >  5:19778/1‑148 (MQ=255)                                                                                                                                   gcCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAgg                   >  5:22129/1‑149 (MQ=255)                                                                                                                                       aaaTTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGtt               <  6:192652/149‑1 (MQ=255)                                                                                                                                           caCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAg            >  7:171114/2‑148 (MQ=255)                                                                                                                                            aCGCTCGAGAGGCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAgg           >  2:175829/1‑148 (MQ=255)                                                                                                                                                tCGAGCCTCTCGAGCGTGGTTTCGGCCAGACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGAc      >  3:216491/1‑149 (MQ=255)                                                                                                                                                  gAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGAt                         >  4:233710/1‑128 (MQ=255)                                                                                                                                                  gAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGAt                         <  3:233710/128‑1 (MQ=255)                                                                                                                                                    gCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTaa  <  2:121165/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCAACGAAAACCTGATTGTCGAGCTCTACTCCAAGTAAGGGCTAGAAAATAGGTGCATCCATGCAGATTTCGGTAAATGAGTTCCTGACTCCCCGCCACAT‑TGATGTGCAGGTAGTCAGTCCGACCCGCGCCAAAATTACGCTCGAGCCTCTCGAGCGTGGTTTCGGCCATACCCTGGGCAACGCGCTGCGCCGCATCCTGTTGTCCTCCATGCCTGGCTGTGCAGTAGTCGAGGCCGAGATCGATGGCGTACTCCATGAGTACTCCGCGATCGAAGGTGTTCAGGAAGACGTAA  >  NC_002947/562858‑563152

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