Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,650,577	C→T	6.4%	A113T (GCA→ACA)	PP_4113 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,650,577	0	C	T	6.4%	182.1 / 2.1	62	A113T (GCA→ACA)	PP_4113	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/25);  new base T (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GCACCGACAACCTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTGGCGCTCGATATCGGCTTCTGA  >  NC_002947/4650435‑4650724                                                                                                                                                |                                                                                                                                                     gCACCGACAACCTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGa                                                                                                                                                 >  1:2479/1‑148 (MQ=255)  cACCGACAACCTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATc                                                                                                                                               <  3:339330/149‑1 (MQ=255)   aCCGACAACCTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATc                                                                                                                                               >  3:22583/1‑148 (MQ=255)     cGACAACCTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATcac                                                                                                                                             >  6:258155/1‑148 (MQ=255)      gACAACCTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACAt                                                                                                                                           <  1:96328/149‑1 (MQ=255)      gACAACCTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACAt                                                                                                                                           >  7:271942/1‑149 (MQ=255)         aaCCTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGAc                                                                                                                                        >  4:93744/1‑149 (MQ=255)           ccTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTg                                                                                                                                      >  8:200225/1‑149 (MQ=255)           ccTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTg                                                                                                                                      >  2:220468/1‑149 (MQ=255)            cTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTgc                                                                                                                                     <  5:258155/149‑1 (MQ=255)                 cgcCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTgg                                                                                                                                >  7:122799/1‑149 (MQ=255)                    cGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCgg                                                                                                                             <  5:296563/149‑1 (MQ=255)                          tACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGc                                                                                                                       >  5:326916/1‑149 (MQ=255)                            tcGGCGAGGCCCGGCGCATCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTgg                                                                                                                                >  6:152846/2‑138 (MQ=255)                            tcGGCGAGGCCCGGCGCATCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTgg                                                                                                                                <  5:152846/137‑1 (MQ=255)                                 gATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCg                                                                                                                <  7:200225/149‑1 (MQ=255)                                  aTGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCgg                                                                                                               <  8:122799/149‑1 (MQ=255)                                   tGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCggg                                                                                                              <  3:93744/149‑1 (MQ=255)                                    gcgcCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCggg                                                                                                              <  2:377465/148‑1 (MQ=255)                                      gcCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCGTACCGGGTTg                                                                                                           <  2:241886/149‑1 (MQ=255)                                      gcCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTg                                                                                                           >  3:315973/1‑149 (MQ=255)                                           gCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGt                                                                                                      <  4:315973/149‑1 (MQ=255)                                            ccTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGTGGCGATCACATGACTGCGCTGGCGGGTCAGCCATAc                                                                                                                  <  6:26294/136‑1 (MQ=255)                                            ccTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGTGGCGATCACATGACTGCGCTGGCGGGTCAGCCATAc                                                                                                                  >  5:26294/1‑136 (MQ=255)                                                 ccGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGgcc                                                                                                <  8:271942/149‑1 (MQ=255)                                                       ccGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑c                                                                                         <  4:22583/149‑1 (MQ=255)                                                        cgtcgtATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑Cg                                                                                        >  2:122021/1‑149 (MQ=255)                                                          tcgtATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTc                                                                                                                          <  7:182277/114‑1 (MQ=255)                                                          tcgtATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTc                                                                                                                          >  8:182277/1‑114 (MQ=255)                                                          tcgtATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGgc                                                                                      >  5:306134/1‑149 (MQ=255)                                                            gtATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGgcca                                                                                    <  5:310051/149‑1 (MQ=255)                                                              aTTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATca                                                                                                                                              >  4:121223/1‑90 (MQ=255)                                                                       cGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGATGGTGCGGCCGCCAT‑CGGCCAGCCAGAGc                                                                            <  1:220468/146‑1 (MQ=255)                                                                        gATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGTGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCgggg                                                                        >  8:325410/1‑149 (MQ=255)                                                                             acaaGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGAGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCAAGCCAGAGCGGGgcgatg                                                                   >  3:192993/1‑149 (MQ=255)                                                                                   tGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTgg                                                             >  4:2803/1‑149 (MQ=255)                                                                                    gACCAGATCGACATAAGTGGTCAGGCCATCGCATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCGGCCAGAGCGGGGCGCTGCGCTgg                                                             >  1:226356/1‑149 (MQ=255)                                                                                    gACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTggg                                                            >  6:171143/1‑149 (MQ=255)                                                                                     aCCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGc                                                           >  7:106308/1‑149 (MQ=255)                                                                                      ccAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGc                                                           <  6:326916/148‑1 (MQ=255)                                                                                              aCATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGc                                                  >  2:366963/1‑149 (MQ=255)                                                                                                  aaGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGt                                              >  6:108266/1‑149 (MQ=255)                                                                                                  aaGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGt                                              >  6:113607/1‑149 (MQ=255)                                                                                                       gTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑Cg                                                                                        <  8:267010/102‑1 (MQ=255)                                                                                                       gTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑Cg                                                                                        >  7:267010/1‑102 (MQ=255)                                                                                                           ggCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGt                                      >  4:349048/1‑148 (MQ=255)                                                                                                            gCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGt                                      >  1:217391/1‑147 (MQ=255)                                                                                                            gCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGt                                      >  5:22851/1‑147 (MQ=255)                                                                                                                        tcttcGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCt                         >  4:41490/1‑149 (MQ=255)                                                                                                                            cGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCg                                <  6:206620/138‑1 (MQ=255)                                                                                                                            cGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCg                                >  5:206620/1‑138 (MQ=255)                                                                                                                            cGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTg                       <  5:108266/147‑1 (MQ=255)                                                                                                                            cGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTg                       <  5:113607/147‑1 (MQ=255)                                                                                                                             gCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTGgcg                    >  5:244874/1‑149 (MQ=255)                                                                                                                             gCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTGgc                     <  7:25845/148‑1 (MQ=255)                                                                                                                               cACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGc                                                                                                                               >  2:177192/1‑41 (MQ=255)                                                                                                                                aCTTCATCAGGCCCTGTGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGCTCACCTTGGCGCTc                 <  7:325410/149‑1 (MQ=255)                                                                                                                                aCTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTGGCGCTc                 <  3:41490/149‑1 (MQ=255)                                                                                                                                         ggCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTGGCGCTCGATATCGGc        >  7:38705/1‑149 (MQ=255)                                                                                                                                          gCCCTGCGCACAACACATGACTGCGCAGGAGGGTCAGCCATACAGGGTTGCTGGTGAGCTCGACAT‑CGCCCAGCAACACCCTCGCTCTGCTCAGGGCAAGTATCGACACTATCAGCGTTAGCCCGTCACCTTGGCGCTCGATATCggat       >  6:177096/1‑147 (MQ=255)                                                                                                                                             cTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTGGCGCTCGATATCTGCTTCt    <  6:306134/149‑1 (MQ=255)                                                                                                                                             cTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTGGCGCTCGATATCGGCTTCt    <  3:349048/149‑1 (MQ=255)                                                                                                                                                cggcgATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCATGCGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTGGCGCTCGATATCGGCTTCTGa  >  1:298949/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     GCACCGACAACCTGGGCGCCGTATTTACCGGCGATGCGCCGGGCCTCGCCGAGGCCGTCGTATTTTTCTTCGATCAACAAGCTGACCAGATCGACATAAGTGGTCAGGCCATC‑CATGTTCTTCGCCACTTCATCAGGCCCTGCGGCGATCACATGACTGCGCTGGCGGGTCAGCCATACCGGGTTGCTGGTGCGGCCGCCAT‑CGGCCAGCCAGAGCGGGGCGCTGCGCTGGGCTTGTCTCGCCAGTTTCAGCGTTAGCCGGTCACCTTGGCGCTCGATATCGGCTTCTGA  >  NC_002947/4650435‑4650724

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