Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,541,734	C→T	6.7%	R322H (CGT→CAT)	dnaB ←	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,541,734	0	C	T	6.7%	219.8 / 3.9	75	R322H (CGT→CAT)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base T (3/2);  total (40/35)
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.81e-01

CAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAACGGACCT  >  NC_002947/5541596‑5541878                                                                                                                                            |                                                                                                                                                   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                                                                                                                                  >  4:270318/1‑148 (MQ=255)          tCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTAc                                                                                                                                  <  3:270318/148‑1 (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                                                                                                              <  5:84029/135‑1 (MQ=255)                                           tcGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAgggg                                                                                                              >  6:84029/1‑135 (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        >  5:36642/1‑149 (MQ=255)                                                                                                                                  gCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAAc        >  1:57464/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  <  5:352071/149‑1 (MQ=255)                                                                                                                                        cTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAACGGACCt  >  3:98460/1‑149 (MQ=255)                                                                                                                                           aCGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑gtcgtcgtcgtc                      <  3:122267/126‑1 (MQ=255)                                                                                                                                           aCGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑gtcgtcgtcgtc                      >  4:122267/1‑126 (MQ=255)                                                                                                                                            |                                                                                                                                                   CAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCA‑GTCGTCGTCGTCCAGCTGGCCGGAACGGACCT  >  NC_002947/5541596‑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: ‑