Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,623,488	G→C	6.8%	A27P (GCG→CCG)	PP_2295 →	antirestriction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,623,488	0	G	C	6.8%	189.5 / 2.9	59	A27P (GCG→CCG)	PP_2295	antirestriction protein
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base C (2/2);  total (29/30)
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.08e-01

CGTTACCACATACGATATATTCACACGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATAAGACCGGCCAGGTGG  >  NC_002947/2623340‑2623635                                                                                                                                                      |                                                                                                                                                     cGTTACCACATACGATATATTCACACGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑acg                                                                                                                                                      <  4:118970/149‑1 (MQ=255)         caTACGATATATTCACACGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑acgcggacgcg                                                                                                                                              >  3:366628/1‑149 (MQ=255)          aTACGATATATTCACACGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑acgcggacgcg                                                                                                                                              >  3:163201/1‑148 (MQ=255)            aCGATATATTCACACGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACgcggacgcggac                                                                                                                                           >  5:97983/1‑149 (MQ=255)              gATATATTCACACGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACgcggacgcggac                                                                                                                                           >  7:202430/1‑147 (MQ=255)              gATATATTCACACGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGa                                                                                                                                         <  8:244060/149‑1 (MQ=255)                    ttCACACGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGa                                                                                                                                    >  7:14722/1‑148 (MQ=255)                        cacGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAg                                                                                                                               <  6:97983/149‑1 (MQ=255)                        cacGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAg                                                                                                                               <  4:203295/149‑1 (MQ=255)                        cacGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAg                                                                                                                               <  2:123391/149‑1 (MQ=255)                                gAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACTCACGCGGACGCGGACGAGCTGAAACAGGCCATTg                                                                                                                        <  1:49472/149‑1 (MQ=255)                                gAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGc                                                                                                                       >  2:43255/1‑149 (MQ=255)                                      aCTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCt                                                                                                                  <  1:43255/148‑1 (MQ=255)                                        tACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACg                                                                                                               >  8:32708/1‑149 (MQ=255)                                          cGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑t                                                                                                             <  8:202430/148‑1 (MQ=255)                                                cacGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTgcg                                                                                                       >  5:365130/1‑148 (MQ=255)                                                   gCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACCCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCAcc                                                                                                   >  2:297735/1‑149 (MQ=255)                                                   gCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACATGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACCCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCAc                                                                                                    >  6:251736/1‑148 (MQ=255)                                                       tACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGcc                                                                                               >  2:321347/1‑149 (MQ=255)                                                         ccTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGccc                                                                                              <  8:14722/148‑1 (MQ=255)                                                         ccTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGccc                                                                                              <  7:32708/148‑1 (MQ=255)                                                          cTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCa                                                                                             >  5:178382/1‑148 (MQ=255)                                                                   cacaATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACGCTCCTGCGCACCAGCCCATACCCTAAc                                                                                    <  1:4672/149‑1 (MQ=255)                                                                    acaATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGt                                                                                  >  5:328181/1‑149 (MQ=255)                                                                     caATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTc                                                                                                                   >  1:23915/1‑116 (MQ=255)                                                                     caATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTc                                                                                                                   <  2:23915/116‑1 (MQ=255)                                                                      aaTGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCa                                                                                <  4:163201/149‑1 (MQ=255)                                                                         gAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTAc                                                                                                                <  4:285094/115‑1 (MQ=255)                                                                               gTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCt                                                                       >  7:88178/1‑149 (MQ=255)                                                                                tACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTg                                                                      <  6:365130/149‑1 (MQ=255)                                                                                        ctgcctgGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTAACTGCCCAGCGt                                                              >  3:341481/1‑149 (MQ=255)                                                                                          gcctgGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTg                                                             >  7:348788/1‑148 (MQ=255)                                                                                           cctgGCCTCCTACAACAACGGCCACCTGCACGGCCGTTTGTTCGACCTGGAGGACT‑ACCCGGACGCGGACGAGCTGAATCAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCa                                                           <  5:251736/149‑1 (MQ=255)                                                                                           cctgGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACCCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCa                                                           <  1:297735/149‑1 (MQ=255)                                                                                                  tcctACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGt                                                      >  3:315127/1‑147 (MQ=255)                                                                                                  tcctACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGt                                                      >  6:308829/1‑147 (MQ=255)                                                                                                    ctACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTaag                                                   >  4:111794/1‑148 (MQ=255)                                                                                                    ctACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTaa                                                    <  5:308829/147‑1 (MQ=255)                                                                                                    ctACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTaa                                                    <  1:129163/147‑1 (MQ=255)                                                                                                       caacaaCGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTaagaa                                                 <  1:321347/147‑1 (MQ=255)                                                                                                       caacaaCGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGc                                               >  2:258657/1‑149 (MQ=255)                                                                                                       caacaaCGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGc                                               >  4:161200/1‑149 (MQ=255)                                                                                                        aacaacGGCCACCTGCACGGCCGTTTGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGcc                                              <  7:251849/149‑1 (MQ=255)                                                                                                         acaacGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGccc                                             >  3:284655/1‑149 (MQ=255)                                                                                                             cGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑ACCTGCGCACCACCCCATACCCTAACGTCACCCTTACCTGCCCAGCGAGCAGCCGTAAGAAGCCCGAg                                          >  8:272143/1‑148 (MQ=255)                                                                                                                ccACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGt                                      <  4:315127/149‑1 (MQ=255)                                                                                                                 cACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGCCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGt                                      >  4:62431/1‑148 (MQ=255)                                                                                                                     tGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACACTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGttcttc                                 >  3:125760/1‑149 (MQ=255)                                                                                                                             cGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGtg                         >  8:61662/1‑149 (MQ=255)                                                                                                                               ttGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGc                              >  6:375391/1‑142 (MQ=255)                                                                                                                               ttGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGc                              <  5:375391/142‑1 (MQ=255)                                                                                                                                 ggTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGtgctgc                     >  4:231922/1‑149 (MQ=255)                                                                                                                                  gTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGtgctgc                     <  3:161200/148‑1 (MQ=255)                                                                                                                                   ttCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTAc                                                                         >  3:121011/1‑95 (MQ=255)                                                                                                                                     cGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGAt                  <  5:354410/148‑1 (MQ=255)                                                                                                                                      gACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATaa                <  7:61662/149‑1 (MQ=255)                                                                                                                                       aCCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGtgctgc                     >  8:4248/1‑143 (MQ=255)                                                                                                                                       aCCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGtgctgc                     <  7:4248/143‑1 (MQ=255)                                                                                                                                       aCCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATAAg               >  7:5925/1‑149 (MQ=255)                                                                                                                                            gaggaCT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATAAGAcc            <  6:328181/147‑1 (MQ=255)                                                                                                                                            gaggaCT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATAAGACCgg          <  8:88178/149‑1 (MQ=255)                                                                                                                                            gaggaCT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATAAGACCg           <  8:5925/148‑1 (MQ=255)                                                                                                                                                  t‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATAAGACCg           <  3:258120/142‑1 (MQ=255)                                                                                                                                                  t‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATAAGACCg           >  4:258120/1‑142 (MQ=255)                                                                                                                                                     cgcggacgcggaCGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATAAGACCGGCCAGGt    <  4:366628/147‑1 (MQ=255)                                                                                                                                                     cgcggacgcggaCGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATAAGACCGGCCAGGTgg  >  2:282445/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CGTTACCACATACGATATATTCACACGGTCGGAAGTGACTACGAACACACGCACTACCTGGAGCGACACAATGAAGATGTACGCAGCCTGCCTGGCCTCCTACAACAACGGCCACCTGCACGGCCGTTGGTTCGACCTGGAGGACT‑ACGCGGACGCGGACGAGCTGAAACAGGCCATTGCGGTCTACG‑TCCTGCGCACCAGCCCATACCCTAACGTCACCGTTACCTGCCCAGCGTGCAGCGGTAAGAAGCCCGAGCAGTTCTTCTGCAGGTGCTGCGATAAGACCGGCCAGGTGG  >  NC_002947/2623340‑2623635

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