Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,574,303	G→T	6.9%	E353D (GAG→GAT)	PP_2257 →	aerotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,574,303	0	G	T	6.9%	186.5 / 3.2	58	E353D (GAG→GAT)	PP_2257	aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base T (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CGAGGCACACACGGCGCACTTGCTGGCCGAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCGGGCACTGGCCGACAAGACCCGCCA  >  NC_002947/2574156‑2574449                                                                                                                                                     |                                                                                                                                                   cgAGGCACACACGGCGCACTTGCTGGCCGAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAgg                                                                                                                                                   <  7:136038/149‑1 (MQ=255)      cacacaCGGCGCACTTGCTGGCCGAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGgccgcc                                                                                                                                              >  5:76189/1‑149 (MQ=255)       acacacGGCGCACTTGCTGGCCGAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGTGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGccgccg                                                                                                                                             >  1:127882/1‑149 (MQ=255)       acacacGGCGCACTTGCTGGCCGAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGGCCTCGGCGATGgcgccg                                                                                                                                             <  1:235632/149‑1 (MQ=255)       acacacGGCGCACTTGCTGGCCGAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGGCCTCGGCGATGgcgccg                                                                                                                                             >  2:235632/1‑149 (MQ=255)        cacacGGCGCACTTGCTGGCCGAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCgg                                                                                                                                            <  4:112174/149‑1 (MQ=255)           acGGCGCACTTGCTGGCCGAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGga                                                                                                                                          <  1:305325/148‑1 (MQ=255)                   cTTGCTGGCCGAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCcgcg                                                                                                                                 <  1:190522/149‑1 (MQ=255)                             gAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGcc                                                                                                                         >  8:330810/1‑147 (MQ=255)                             gAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGa                                                                                                                       <  5:105534/149‑1 (MQ=255)                             gAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGa                                                                                                                       <  2:249741/149‑1 (MQ=255)                               gCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAGCACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAg                                                                                                                      <  1:264179/149‑1 (MQ=255)                               gCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCTGAGGCCGCCGGGATGATCCGCGGCATTGCCGAg                                                                                                                      >  1:178306/1‑148 (MQ=255)                               gCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGc                                                                                                                     <  7:293133/149‑1 (MQ=255)                                        ccAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCCCCTGGCGGGTCAG‑ACCGGCGCCAGCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGCGCAaaccaccc                                                                                                            >  4:325040/1‑149 (MQ=255)                                        ccAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGgatg                                                                                                                                        >  7:337039/1‑121 (MQ=255)                                        ccAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGgatg                                                                                                                                        <  8:337039/121‑1 (MQ=255)                                          aGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGCGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACctg                                                                                                          >  7:89344/1‑149 (MQ=255)                                                  gATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCt                                                                                                  <  6:64707/149‑1 (MQ=255)                                                      ccAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGcc                                                                                                                         <  4:229251/122‑1 (MQ=255)                                                      ccAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGcc                                                                                                                         >  3:229251/1‑122 (MQ=255)                                                       cAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAAcg                                                                                             >  7:194263/1‑149 (MQ=255)                                                        aGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAAcg                                                                                             <  5:140798/148‑1 (MQ=255)                                                        aGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAAcg                                                                                             >  6:140798/1‑148 (MQ=255)                                                                ggCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCg                                                                                    <  7:330810/149‑1 (MQ=255)                                                                 gCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATc                                                                                     >  5:254335/1‑147 (MQ=255)                                                                  ccATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAg                                                                                  <  6:76189/149‑1 (MQ=255)                                                                        gTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCACCGCCGCCAGCGAGGCTGcc                                                                            >  8:262063/1‑149 (MQ=255)                                                                           aCCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGt                                                                         >  3:243625/1‑149 (MQ=255)                                                                                gCCAGCACCGTGACGCAGCTCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCgg                                                                     <  3:325040/148‑1 (MQ=255)                                                                                 ccAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACgccgcc                                                                                        <  2:220055/128‑1 (MQ=255)                                                                                 ccAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACgccgcc                                                                                        >  1:220055/1‑128 (MQ=255)                                                                                    gaaCCGGGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑CCCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCACCCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAg                                                                <  8:89344/147‑1 (MQ=255)                                                                                               cAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCgggg                                                     >  4:206388/1‑149 (MQ=255)                                                                                                aGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCTGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGc                                                    <  2:178306/149‑1 (MQ=255)                                                                                                 gATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCt                                                   >  2:76847/1‑149 (MQ=255)                                                                                                       ccAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAACCCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCg                                             >  7:230883/1‑149 (MQ=255)                                                                                                                  gCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCcgac                                  >  1:191720/1‑149 (MQ=255)                                                                                                                   ccACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCggg                                                                    <  1:160482/114‑1 (MQ=255)                                                                                                                   ccACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCggg                                                                    >  2:160482/1‑114 (MQ=255)                                                                                                                   ccACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGcc                                              <  8:86194/136‑1 (MQ=255)                                                                                                                   ccACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGcc                                              >  7:86194/1‑136 (MQ=255)                                                                                                                     aCCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCgggg                                                     <  7:260418/127‑1 (MQ=255)                                                                                                                     aCCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCgggg                                                     >  8:260418/1‑127 (MQ=255)                                                                                                                        tGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTg                            >  7:33501/1‑149 (MQ=255)                                                                                                                         ggCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCACACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGc                           >  7:213597/1‑149 (MQ=255)                                                                                                                          gCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGc                           <  2:127882/148‑1 (MQ=255)                                                                                                                            gggTCCG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGTCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCggg                        <  7:228987/149‑1 (MQ=255)                                                                                                                              gTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCGGGCa                      >  2:233586/1‑149 (MQ=255)                                                                                                                               tCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCGGGCAc                     >  4:313375/1‑149 (MQ=255)                                                                                                                                 aG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCg                          <  5:313200/142‑1 (MQ=255)                                                                                                                                 aG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCg                          >  6:313200/1‑142 (MQ=255)                                                                                                                                 aG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCGGGCACTg                   <  2:191720/149‑1 (MQ=255)                                                                                                                                   tcccGGCGGCCTCGGCGATGGCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCGGGCACTgg                  >  7:279939/3‑149 (MQ=25)                                                                                                                                     ccGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCGGGCACTGGcc                <  4:200859/149‑1 (MQ=255)                                                                                                                                          ggcctCGGCGATGGCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCGGGCACTGGCCGACaa           <  8:279939/145‑1 (MQ=255)                                                                                                                                               cGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGATCGCCGACGAGGTGCGGGCACTGGCCGACAAGAccc      >  5:367948/1‑149 (MQ=255)                                                                                                                                                 ccGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCGGGCACTGGCCGACAAGACCCGc    <  6:254335/149‑1 (MQ=255)                                                                                                                                                   gAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCGGGCACTGGCCGACAAGACCCGCCa  <  8:194263/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CGAGGCACACACGGCGCACTTGCTGGCCGAGCAAGGTAGCCAGGTGGCCGATGCCAGCGGTGCGGCCATTCGTAACCTGGCCAGCACCGTGACGCAGATCAACCAGGCGGTCAGCCACCTGGCGGGTCAG‑ACCGGCGCCATCGCCGAGGCCGCCGGGATGATCCGCGGCATTGCCGAGCAAACCAACCTGCTGGCCCTCAACGCCGCCATCGAGGCTGCCCGTGCGGGCGAGCAGGGCCGGGGCTTTGCCGTGGTCGCCGACGAGGTGCGGGCACTGGCCGACAAGACCCGCCA  >  NC_002947/2574156‑2574449

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