Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,541,810	C→T	23.5%	D297N (GAC→AAC)	dnaB ←	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,541,810	0	C	T	23.5%	120.1 / 36.0	68	D297N (GAC→AAC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (22/30);  new base T (9/7);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

TCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCAGCGACTCACCTGGCATCTCGAGGGAGAA  >  NC_002947/5541662‑5541951                                                                                                                                                     |                                                                                                                                              tCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTc                                                                                                                                               <  3:20944/149‑1 (MQ=255)             cccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGtt                                                                                                                                   <  3:242832/149‑1 (MQ=255)             cccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCCATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGtt                                                                                                                                   <  4:182463/149‑1 (MQ=255)              ccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGtt                                                                                                                                   >  5:200877/1‑148 (MQ=255)              ccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGtt                                                                                                                                   >  6:318405/1‑148 (MQ=255)                      tGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGcc                                                                                                                            >  7:40312/1‑147 (MQ=255)                        cATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTg                                                                                                                              >  4:125512/1‑143 (MQ=255)                        cATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTg                                                                                                                              <  3:125512/143‑1 (MQ=255)                        cATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGa                                                                                                                          <  5:164773/147‑1 (MQ=255)                               tGCAGGTAGTCTACCTTGATCATGGCTTTTTCGCTATGCTCACGGGCCAGGTGGGGGGTACGCTCGCGCTTTTCCGAGGGGCTGATCCCGGCGTTATCGTCTATTAACAATTTTCGTTCGTTGAGCAGGTTCACTGCCGAGTTCTGgc                                                                                                                  <  8:107153/148‑1 (MQ=255)                                gCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgcgc                                                                                                                >  7:71963/1‑149 (MQ=255)                                 cAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgcgc                                                                                                                <  4:277842/148‑1 (MQ=255)                                 cAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGcgcg                                                                                                               >  2:205018/1‑149 (MQ=255)                                  aGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAg                                                                                                                         <  4:205107/138‑1 (MQ=255)                                  aGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAg                                                                                                                         >  3:205107/1‑138 (MQ=255)                                   ggTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGc                                                                                                             <  8:71963/149‑1 (MQ=255)                                       gTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAg                                                                                                          <  2:235789/148‑1 (MQ=255)                                           aCCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc                                                                                                     >  3:50993/1‑149 (MQ=255)                                           aCCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc                                                                                                     >  3:115514/1‑149 (MQ=255)                                           aCCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCa                                                                                                           <  4:220339/143‑1 (MQ=255)                                           aCCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCa                                                                                                           >  3:220339/1‑143 (MQ=255)                                              aTGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtc                                                                                                  <  2:102114/149‑1 (MQ=255)                                              aTGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtc                                                                                                  <  7:26472/149‑1 (MQ=255)                                                 ctCATGGCGAATTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtcgt                                                                                                <  7:301455/147‑1 (MQ=255)                                                  tCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtcgtc                                                                                               >  7:242493/1‑148 (MQ=255)                                                    aTGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCa                                                                                             <  3:122766/148‑1 (MQ=255)                                                        cGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTgg                                                                                        <  4:50993/149‑1 (MQ=255)                                                         gATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGc                                                                                       >  1:269750/1‑149 (MQ=255)                                                          aTTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTcc                                                                                              >  5:93480/1‑141 (MQ=255)                                                          aTTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTcc                                                                                              <  6:93480/141‑1 (MQ=255)                                                              cGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAg                                                                                                                                      <  7:272009/97‑1 (MQ=255)                                                              cGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAg                                                                                                                                      >  8:272009/1‑97 (MQ=255)                                                                  aTGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGa                                                                              >  4:209218/1‑149 (MQ=255)                                                                      tCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCgg                                                                                    >  5:254088/1‑139 (MQ=255)                                                                      tCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCgg                                                                                    <  6:254088/139‑1 (MQ=255)                                                                      tCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAAc                                                                                 <  2:22007/142‑1 (MQ=255)                                                                      tCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAAc                                                                                 >  1:22007/1‑142 (MQ=255)                                                                         cGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCACCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTtggt                                                                       <  2:288400/149‑1 (MQ=255)                                                                            gCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTtggtctg                                                                    >  8:16246/1‑149 (MQ=255)                                                                             ccAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTtggtctgg                                                                   <  5:318405/149‑1 (MQ=255)                                                                              cAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTtggtctggt                                                                  <  4:115514/149‑1 (MQ=255)                                                                               aGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCTGCTGGCCGGAACGGACCTtggtctggtc                                                                 <  1:205018/149‑1 (MQ=255)                                                                                   ggcggGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATa                                                             <  8:40312/149‑1 (MQ=255)                                                                                       ggTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGa                                                                              >  1:101419/1‑128 (MQ=255)                                                                                       ggTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGa                                                                              <  2:101419/128‑1 (MQ=255)                                                                                            gcgcgcgcATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAgg                                                    <  3:209218/149‑1 (MQ=255)                                                                                                  gcATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTtggtct                                                                     <  4:336207/126‑1 (MQ=255)                                                                                                  gcATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTtggtct                                                                     >  3:336207/1‑126 (MQ=255)                                                                                                      ttCAGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGACGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGc                                           >  2:342454/1‑148 (MQ=255)                                                                                                        ccGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATa                                        >  4:25951/1‑149 (MQ=255)                                                                                                         cGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGACCTGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATAc                                       >  5:121495/1‑149 (MQ=255)                                                                                                           aGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGc                                     <  7:16246/149‑1 (MQ=255)                                                                                                             gggCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCCTACGCAt                                   >  7:252088/1‑149 (MQ=255)                                                                                                               gCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGa                                 >  3:172136/1‑149 (MQ=255)                                                                                                                   aTACCCGCCGTATCGTCGATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCa                              <  6:200877/148‑1 (MQ=255)                                                                                                                       ccgccgTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCg                                                                >  2:57192/1‑110 (MQ=255)                                                                                                                       ccgccgTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCg                                                                <  1:57192/110‑1 (MQ=255)                                                                                                                       ccgccgTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCAGCGAc                         >  2:139381/1‑149 (MQ=255)                                                                                                                            gTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGACGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCAGCGACTCAcc                    <  1:342454/149‑1 (MQ=255)                                                                                                                                 gtcgATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAggg                                                   >  5:189741/1‑113 (MQ=255)                                                                                                                                 gtcgATGAACAATTTTCGGTTGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAggg                                                   <  6:189741/113‑1 (MQ=255)                                                                                                                                 gtcgATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCgg                                                                                    >  2:236482/1‑80 (MQ=255)                                                                                                                                 gtcgATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGACCTGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCAGCGACTCACCTGGc                <  6:121495/148‑1 (MQ=255)                                                                                                                                 gtcgATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGACAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCgg                                                                                    <  1:236482/80‑1 (MQ=255)                                                                                                                                     aTGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCAGCGACTCACCTGGCAtctc           <  2:269750/149‑1 (MQ=255)                                                                                                                                        aaCAATTTTCGGTCGTTGAGCAGGTTTACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCAGCGACTCACCTGGc                <  2:169281/141‑1 (MQ=255)                                                                                                                                        aaCAATTTTCGGTCGTTGAGCAGGTTTACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCAGCGACTCACCTGGc                >  1:169281/1‑141 (MQ=255)                                                                                                                                             ttttCGGTTGTTGAGCAGGTTCACTGCCGACGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCAGCGACTCACCTGGCATCTCGAGGgag    >  1:28821/1‑148 (MQ=255)                                                                                                                                               ttCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCAGCGACTCACCTGGCATCTCGAGGGAGaa  <  1:139381/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              TCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTGGCCGGAACGGACCTTGGTCTGGTCGATACGCCCCAGGGACGAAAGCATACGCATGATCAGCGACTCACCTGGCATCTCGAGGGAGAA  >  NC_002947/5541662‑5541951

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