Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,555,061	T→C	12.5%	L378P (CTG→CCG)	sbcB →	exodeoxyribonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,555,061	0	T	C	12.5%	157.8 / 13.0	64	L378P (CTG→CCG)	sbcB	exodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base C (4/4);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

TTGGCCAATCAACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAA  >  NC_002947/1554913‑1555200                                                                                                                                                     |                                                                                                                                            ttGGCCAATCAACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCt                                                                                                                                             >  8:304217/1‑149 (MQ=255)     ccAATCAACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGaa                                                                                                                                         <  6:207532/149‑1 (MQ=255)       aaTCAACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTGGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGACGGTTTTATAGGAGACCGCGACCGGCGCTTAAGCGAGCAAGTTCGTGCGCTGGAAcc                                                                                                                                       >  1:255970/1‑149 (MQ=255)       aaTCAACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCCggttc                                                                                                                                        >  5:35522/1‑145 (MQ=255)           aaaaGGAGCAATGGCAAGAAAAGCTGGAGAACATCTATGGCAAGGACGAGTTAGCCCATAGCGAGGACCCGGAACAGCAGTTTTATGAGGGAGTTATAGGTGACCGAGACCGTCTCTAATTCGATCAAATTCGTGCGCTAGAACCGgct                                                                                                                                   >  8:69621/4‑148 (MQ=255)            aCAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGgcg                                                                                                                                   >  5:287076/1‑148 (MQ=255)             cAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCCGGTTCCAGCGCa                                                                                                                                 <  6:35522/149‑1 (MQ=255)             cAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCCGGTTCCAGCGCa                                                                                                                                 <  7:286359/149‑1 (MQ=255)              aGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCATGCGCCGGTTCCAGCGCac                                                                                                                                <  1:213867/149‑2 (MQ=255)                gcgcAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGtt                                                                                                                              <  5:112207/149‑1 (MQ=255)                  gcAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGa                                                                                                                            <  5:210264/149‑1 (MQ=255)                    aaTGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGc                                                                                                                          <  7:50695/149‑1 (MQ=255)                        gCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACCGCAGTTGTATGAGGGTTTTATAGGCGACCGCGCCCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATg                                                                                                                      <  4:21074/149‑1 (MQ=255)                         cAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATgg                                                                                                                     >  7:173899/1‑149 (MQ=255)                          aaGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGc                                                                                                                    >  3:118730/1‑149 (MQ=255)                          aaGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGc                                                                                                                    <  4:38581/149‑1 (MQ=255)                                 gCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGa                                                                                                             <  7:304217/149‑1 (MQ=255)                                           aTCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCgatga                                                                                                   >  5:306486/1‑149 (MQ=255)                                             cTATGGCAAGGACGAGTTCGCCACGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACAGTCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGCCCAGAGGATGtt                                                                                                         >  2:236521/1‑141 (MQ=255)                                             cTATGGCAAGGACGAGATCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCATGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGtt                                                                                                         <  1:236521/141‑1 (MQ=255)                                                 ggCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCcgcg                                                                                             >  3:187998/1‑149 (MQ=255)                                                    aaGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGttt                                                                                          <  6:287076/149‑1 (MQ=255)                                                      ggACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGc                                                                                        >  3:244738/1‑149 (MQ=255)                                                      ggACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGc                                                                                        >  1:34092/1‑149 (MQ=255)                                                        aCGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCg                                                                                      <  6:306486/149‑1 (MQ=255)                                                           aGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGa                                                                                                                            <  6:125710/108‑1 (MQ=255)                                                           aGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGa                                                                                                                            >  5:125710/1‑108 (MQ=255)                                                                gCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCtgtt                                                                              >  4:301517/1‑149 (MQ=255)                                                                   ccGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTGTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTgttgtt                                                                           <  2:255970/149‑1 (MQ=255)                                                                    cgagcgagGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTc                                                                          <  4:138027/149‑1 (MQ=255)                                                                               ccGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCggg                                                                >  1:264870/1‑148 (MQ=255)                                                                                   aaCAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAAc                                                                                                                                        <  4:153520/72‑1 (MQ=255)                                                                                   aaCAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAAc                                                                                                                                        >  3:153520/1‑72 (MQ=255)                                                                                      agcagTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGATGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTacc                                                        >  8:71932/1‑147 (MQ=255)                                                                                          gTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                    <  4:89648/149‑1 (MQ=255)                                                                                             gTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCt                                                  >  3:18680/1‑148 (MQ=255)                                                                                              tATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTg                                                 >  7:27238/1‑148 (MQ=255)                                                                                              tATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGa                                                <  2:150269/149‑1 (MQ=255)                                                                                               aTGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGCGCCATGGCCAGTGGATGTTCGATGATCCGCGTTGGCCGGAGATGTTGTTCCGCTACCGGGCCGGTAGCTTccaggaa                                               >  8:54898/1‑144 (MQ=255)                                                                                                    ggTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCtt                                          >  5:159660/1‑149 (MQ=255)                                                                                                      ttttATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGa                                        <  4:24574/149‑1 (MQ=255)                                                                                                           taGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCgg                                    <  3:301517/148‑1 (MQ=255)                                                                                                            aGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGc                                   >  8:19038/1‑148 (MQ=255)                                                                                                               tGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCg                                  >  8:189435/1‑146 (MQ=255)                                                                                                                   cgcgACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCg                           >  8:299418/1‑149 (MQ=255)                                                                                                                     cgACCGGCGCTTATGCGAGCAAGGTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGc                         <  8:173899/149‑1 (MQ=255)                                                                                                                         cGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCAAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACg                     >  1:86139/1‑149 (MQ=255)                                                                                                                          ggCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAAc                                             <  3:170195/124‑1 (MQ=255)                                                                                                                          ggCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAAc                                             >  4:170195/1‑124 (MQ=255)                                                                                                                            cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCTAAGAGCGGCAACGTTg                  <  4:244738/149‑1 (MQ=255)                                                                                                                            cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                  <  4:118730/149‑1 (MQ=255)                                                                                                                            cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                  <  2:34092/149‑1 (MQ=255)                                                                                                                            cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                  <  7:189435/149‑1 (MQ=255)                                                                                                                            cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                  <  5:246594/149‑1 (MQ=255)                                                                                                                            cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGAGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                  <  7:299418/149‑1 (MQ=255)                                                                                                                             gcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTgg                 <  4:187998/149‑1 (MQ=255)                                                                                                                              cTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGt                <  4:18680/149‑1 (MQ=255)                                                                                                                                           tcaactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCg                           <  2:192472/120‑1 (MQ=255)                                                                                                                                           tcaactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGcc    >  6:135481/6‑148 (MQ=255)                                                                                                                                           tcaactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGcc    >  2:191317/6‑148 (MQ=255)                                                                                                                                           gTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                    >  1:279495/1‑100 (MQ=255)                                                                                                                                           gTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                    <  2:279495/100‑1 (MQ=255)                                                                                                                                            caactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCaa  >  7:230324/5‑149 (MQ=255)                                                                                                                                                tGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCt                                                  <  8:284523/97‑1 (MQ=255)                                                                                                                                                tGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCt                                                  >  7:284523/1‑97 (MQ=255)                                                                                                                                                     |                                                                                                                                            TTGGCCAATCAACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAA  >  NC_002947/1554913‑1555200

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