Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,541,730	C→A	8.9%	E323D (GAG→GAT)	dnaB ←	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,541,730	0	C	A	8.9%	228.0 / 9.8	75	E323D (GAG→GAT)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (35/32);  new base A (4/3);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

GGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAACGGACCT  >  NC_002947/5541583‑5541878                                                                                                                                                     |                                                                                                                                                       ggCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCt                                                                                                                                                       <  1:307797/149‑1 (MQ=255)   cAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCa                                                                                                                                                     <  5:16080/149‑1 (MQ=255)   cAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCa                                                                                                                                                     >  7:155726/1‑149 (MQ=255)     aTGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAAATCGCCATGATCatg                                                                                                                                                   <  2:69172/149‑3 (MQ=255)              cAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCCGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAggcg                                                                                                                                          <  3:165143/149‑1 (MQ=255)                gTTGTATTCTTTGGCCTGGGCCTAGTGCGTGCTGTAAATCTCGGAAATTTC‑GTTGGTACGGTTGTGACCGTCCGACCCGGGATTCTGCATCAGCTGCTGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAggcggc                                                                                                                                        <  7:141794/149‑1 (MQ=255)                 ttGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGgcggcg                                                                                                                                       >  1:270778/1‑149 (MQ=255)                   gAATTCTTTGGCCAGGGCCTTGTGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                                      <  2:326303/148‑1 (MQ=255)                   gAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                                      <  2:270778/148‑1 (MQ=255)                   gAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCACGCGGCggg                                                                                                                                     >  7:259419/1‑149 (MQ=255)                       tCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTAc                                                                                                                                  <  3:270318/148‑1 (MQ=255)                       tCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTAc                                                                                                                                  >  4:270318/1‑148 (MQ=255)                        cTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTAcgc                                                                                                                                >  7:303884/1‑149 (MQ=255)                              ccAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTAcgcgcgcg                                                                                                                           >  7:203161/1‑148 (MQ=255)                               cAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTCGGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                                      >  1:238959/1‑137 (MQ=255)                               cAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                                      <  2:238959/136‑1 (MQ=255)                               cAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCa                                                                                                                         >  3:117342/1‑149 (MQ=255)                                     cTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTcc                                                                                                                    >  1:60526/1‑148 (MQ=255)                                     cTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTcc                                                                                                                    >  6:346463/1‑148 (MQ=255)                                        gagcgagcGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAgg                                                                                                                >  3:268658/1‑149 (MQ=255)                                          gcgagcgGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAgggg                                                                                                              >  7:32760/1‑149 (MQ=255)                                                 ggAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑c                                                                                                      <  8:203161/149‑1 (MQ=255)                                                  gAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACAATGATCATGGCGAATTCGCCATGCTCACGGGCCAGGCGCAGGGTACTCGCGCGCATTTCCTAGTGTCTGATA‑c                                                                                                      >  5:156643/1‑148 (MQ=255)                                                        tcGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAgggg                                                                                                              >  6:84029/1‑135 (MQ=255)                                                        tcGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAgggg                                                                                                              <  5:84029/135‑1 (MQ=255)                                                         cGGAAATTTC‑GTTGGTACGGTTGTCACCGGCAGACCCCGGAAACTGAAAAAGCTGCAGATAGACGACCATAATAATGGCAATTTAGCCATGATCACGGGCCAGGAGGCGGGTACGCGCGCGCATTTCAGAGGGGCTGAAA‑AACGCCGt                                                                                               >  7:224863/1‑148 (MQ=255)                                                           gAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTAtc                                                                                            <  7:269517/149‑1 (MQ=255)                                                            aaaTTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAggcggc                                                                                                                                        <  5:195828/105‑1 (MQ=255)                                                            aaaTTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAggcggc                                                                                                                                        >  6:195828/1‑105 (MQ=255)                                                                 tC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGa                                                                                       <  6:53705/148‑1 (MQ=255)                                                                 tC‑GTTGGTACGGTTGGCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGa                                                                                       <  6:221177/148‑1 (MQ=255)                                                                  c‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGAt                                                                                      >  2:3139/1‑148 (MQ=255)                                                                             gTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTc                                                                           <  8:155726/149‑1 (MQ=255)                                                                              ttGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGa                                                                                       <  7:275090/136‑1 (MQ=255)                                                                              ttGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGa                                                                                       >  8:275090/1‑136 (MQ=255)                                                                              ttGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCg                                                                          <  4:268658/149‑1 (MQ=255)                                                                                 tCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTc                                                                       <  4:117342/149‑1 (MQ=255)                                                                                 gcACCGGCCGCCCCCGGAATCTGCATCAGCTGCTGGTAGTCGACCATGATCATGGCGATTTCGCGATGGTCACGGGCGAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTc                                                                       <  8:259419/148‑1 (MQ=255)                                                                                    ccGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATGGTCGATGAACAATTTTCGGTCGtt                                                                    <  2:60526/149‑1 (MQ=255)                                                                                    ccGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGtt                                                                    <  5:346463/149‑1 (MQ=255)                                                                                     cGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGtt                                                                    >  5:86899/1‑148 (MQ=255)                                                                                      ggCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGtt                                                                    >  3:143786/1‑147 (MQ=255)                                                                                         cGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAg                                                                 >  4:14222/1‑147 (MQ=255)                                                                                          gACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAg                                                              >  8:2079/1‑149 (MQ=255)                                                                                             cccGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGtt                                                           >  3:244565/1‑149 (MQ=255)                                                                                               cGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCa                                                         >  8:204147/1‑149 (MQ=255)                                                                                                 gAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCCGCCGTATCGTCGATGAACAAtttt                                                                            <  1:185190/129‑1 (MQ=255)                                                                                                 gAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAAtttt                                                                            >  2:185190/1‑128 (MQ=255)                                                                                                  aaTCTGCATCAGCTGGAGGTAGTCGACCATGATCATGGCGATTTCGAACTGCTGACGGCCCAGGCGGCGGGTACGCGCGCGCATTTGCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTg                                                      >  4:53495/1‑149 (MQ=255)                                                                                                        cATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGcc                                                    <  8:303884/145‑1 (MQ=255)                                                                                                        cATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGa                                                  <  7:204147/147‑1 (MQ=255)                                                                                                        cATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAgg                                                <  6:86899/149‑1 (MQ=255)                                                                                                        cATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAg                                                 <  3:14222/148‑1 (MQ=255)                                                                                                             gCTGCATGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAg                                                              <  3:265131/130‑1 (MQ=255)                                                                                                             gCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAg                                                              >  4:265131/1‑130 (MQ=255)                                                                                                              cTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGTGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgc                                          <  5:99818/149‑1 (MQ=255)                                                                                                                tcatggaGTCGACCATGATCATGGCGATTTCGCCATGCGCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgc                                          <  4:341560/141‑1 (MQ=255)                                                                                                                  aGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATATTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgcgag                                       >  3:341560/1‑146 (MQ=255)                                                                                                                    gTAGTCGACCATGATCATGGCCATTTCGCCATCCTCACGGGCCAGGCGGCGGGTACCCGCGCGCATTTCCGAGGGGCTCATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCCTTGAGCAGGTTCCCTGCCGACGTCAGGCGCGGc                                     >  6:355709/1‑148 (MQ=255)                                                                                                                         cGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑gtc                               <  8:32760/148‑1 (MQ=255)                                                                                                                            ccgtGAGCATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑gtcgtcg                           >  1:99937/4‑149 (MQ=255)                                                                                                                             cATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑gtcgtcgt                          >  2:37610/1‑149 (MQ=255)                                                                                                                               tGAGCATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAAc                                                                                  >  4:133700/1‑92 (MQ=38)                                                                                                                               tGAGCATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAAc                                                                                  <  3:133700/92‑1 (MQ=38)                                                                                                                                 agcATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑gtcgtcgtcgtc                      <  2:99937/147‑1 (MQ=255)                                                                                                                                  tCATGGCTATTTCGCCATGCTCAGGGGCCAGCCGGCGGGTACGCGCGCGCATTTCGGTGGGGCTGATA‑CGCGCCGTATCGTCTATGAAc                                                                                  <  4:178456/89‑1 (MQ=255)                                                                                                                                  tCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAAc                                                                                  >  3:178456/1‑89 (MQ=255)                                                                                                                                  tCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTcc                     >  6:352071/1‑149 (MQ=255)                                                                                                                                  gcATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑gtcgtcgtcgtc                      >  6:309955/2‑148 (MQ=255)                                                                                                                                       gCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTg                >  1:179298/1‑149 (MQ=255)                                                                                                                                               gCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAAc        >  1:57464/1‑149 (MQ=255)                                                                                                                                               gCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAAc        >  5:36642/1‑149 (MQ=255)                                                                                                                                                ccATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCACGTCGTCGTCGTCCAGCTGGCCGGAAc        <  1:3139/149‑1 (MQ=255)                                                                                                                                                  aTGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAACGGa     <  1:205456/149‑1 (MQ=255)                                                                                                                                                    gCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAACGGAcc   >  5:301886/1‑149 (MQ=255)                                                                                                                                                     cTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAACGGACCt  >  3:98460/1‑149 (MQ=255)                                                                                                                                                     cTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAACGGACCt  <  5:352071/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAACGGACCT  >  NC_002947/5541583‑5541878

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑