Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,125,301	A→T	7.4%	intergenic (+30/+51)	PP_3628 → / ← PP_3629	choline/carnitine/betaine transporter family protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,125,301	0	A	T	7.4%	199.3 / 5.1	67	intergenic (+30/+51)	PP_3628/PP_3629	choline/carnitine/betaine transporter family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/31);  new base T (2/3);  total (33/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.40e-01

ATCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTG  >  NC_002947/4125155‑4125447                                                                                                                                                    |                                                                                                                                                      aTCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTaa                                                                                                                                                      <  3:91293/148‑1 (MQ=255)       aaGGCCGTACTGTATGGCTCGGTGAAGCTTGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGccc                                                                                                                                               <  7:31642/149‑1 (MQ=255)           ccGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCG‑‑tgg                                                                                                                                         <  4:95114/149‑3 (MQ=255)                 tGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATg                                                                                                                                  <  2:278480/149‑1 (MQ=255)                      ggCGCGTTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                                 <  5:220531/132‑1 (MQ=255)                      ggCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                                 >  6:220531/1‑132 (MQ=255)                       gcgcGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGc                                                                                                                                 <  7:343973/144‑1 (MQ=255)                       gcgcGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGc                                                                                                                                 >  8:343973/1‑144 (MQ=255)                             tGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGcc                                                                                                                      >  8:262550/1‑149 (MQ=255)                                  cATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCctgct                                                                                                                  >  6:224542/1‑148 (MQ=255)                                      gggAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGCTGCGAGACCGCGCCTGCTGaaa                                                                                                              <  7:56209/148‑1 (MQ=255)                                             gTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGCCGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATc                                                                                                      <  6:126272/149‑1 (MQ=255)                                             gTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCACCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATc                                                                                                      >  5:327752/1‑149 (MQ=255)                                              ttCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑cgcg                                                                                                                                     >  2:131586/1‑117 (MQ=255)                                              ttCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑cgcg                                                                                                                                     <  1:131586/117‑1 (MQ=255)                                              ttCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCt                                                                                                     >  6:144557/1‑149 (MQ=255)                                                cTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGc                                                                                                   >  2:307350/1‑149 (MQ=255)                                                  aCATTGTGGTGCAGGAGAAGAAGGCATCGGCCAATACAGCGGCAGA‑GGCCGAAGTAACCCGAAACTGACGGCGTCAACCCGAAAACCACGGTCTGCTAAGAAGCCCG‑‑TT‑CGCGATTCGAGACCGAGCCTGCTGAAACATGGATCTgat                                                                                                  >  1:110616/1‑146 (MQ=255)                                                      tgtgGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAAc                                                                                             >  7:129930/1‑149 (MQ=255)                                                       gtggtgCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCTGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTg                                                                                                    >  3:29439/1‑141 (MQ=255)                                                       gtggtgCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTg                                                                                                    <  4:29439/141‑1 (MQ=255)                                                              aGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAc                                                                                     >  2:230341/1‑149 (MQ=255)                                                                  gaagaaGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgccg                                                                                 >  3:288426/1‑149 (MQ=255)                                                                      aagGCATCGGCCAAGACAGCGGCAGAGGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGCCCGCGCCTGCTGAAACAt                                                                                                           >  1:198511/1‑120 (MQ=255)                                                                      aagGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACAt                                                                                                           >  6:144662/1‑119 (MQ=255)                                                                      aagGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACAt                                                                                                           <  2:198511/119‑1 (MQ=255)                                                                      aagGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACAt                                                                                                           <  5:144662/119‑1 (MQ=255)                                                                      aagGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgccgccgc                                                                             <  8:129930/149‑1 (MQ=255)                                                                       agGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgccgccgc                                                                             >  7:144442/1‑148 (MQ=255)                                                                               gCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGc                                                                    <  6:327752/149‑1 (MQ=255)                                                                               gCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGc                                                                    <  7:262550/149‑1 (MQ=255)                                                                               gCCAAGACAGCGGCAGA‑GCCCGAAGTAACCCGCAACTGACGGCGTCAACCAGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGc                                                                    <  6:267878/149‑1 (MQ=255)                                                                                ccAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCctg                                                                                                                    >  7:78112/1‑100 (MQ=255)                                                                                ccAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCctg                                                                                                                    <  8:78112/100‑1 (MQ=255)                                                                                  aaGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTCTTACGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGc                                                                    <  1:307350/149‑1 (MQ=255)                                                                                      cAGCGGCAGA‑GGCCGAAGTAACCCGCCACTGACGGCGTCAACCAGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATc                                                              <  2:110616/148‑1 (MQ=255)                                                                                      cAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCg                                                             <  4:288426/149‑1 (MQ=255)                                                                                       aGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGa                                                            >  8:278047/1‑149 (MQ=255)                                                                                         cGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATc                                                          <  8:144442/149‑1 (MQ=255)                                                                                            cAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTg                                                       <  5:144557/149‑1 (MQ=255)                                                                                                        gTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGcat                                            <  5:224542/149‑1 (MQ=255)                                                                                                          aaCCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGTGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGcatca                                          >  3:351134/1‑149 (MQ=255)                                                                                                          aaCCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATgg                                                                                                         <  1:211223/86‑1 (MQ=255)                                                                                                          aaCCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATgg                                                                                                         >  2:211223/1‑86 (MQ=255)                                                                                                          aaCCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGcc                                                                   <  8:348314/124‑1 (MQ=255)                                                                                                          aaCCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGcc                                                                   >  7:348314/1‑124 (MQ=255)                                                                                                             ccGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTg                                                                      <  3:45942/118‑1 (MQ=255)                                                                                                             ccGCAACTGAAGGCGTCAACCAGAAAACCACGGTCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAAACTTCCAACGCCGACGCTATCttc                                                                      >  4:45942/1‑117 (MQ=255)                                                                                                                      aCGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGTGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAac                               <  4:351134/148‑1 (MQ=255)                                                                                                                      aCGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTc                                                                                         >  7:353516/1‑90 (MQ=255)                                                                                                                      aCGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTc                                                                                         <  8:353516/90‑1 (MQ=255)                                                                                                                      aCGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAac                               <  1:230341/148‑1 (MQ=255)                                                                                                                         gCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAacgcac                           >  2:119340/1‑149 (MQ=255)                                                                                                                           gTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAacgcacgc                         >  7:63142/1‑149 (MQ=255)                                                                                                                              aaCCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTg                      <  1:119340/149‑1 (MQ=255)                                                                                                                              aaCCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTg                      >  7:114798/1‑149 (MQ=255)                                                                                                                                 ccGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGgcc                   >  4:109537/1‑149 (MQ=255)                                                                                                                                  cGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTg                      >  7:312692/1‑145 (MQ=255)                                                                                                                                  cGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTg                      <  8:312692/145‑1 (MQ=255)                                                                                                                                     aaaCCACGGGCTTCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCTCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGccgccg               <  7:79353/149‑1 (MQ=255)                                                                                                                                        gaaCGGGCTGCTTAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATg                                        >  6:373124/3‑121 (MQ=255)                                                                                                                                        gaaCGGGCTGCTTAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATg                                        <  5:373124/119‑1 (MQ=255)                                                                                                                                         cACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACg           >  2:319098/1‑149 (MQ=255)                                                                                                                                              gctgctTAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGAt                                                           >  1:83949/1‑96 (MQ=255)                                                                                                                                              gctgctTAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGAt                                                           <  2:83949/96‑1 (MQ=255)                                                                                                                                                  ctAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGgccg                  >  7:348771/1‑133 (MQ=255)                                                                                                                                                  ctAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGgccg                  <  8:348771/133‑1 (MQ=255)                                                                                                                                                  ctAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTg  >  5:128335/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      ATCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGA‑GGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCG‑‑TT‑CGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTG  >  NC_002947/4125155‑4125447

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