Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,541,720	T→A	8.8%	I327F (ATC→TTC)	dnaB ←	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,541,720	0	T	A	8.8%	217.1 / 8.3	77	I327F (ATC→TTC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base A (4/3);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

GCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTG  >  NC_002947/5541573‑5541864                                                                                                                                                     |                                                                                                                                                  gCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAtt                                                                                                                                                  >  5:53705/1‑149 (MQ=255)  cTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAttt                                                                                                                                                 >  2:205456/1‑149 (MQ=255)      gAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGcc                                                                                                                                             <  5:4827/149‑1 (MQ=255)       aaaGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGTATCTTCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCa                                                                                                                                            <  4:170264/149‑1 (MQ=255)       aaaGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCa                                                                                                                                            <  3:15520/149‑1 (MQ=255)        aaGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCAt                                                                                                                                           <  7:127704/149‑1 (MQ=255)        aaGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCAt                                                                                                                                           <  1:90436/149‑1 (MQ=255)           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)                              aaTTCTTTGGCCAGGGCCTTGAGCGATCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTc                                                                                                                                                <  4:86876/122‑1 (MQ=255)                              aaTTCTTTGGCCAGGGCCTTGAGCGAGCGTGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTc                                                                                                                                                >  3:86876/1‑122 (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                                                                                                     >  6:346463/1‑148 (MQ=255)                                               cTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTcc                                                                                                     >  1:60526/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                                                                                                                         >  6:195828/1‑105 (MQ=255)                                                                      aaaTTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAggcggc                                                                                                                         <  5:195828/105‑1 (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                                                                        >  8:275090/1‑136 (MQ=255)                                                                                        ttGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGa                                                                        <  7:275090/136‑1 (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‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtc                 <  8:32760/148‑1 (MQ=255)                                                                                                                                      ccgtGAGCATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcg             >  1:99937/4‑149 (MQ=255)                                                                                                                                       cATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgt            >  2:37610/1‑149 (MQ=255)                                                                                                                                         tGAGCATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAAc                                                                   >  4:133700/1‑92 (MQ=38)                                                                                                                                         tGAGCATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAAc                                                                   <  3:133700/92‑1 (MQ=38)                                                                                                                                           agcATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtcgtc        <  2:99937/147‑1 (MQ=255)                                                                                                                                            tCATGGCTATTTCGCCATGCTCAGGGGCCAGCCGGCGGGTACGCGCGCGCATTTCGGTGGGGCTGATA‑CGCGCCGTATCGTCTATGAAc                                                                   <  4:178456/89‑1 (MQ=255)                                                                                                                                            tCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAAc                                                                   >  3:178456/1‑89 (MQ=255)                                                                                                                                            tCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTcc       >  6:352071/1‑149 (MQ=255)                                                                                                                                            gcATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtcgtc        >  6:309955/2‑148 (MQ=255)                                                                                                                                                 gCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTg  >  1:179298/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCCTTGAGCGAGCGGGAAATCTCGGAAATTTC‑GTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATA‑CCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTCCAGCTG  >  NC_002947/5541573‑5541864

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