Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,709,381	C→T	5.9%	V1133V (GTC→GTT)	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,709,381	0	C	T	5.9%	256.2 / 3.8	85	V1133V (GTC→GTT)	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/41);  new base T (3/2);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-01

CAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGACGTATTC  >  NC_002947/5709236‑5709527                                                                                                                                                  |                                                                                                                                                    cAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCcgc                                                                                                                                                  <  5:299130/149‑1 (MQ=255) cAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCcg                                                                                                                                                   <  1:351104/148‑1 (MQ=255) cAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCcg                                                                                                                                                   <  7:88076/148‑1 (MQ=255)  aGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCcgcg                                                                                                                                                 <  8:337128/149‑1 (MQ=255)  aGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCcgc                                                                                                                                                  >  8:168071/1‑148 (MQ=255)    ccATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGAt                                                                                                                                               >  2:265459/1‑149 (MQ=255)      aTCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGGCGACCACCACGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTc                                                                                                                                             <  7:111223/149‑1 (MQ=255)          aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                         <  1:187614/149‑1 (MQ=255)          aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                         >  5:174708/1‑149 (MQ=255)          aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                         >  8:263922/1‑149 (MQ=255)          aaaaCGACATGCTGCTGGGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGGACCCGTTCGCGAtaatcga                                                                                                                                         <  2:269833/149‑7 (MQ=255)            aaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATTCAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagcca                                                                                                                                       <  1:265459/149‑1 (MQ=255)            aaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                         >  6:129664/1‑147 (MQ=255)             ccGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagccag                                                                                                                                      <  4:43286/148‑1 (MQ=255)             aCGAAATGCTGCTGGGCCGCCGGTTGCTGGAGCAGGTCGACCACCAGGCGCTGCGCGAACTGTCCATCGGGCTGCCCCGCGGTCCAGCGACTGCCGAGCATATACAGCAAAACCTCTATTATCGCGAACGGGTCCGCGATTccagccag                                                                                                                                      <  2:268911/149‑1 (MQ=255)               gACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAgccaggc                                                                                                                                    >  1:285920/1‑149 (MQ=255)                aCATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGa                                                                                                                                                <  4:135830/136‑1 (MQ=255)                aCATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGa                                                                                                                                                >  3:135830/1‑136 (MQ=255)                aCATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAgccaggc                                                                                                                                    >  1:96120/1‑148 (MQ=255)                        ctgTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAgccgcc                                                                                                                            <  5:129664/148‑1 (MQ=255)                           tgCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAgccgccgc                                                                                                                          <  2:258638/147‑1 (MQ=255)                             ccgccgCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGc                                                                                                                      >  1:185200/1‑149 (MQ=255)                             ccgccgCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGc                                                                                                                      <  7:263922/149‑1 (MQ=255)                                       cttgAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCTGGTTCCAGAGACTGCCGAGAAGATACAGCAATACCTCTATTATCTCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCgtca                                                                                                             <  1:157116/145‑1 (MQ=255)                                       cTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCgtca                                                                                                             <  4:80498/148‑1 (MQ=255)                                       cTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCgtca                                                                                                             <  5:370411/148‑1 (MQ=255)                                        tGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccag                                                                                                                                          <  6:343365/118‑1 (MQ=255)                                        tGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccag                                                                                                                                          >  5:343365/1‑118 (MQ=255)                                            gCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGc                                                                                                       >  8:150787/1‑149 (MQ=255)                                                gTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCt                                                                                                   <  4:9429/149‑1 (MQ=255)                                                   gACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTTCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTgga                                                                                                <  7:168071/149‑1 (MQ=255)                                                   gACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTgg                                                                                                 <  6:174708/148‑1 (MQ=255)                                                    accaccAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTggat                                                                                               <  2:185200/149‑1 (MQ=255)                                                     ccaccaGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTggatc                                                                                              <  8:342355/149‑1 (MQ=255)                                                      caccaGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCCCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTAGTCGTCAGTCAGCGGCTggatcg                                                                                             >  7:37354/1‑149 (MQ=255)                                                          aGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGTTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTggatcggatc                                                                                         <  5:187655/149‑1 (MQ=255)                                                          aGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCATCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTggatcggatc                                                                                         <  2:47844/149‑1 (MQ=255)                                                           tgcgcTGATCGCATAGTCCTTCGGGATGCCCTGCGTTCCTTAGACTGCTGAGCAGATACAGCAATACCTCGATTATCGCGTAGGGGTGGTCGTTTCGAGCGAGGGCTGGGGCCGCGTGCTTGTCGTCATTCTGCGGGTGGATCGGATcc                                                                                        <  2:251162/148‑1 (MQ=255)                                                                   gcgcACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGc                                                                                    >  3:120343/1‑145 (MQ=255)                                                                   gcgcACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGc                                                                                    <  4:120343/145‑1 (MQ=255)                                                                          gTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGAGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑Tg                                                                         >  4:12707/1‑148 (MQ=255)                                                                           tCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑Tg                                                                         >  7:31634/1‑147 (MQ=255)                                                                             cATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAg                                                                      >  8:39166/1‑148 (MQ=255)                                                                             cATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAg                                                                      >  8:45025/1‑148 (MQ=255)                                                                              aTCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGa                                                                    >  2:292407/1‑149 (MQ=255)                                                                                 gggCTGCCCCGCGTGCCAGCGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGAcc                                                                  <  5:112700/148‑1 (MQ=255)                                                                                  ggCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGAccc                                                                 >  2:315028/1‑148 (MQ=255)                                                                                     tGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCa                                                                       <  1:128759/139‑1 (MQ=255)                                                                                     tGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCa                                                                       >  2:128759/1‑139 (MQ=255)                                                                                     tGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGAcccc                                                                >  4:95769/1‑146 (MQ=255)                                                                                     tGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGAcccc                                                                >  5:259024/1‑146 (MQ=255)                                                                                          cgcgTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGa                                                        >  8:338814/1‑149 (MQ=255)                                                                                               tCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTc                                                   >  4:136998/1‑149 (MQ=255)                                                                                               tCCAGAGACTGCCCAGCAGATACAGCAATACCGCGATTATCGCGAACGGGGCCGCGATTCCCGCCAGGCCAGCCGCCGCCTGCTTGTCGGCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTc                                                   <  6:315129/149‑1 (MQ=255)                                                                                                 cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTccc                                                 <  2:285920/149‑1 (MQ=255)                                                                                                  agagACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTccc                                                 >  8:317883/1‑148 (MQ=255)                                                                                                    agaCTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAAt                                              <  1:315028/149‑1 (MQ=255)                                                                                                        tGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTc                                                   <  5:262251/140‑1 (MQ=255)                                                                                                        tGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTc                                                   >  6:262251/1‑140 (MQ=255)                                                                                                          ccGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTg                                                                                                                   >  5:69573/1‑75 (MQ=255)                                                                                                          ccGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTg                                                                                                                   <  6:69573/75‑1 (MQ=255)                                                                                                             aGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACg                                     <  7:317883/149‑1 (MQ=255)                                                                                                                aGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGaaa                                  <  2:96120/149‑1 (MQ=255)                                                                                                                 gATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAAc                                 >  2:328004/1‑149 (MQ=255)                                                                                                                    aCAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTggatcggatc                                                                                         <  1:85848/91‑1 (MQ=255)                                                                                                                    aCAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTggatcggatc                                                                                         >  2:85848/1‑91 (MQ=255)                                                                                                                     cAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGc                             <  7:338814/149‑1 (MQ=255)                                                                                                                      aGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCACCCCCCGCCAGCTTGTCGGCAGGCAGCGGCTCGAACGGATCCGGGCCG‑AACACAAGCAGGACCCCCAAAACCATGTACCCCAAAAAGCCACCaaccgacgc                            >  8:237091/1‑142 (MQ=255)                                                                                                                        cAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGcc                            >  8:299658/1‑147 (MQ=255)                                                                                                                        cAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGcc                            >  5:152630/1‑147 (MQ=255)                                                                                                                        cAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGcc                            >  6:242744/1‑147 (MQ=255)                                                                                                                              cTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCACAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCa                     >  2:171286/1‑148 (MQ=255)                                                                                                                              cTCGATTATAGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAt                    <  7:150787/149‑1 (MQ=255)                                                                                                                               tCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGa                                                        <  3:230588/112‑1 (MQ=255)                                                                                                                               tCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGa                                                        >  4:230588/1‑112 (MQ=255)                                                                                                                               tCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                          >  3:15693/1‑142 (MQ=255)                                                                                                                               tCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                          <  4:15693/142‑1 (MQ=255)                                                                                                                                 tggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTg                  >  3:140962/5‑148 (MQ=255)                                                                                                                                 tggaaTCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTg                  >  3:146504/5‑148 (MQ=255)                                                                                                                                 tggaaTCGCGAACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCa                                                                       >  7:134787/5‑95 (MQ=255)                                                                                                                                 tggaaTCGCGAACCCGTTCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCa                                                                       <  8:134787/91‑1 (MQ=255)                                                                                                                                 gATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt                   <  6:259024/147‑1 (MQ=255)                                                                                                                                   ttATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAATATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGaa                <  3:12707/148‑1 (MQ=255)                                                                                                                                                gTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGACGTATTc  >  4:185622/1‑149 (MQ=255)                                                                                                                                                gTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGACGTATTc  >  5:230686/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACAC‑TGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGACGTATTC  >  NC_002947/5709236‑5709527

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