Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,431,918	A→G	7.4%	A187A (GCA→GCG)	PP_3047 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,431,918	0	A	G	7.4%	185.9 / 2.8	66	A187A (GCA→GCG)	PP_3047	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base G (3/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTACGAAGGTACGCAATACGCCGCTGCGCTGCGTGCCGATGGTCGC  >  NC_002947/3431773‑3432065                                                                                                                                                     |                                                                                                                                                    tGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGcacaa                                                                                                                                                  >  3:265415/1‑148 (MQ=255) tGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGcac                                                                                                                                                    <  4:265415/147‑1 (MQ=255)   aaCGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACag                                                                                                                                                  <  2:168031/147‑1 (MQ=255)     cGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACaggagg                                                                                                                                              <  7:391994/149‑1 (MQ=255)                  gATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTggggg                                                                                                                                  <  3:290169/148‑1 (MQ=255)                  gATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTggg                                                                                                                                    <  7:316862/146‑1 (MQ=255)                  gATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGc                                                                                                                                 <  3:367695/149‑1 (MQ=255)                   aTGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTCGGGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGc                                                                                                                                 >  1:314551/1‑149 (MQ=255)                   aTGCGGACGGTAGCAACGTGCTGCTCAATATCTATGTTCGTTTC‑TGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTTCCGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGc                                                                                                                                 <  1:320808/149‑1 (MQ=255)                      cGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGAc                                                                                                                              >  4:348686/1‑148 (MQ=255)                       ggACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCGCAGGAGGAGGCGCTGGGGGCGACTa                                                                                                                            >  3:215150/1‑149 (MQ=255)                       ggACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCGCAGGAGGAGGCGCTGGGGGCGACTa                                                                                                                            >  3:208228/1‑149 (MQ=255)                         aCGGTAGCGCCGTGCTGCTGAATGTGTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGAGGGAGCT‑GGGTGCCGCAGGTACCGAGGTT‑CGGGTGGAGTGTGGTGCCGCGCTGGGTGCGGAGGAGGAGGCGCTGGGGGCGACTAcc                                                                                                                          >  4:77918/1‑149 (MQ=255)                         aCGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTAcc                                                                                                                          >  4:39041/1‑149 (MQ=255)                               gCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAAc                                                                                                                    <  6:79817/149‑1 (MQ=255)                                cACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACg                                                                                                                   <  3:348686/149‑1 (MQ=255)                                    gTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGc                                                                                                               <  4:157489/149‑1 (MQ=255)                                        tactCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGccc                                                                                                             >  6:138041/3‑147 (MQ=255)                                          ctCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGTCGCTGGTGGCGACTACCAGCAACGGTGCCCGAGc                                                                                                         <  7:207640/149‑1 (MQ=255)                                          ctCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGc                                                                                                         <  4:228426/149‑1 (MQ=255)                                           tCAATATCTATGATCGTTTC‑GGCCTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGCAGTGTGTTGCCGCACTGGCTGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCg                                                                                                        >  7:90419/1‑149 (MQ=255)                                            cAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGt                                                                                                       >  1:160104/1‑149 (MQ=255)                                                 tCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGTGCGACTACCAGCAACGGTGCCCGAGCGTTCTg                                                                                                   <  1:271328/148‑1 (MQ=255)                                                  cTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCTAGGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGc                                                                                                  >  1:308711/1‑149 (MQ=255)                                                   tATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCgg                                                                                                >  4:249265/1‑149 (MQ=255)                                                       aTCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGt                                                                                               >  1:263945/1‑146 (MQ=255)                                                       aTCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTaa                                                                                             <  5:47177/148‑1 (MQ=255)                                                             tC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACtt                                                                                       <  3:77918/148‑1 (MQ=255)                                                                ggACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCt                                                                                     <  2:246555/148‑1 (MQ=255)                                                                 gACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTgg                                                                                   <  2:361682/149‑1 (MQ=255)                                                                   cTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTggcg                                                                                 >  1:103248/1‑149 (MQ=255)                                                                    ttCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTggcgg                                                                                >  2:272317/1‑149 (MQ=255)                                                                      cAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGaa                                                                                          <  2:221885/137‑1 (MQ=255)                                                                      cAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGaa                                                                                          >  1:221885/1‑137 (MQ=255)                                                                        gCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGgcgc                                                                            <  2:103248/149‑1 (MQ=255)                                                                         ccAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGAGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGgcgc                                                                            >  7:341456/1‑148 (MQ=255)                                                                               gagatcgagatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGAtctc                                                                     >  6:401558/1‑149 (MQ=255)                                                                               gagatcgagatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGATGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGAGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGAtctc                                                                     >  6:284975/1‑149 (MQ=255)                                                                                agatcgagatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCg                                                                    <  3:110871/149‑1 (MQ=255)                                                                                agatcgagatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCg                                                                    >  3:288499/1‑149 (MQ=255)                                                                                  atcgagatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCGCAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCa                                                                  <  4:215150/149‑1 (MQ=255)                                                                                  atcgagatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCGCAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCa                                                                  <  4:208228/149‑1 (MQ=255)                                                                                    cgagatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCt                                                                           >  4:280087/1‑138 (MQ=255)                                                                                    cgagatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCt                                                                           <  3:280087/138‑1 (MQ=255)                                                                                     gagatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATc                                                                <  8:19450/148‑1 (MQ=255)                                                                                      agatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCACTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCAt                                                                 >  2:100827/1‑146 (MQ=255)                                                                                       gatGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGAAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCagc                                                             >  3:231416/1‑148 (MQ=255)                                                                                              cT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTaa                                                      <  2:436364/149‑1 (MQ=255)                                                                                                  ggTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTaagaa                                                   <  3:249265/149‑1 (MQ=255)                                                                                                     aCCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGAc                                                >  6:3525/1‑149 (MQ=255)                                                                                                      ccGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGAcc                                               >  7:150443/1‑149 (MQ=255)                                                                                                      ccGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGAc                                                <  5:3525/148‑1 (MQ=255)                                                                                                         cAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCaga                                                             >  8:261795/1‑132 (MQ=255)                                                                                                         cAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCaga                                                             <  7:261795/132‑1 (MQ=255)                                                                                                          aGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTAc                                            >  3:184261/1‑148 (MQ=255)                                                                                                                 gATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGc                                                                                                                                 <  4:312884/56‑1 (MQ=255)                                                                                                                 gATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGc                                                                                                                                 >  3:312884/1‑56 (MQ=255)                                                                                                                 gATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTACGAAGGTAc                                    <  4:272653/149‑1 (MQ=255)                                                                                                                         gggTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTACGAAGGTACGCAATAc                             <  4:288499/149‑1 (MQ=255)                                                                                                                                gtgtgtTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGtt                                                        <  7:319795/115‑1 (MQ=255)                                                                                                                                gtgtgtTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGtt                                                        >  8:319795/1‑115 (MQ=255)                                                                                                                                  gtgtTGCCGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTACGAAGGTACGCAATACGCcgctgcg                    >  5:342829/1‑149 (MQ=255)                                                                                                                                   tgtTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTACGAAGGTACGCAATACGCcgctgcgc                   >  6:341826/1‑149 (MQ=255)                                                                                                                                       gACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGAtctc                                                                     >  7:301907/1‑95 (MQ=255)                                                                                                                                       gACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGAtctc                                                                     <  8:301907/95‑1 (MQ=255)                                                                                                                                          gCACTGGACGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTACGAAGGTACGCAATACGCCGCTGCGCTGCGTGc            <  4:293530/149‑1 (MQ=255)                                                                                                                                             cTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTACGAAGGTACGCAATACGCCGCTGCGCTGCGTGCCGa         >  4:139716/1‑149 (MQ=255)                                                                                                                                              tGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTACGAAGGTACGCAATACGCCGCTGCGCTGCGTGCCGAt        >  7:291622/1‑149 (MQ=255)                                                                                                                                                    cacaGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTACGAAGGTACGCAATACGCCGCTGCGCTGCGTGCCGATGGTcgc  <  2:160104/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGAACGGCGTTGTTCTGGATGCGGACGGTAGCACCGTGCTGCTCAATATCTATGATCGTTTC‑GGACTTCAGCCAATCGAGATCGAGATGGAGCT‑GGGTACCGCAGGTACCGATGTT‑CGGGTGAAGTGTGTTGACGCACTGGATGCACAGGAGGAGGCGCTGGGGGCGACTACCAGCAACGGTGCCCGAGCGTTCTGCGGTAAGAACTTCTGGCGGGCGCTGATCTCGCATCAGAGCGTTAAGAAGACCTACGAAGGTACGCAATACGCCGCTGCGCTGCGTGCCGATGGTCGC  >  NC_002947/3431773‑3432065

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