Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,135,888	G→C	10.9%	L228V (CTC→GTC)	PP_4521 ←	aerotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,135,888	0	G	C	10.9%	231.2 / 17.7	83	L228V (CTC→GTC)	PP_4521	aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (34/40);  new base C (4/5);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

TCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCATCAG  >  NC_002947/5135740‑5136035                                                                                                                                                     |                                                                                                                                                    tcGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAg                                                                                                                                                     >  6:358023/1‑149 (MQ=255)     tgcgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTgg                                                                                                                                                 >  8:168227/1‑149 (MQ=255)       cgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGGCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTgg                                                                                                                                                 <  4:357393/147‑1 (MQ=255)       cgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcc                                                                                                                                               <  3:349978/149‑1 (MQ=255)        gctgcCGGTCGAGTATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGccc                                                                                                                                              <  5:299024/149‑1 (MQ=255)        gctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcca                                                                                                                                              <  8:25363/149‑2 (MQ=255)        gctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcca                                                                                                                                              <  4:318813/149‑2 (MQ=255)         ctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCt                                                                                                                                             >  2:158742/1‑149 (MQ=255)         ctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCt                                                                                                                                             <  6:283907/149‑1 (MQ=255)         ctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGccat                                                                                                                                             <  1:366039/149‑3 (MQ=255)                  gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAg                                                                                                                                         >  6:47550/1‑144 (MQ=255)                  gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAg                                                                                                                                         <  5:47550/144‑1 (MQ=255)                    gCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTggg                                                                                                                                  <  1:100912/149‑1 (MQ=255)                    gCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTggg                                                                                                                                  >  4:13020/1‑149 (MQ=255)                     cATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGt                                                                                                                                 <  5:251622/149‑1 (MQ=255)                     cATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcca                                                                                                                                              >  7:339789/1‑135 (MQ=255)                     cATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcca                                                                                                                                              <  8:339789/136‑2 (MQ=255)                        cTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGt                                                                                                                              >  4:77843/1‑149 (MQ=255)                           cgcCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTc                                                                                                                             <  3:304565/147‑1 (MQ=255)                           cgcATGGGCATCGGAGAAGTCGGTAGATTGAGCCGCCCGCGACTTTACCCCGACCGCGGCATCCACCAGGGGTGTCAGGGGCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCTAGCTGGCCCTTCTGGCCTGGGTTGTc                                                                                                                             <  8:199008/147‑1 (MQ=255)                             ccTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTg                                                                                                                         >  8:87787/1‑149 (MQ=255)                             ccTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTAGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTg                                                                                                                         >  4:29758/1‑149 (MQ=255)                              cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTg                                                                                                                         >  8:69424/1‑148 (MQ=255)                                gggCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGCCCCCGGCAGCAACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTGCTTCGCGGATGATGCCCCGGGGGAGCTGGCCGTGCGGGCGTGGGTTggggatac                                                                                                                        >  4:104192/1‑140 (MQ=255)                                 ggCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAg                                                                                                                      <  3:124891/148‑1 (MQ=255)                                  gCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGt                                                                                                                    <  7:168227/149‑1 (MQ=255)                                    aTCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGt                                                                                                                    >  3:241776/1‑147 (MQ=255)                                    aTCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCa                                                                                                                  <  3:77843/149‑1 (MQ=255)                                        gcccaGTCGGCAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCCGGCGGGTCAGGGCCGTTGTCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGGTGTCGCTGTAGGTCAgcgc                                                                                                              <  2:259461/145‑1 (MQ=255)                                           aaGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGccac                                                                                                           >  5:67868/1‑149 (MQ=255)                                             gTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGc                                                                                                                           >  4:205658/1‑131 (MQ=255)                                             gTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGc                                                                                                                           <  3:205658/131‑1 (MQ=255)                                             gTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGccacc                                                                                                          >  1:278593/1‑148 (MQ=255)                                                ggAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACg                                                                                                      >  8:271290/1‑149 (MQ=255)                                                ggAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACg                                                                                                      >  1:85733/1‑149 (MQ=255)                                                   aGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGAc                                                                                                   <  3:13020/149‑1 (MQ=255)                                                      ttGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTgg                                                                                                <  6:290058/149‑1 (MQ=255)                                                         cgccgcCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTa                                                                                              <  2:48838/148‑1 (MQ=255)                                                         cgccgcCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTa                                                                                              >  4:130449/1‑148 (MQ=255)                                                          gccgccCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTaaa                                                                                            <  7:271290/149‑1 (MQ=255)                                                          gccgccCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTa                                                                                              >  8:66406/1‑147 (MQ=255)                                                           ccgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAg                                                                                           <  7:192982/149‑1 (MQ=255)                                                            cgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAgg                                                                                          >  1:160629/1‑149 (MQ=255)                                                               ccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTAAGCGCCACCAACGGCCTGGTAACGGCCt                                                                                       >  8:327736/1‑149 (MQ=255)                                                                cgcgACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCtt                                                                                      <  2:58665/149‑1 (MQ=255)                                                                   gACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTggg                                                                                   <  4:325498/149‑1 (MQ=255)                                                                       ttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGc                                                                               >  3:186118/1‑149 (MQ=255)                                                                           cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTa                                                                                                                       <  4:259558/105‑1 (MQ=255)                                                                           cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTa                                                                                                                       >  3:259558/1‑105 (MQ=255)                                                                           cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAg                                                                           >  1:333222/1‑149 (MQ=255)                                                                           cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAg                                                                           >  7:205048/1‑149 (MQ=255)                                                                             cgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGACGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGCCTGGTAAAGGCCTTGGGATGCTCAGc                                                                          >  8:342039/1‑148 (MQ=255)                                                                                ccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAAc                                                                      <  2:316170/149‑1 (MQ=255)                                                                                   cgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTg                                                                   <  3:29758/149‑1 (MQ=255)                                                                                    ggCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGCTGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGc                                                                  <  5:358023/149‑1 (MQ=255)                                                                                     gCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGcg                                                                 <  6:67868/149‑1 (MQ=255)                                                                                       aTCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGcgac                                                               <  1:158742/149‑1 (MQ=255)                                                                                          caccaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGaa                                                            >  4:22469/1‑149 (MQ=255)                                                                                           accaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGc                                                                  <  8:117593/142‑1 (MQ=255)                                                                                           accaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGc                                                                  >  7:117593/1‑142 (MQ=255)                                                                                           accaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAAt                                                           >  3:185687/1‑149 (MQ=255)                                                                                            ccaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATg                                                          <  6:165160/149‑1 (MQ=255)                                                                                                gCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCt                                                      <  7:66406/149‑1 (MQ=255)                                                                                                        gggCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGc                                              <  4:185687/149‑1 (MQ=255)                                                                                                           ccGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGGCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCgg                                           <  7:327736/149‑1 (MQ=255)                                                                                                              tttGCCGACGGGCTTCTTCGCGTATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGc                                        <  8:73813/149‑1 (MQ=255)                                                                                                              tttGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGc                                        >  5:68634/1‑149 (MQ=255)                                                                                                                 gCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTc                                      <  2:333222/148‑1 (MQ=255)                                                                                                                    gACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGAcc                                  <  7:87787/149‑1 (MQ=255)                                                                                                                      cGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCa                                <  3:130449/149‑1 (MQ=255)                                                                                                                         gCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTa                             >  6:120382/1‑149 (MQ=255)                                                                                                                          cttcttCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTAc                            >  6:254391/1‑149 (MQ=255)                                                                                                                              ttcGATGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCAGCGTTGACGCTGTAGCTCAGCGCCACCAACGGACTGGTAAAGTCCTTGGGATGCTCAGCCAACGTGc                                                                  >  7:315429/1‑107 (MQ=255)                                                                                                                                  cTGATGACGGCCAGGTCGAGCTGGCCCTTCATGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCATCTGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcag                    <  7:342039/149‑1 (MQ=255)                                                                                                                                     tgaagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcag                 >  7:204792/6‑149 (MQ=255)                                                                                                                                     tgaagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGACGCTGTAGGTCAGCGCCACCACCGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcaccag                 >  7:139386/6‑149 (MQ=255)                                                                                                                                     atgatgGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcag                 <  4:241776/149‑1 (MQ=255)                                                                                                                                     atgatgGCCAGGGCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcag                 <  2:278593/149‑1 (MQ=255)                                                                                                                                      gaagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcag                 >  1:20645/5‑148 (MQ=255)                                                                                                                                           gCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCGCGACCCAGTACCAAGCCAGCAGCAGa             <  2:85733/147‑1 (MQ=255)                                                                                                                                            ccAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACtt          >  5:221862/1‑149 (MQ=255)                                                                                                                                             cAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTAGATGTAGGTCAGCGACACCAACGGACTGGTAAAGGCCTTGGGATGCTAAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAATCCAGCAGCAGACTTg         >  3:71213/1‑149 (MQ=255)                                                                                                                                               ggTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGcc       <  5:254391/149‑1 (MQ=255)                                                                                                                                                gTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCa      >  1:105130/1‑149 (MQ=255)                                                                                                                                                     gCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCATCAg  >  3:135916/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCATCAG  >  NC_002947/5135740‑5136035

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