Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,284,431	G→A	6.4%	R279H (CGT→CAT)	PP_4658 →	methyl‑accepting chemotaxis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,284,431	0	G	A	6.4%	187.1 / 2.1	62	R279H (CGT→CAT)	PP_4658	methyl‑accepting chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base A (2/2);  total (33/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 = 8.35e-01

GAAAACGCAGTACTGGCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTG  >  NC_002947/5284283‑5284580                                                                                                                                                      |                                                                                                                                                      gAAAACGCAGTACTGGCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCg                                                                                                                                                       >  8:255745/1‑149 (MQ=255)  aaaaCGCAGTACTGGCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGt                                                                                                                                                      <  3:115201/149‑1 (MQ=255)      cGCAGTACTGGCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACa                                                                                                                                                  >  6:247440/1‑149 (MQ=255)      cGCAGTACTGGCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACa                                                                                                                                                  <  7:255745/149‑1 (MQ=255)             cTGGCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGCCCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGagag                                                                                                                                           <  8:316218/149‑1 (MQ=255)             cTGGCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGagag                                                                                                                                           <  5:247440/149‑1 (MQ=255)              tGGCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGagag                                                                                                                                           >  3:276174/1‑148 (MQ=255)                gCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCa                                                                                                                                         <  7:108455/148‑1 (MQ=255)                gCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCAc                                                                                                                                        >  2:133066/1‑149 (MQ=255)                gCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTGGTACCAGATCCAGCTGGAGCGCCATGACATCga                                                                                                                                              <  6:17262/143‑1 (MQ=255)                gCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTGGTACCAGATCCAGCTGGAGCGCCATGACATCga                                                                                                                                              >  5:17262/1‑143 (MQ=255)                   tATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCCCCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACTTCGAGAGCAcc                                                                                                                                       >  5:282369/1‑147 (MQ=255)                   tATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCAcc                                                                                                                                       >  1:188989/1‑147 (MQ=255)                    aTCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAg                                                                                                                                    >  3:131601/1‑149 (MQ=255)                     tCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAgg                                                                                                                                   <  2:188989/149‑1 (MQ=255)                        tGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGccc                                                                                                                                <  1:133066/149‑1 (MQ=255)                         gACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCATGACATCGAGAGCACCCAGGCCCg                                                                                                                               >  6:337741/1‑149 (MQ=255)                             ccGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGc                                                                                                                           >  8:164926/1‑149 (MQ=255)                                   gCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAg                                                                                                                     <  8:31937/149‑1 (MQ=255)                                     cAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGc                                                                                                                   >  6:66903/1‑149 (MQ=255)                                      aGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGc                                                                                                                   >  5:237239/1‑148 (MQ=255)                                         tCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCAt                                                                                                                <  7:164926/148‑1 (MQ=255)                                           cgcgACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCTCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGc                                                                                                             <  6:141538/149‑1 (MQ=255)                                             cgACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCca                                                                                                           <  7:317678/149‑1 (MQ=255)                                                        ccAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGgcg                                                                                                >  1:219153/1‑149 (MQ=255)                                                         cAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGgcgc                                                                                               >  7:274800/1‑149 (MQ=255)                                                          aaCATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGgcg                                                                                                >  7:26974/1‑147 (MQ=255)                                                          aaCATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGgcg                                                                                                <  8:26974/147‑1 (MQ=255)                                                          aaCATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCgct                                                                                              >  6:326904/1‑149 (MQ=255)                                                            cATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCgctgc                                                                                            <  4:131601/149‑1 (MQ=255)                                                              tCACCACCTCGCGCGCCGAAATGACCGTGGCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCgctgct                                                                                           >  1:21929/1‑149 (MQ=255)                                                              tCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGctgctg                                                                                          <  2:21929/149‑1 (MQ=255)                                                                  caccTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGctgctg                                                                                          >  1:97424/1‑145 (MQ=255)                                                                  caccTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGctgctg                                                                                          <  2:97424/145‑1 (MQ=255)                                                                  caccTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCg                                                                                      >  6:99261/1‑149 (MQ=255)                                                                            gcCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTggc                                                                            >  6:11867/1‑149 (MQ=255)                                                                                   tGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTg                                                                     <  4:276174/149‑1 (MQ=255)                                                                                       cGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGa                                                                 <  2:219153/149‑1 (MQ=255)                                                                                                    cGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGc                                                         >  5:375328/1‑145 (MQ=255)                                                                                                        aTCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCAc                                                 >  5:190610/1‑149 (MQ=255)                                                                                                           gTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCAcgcg                                              >  5:250463/1‑149 (MQ=255)                                                                                                           gTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCAcgcg                                              >  3:92424/1‑149 (MQ=255)                                                                                                              aaGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGccc                                           >  6:214557/1‑149 (MQ=255)                                                                                                                 cgcAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGtt                                        <  4:251902/149‑1 (MQ=255)                                                                                                                 cgcAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGtt                                        <  6:250463/149‑1 (MQ=255)                                                                                                                     gTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCaa                                    <  5:99261/149‑1 (MQ=255)                                                                                                                      tGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAg                                   <  6:375328/149‑1 (MQ=255)                                                                                                                      tGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCT‑CTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGa                                  >  8:80289/1‑149 (MQ=255)                                                                                                                           cgTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACAcc                              <  5:66903/149‑1 (MQ=255)                                                                                                                                tACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCtggt                         >  1:100187/1‑149 (MQ=255)                                                                                                                                 aCCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGc                                                         >  2:14651/1‑116 (MQ=255)                                                                                                                                  ccAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATc                                                               >  4:231135/1‑109 (MQ=255)                                                                                                                                  ccAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATc                                                               <  3:231135/109‑1 (MQ=255)                                                                                                                                   cAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGa                                                                 <  2:190110/106‑1 (MQ=255)                                                                                                                                   cAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGa                                                                 >  1:190110/1‑106 (MQ=255)                                                                                                                                      aTCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGt                   >  1:190581/1‑149 (MQ=255)                                                                                                                                       tCCAGCTGGAGCGCCATGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTg                  <  5:337741/149‑1 (MQ=255)                                                                                                                                        ccAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGt                   <  4:145953/147‑1 (MQ=255)                                                                                                                                        ccAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGt                   >  3:145953/1‑147 (MQ=255)                                                                                                                                               gAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAAt          >  6:322979/1‑149 (MQ=255)                                                                                                                                               gAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAAt          >  7:152356/1‑149 (MQ=255)                                                                                                                                                   gcCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGcc      <  5:326904/149‑1 (MQ=255)                                                                                                                                                    ccGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCa     <  5:322979/149‑1 (MQ=255)                                                                                                                                                      gTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCT‑CTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTg  <  7:80289/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GAAAACGCAGTACTGGCCTATCGTGACGCCGTACGCCAGTTCCGCGACGCGGTTGCCAACATCACCACCTCGCGCGCCGAAATGACCGT‑GCAGGGCGCCGACATCGTCAAGCGCAGTGACGCGTTGTACCAGATCCAGCTGGAGCGCCGTGACATCGAGAGCACCCAGGCCCGCAGCCTGCAGGCCATCGCCACACTGTTGGCGCTGCTGGTCGGCGTGCTGGCGGCGGTGCTGATCACCCGCCAGATCACGCGCCCGTTGCAAGACACCCTGGTGGCCGTGGAGAAAATCGCCAGTG  >  NC_002947/5284283‑5284580

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