Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,650,014	A→T	8.6%	intergenic (+1491/+200)	PP_3214 → / ← PP_3216	hypothetical protein/cysteine dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,650,014	0	A	T	8.6%	159.7 / 5.6	58	intergenic (+1491/+200)	PP_3214/PP_3216	hypothetical protein/cysteine dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base T (2/3);  total (28/30)
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.06e-01

AAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCT  >  NC_002947/3649868‑3650159                                                                                                                                                     |                                                                                                                                                  aaaGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATc                                                                                                                                                 <  5:94147/149‑1 (MQ=255) aaaGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATc                                                                                                                                                 >  5:86639/1‑149 (MQ=255)             ggCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCTACGATGACGGTGa                                                                                                                                      <  5:175172/148‑1 (MQ=255)                aCGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCTACGATGACGGTGAggg                                                                                                                                   <  1:286024/148‑3 (MQ=255)                aCGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCTACGATGACGGTGAggg                                                                                                                                   >  2:286024/1‑146 (MQ=255)                 cGTTGACATCACGAAACAACTCGCCG‑CCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATc                                                                                                                                                 <  8:70537/132‑1 (MQ=255)                 cGTTGACATCACGAAACAACTCGCCG‑CCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCACTCACCGTCATc                                                                                                                                                 >  7:70537/1‑132 (MQ=255)                         tCACGAAACAACTCGCCGCCCGCAGCGCACCAGCCGGCAAGGCCGCTATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCccg                                                                                                                         <  6:210983/149‑1 (MQ=255)                             gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCgg                                                                                                                     >  3:298449/1‑149 (MQ=255)                             gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCgg                                                                                                                     >  2:28997/1‑149 (MQ=255)                              aaaCAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCggg                                                                                                                    >  7:225120/1‑149 (MQ=255)                                  aacTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAc                                                                                                                >  3:32290/1‑149 (MQ=255)                                  aacTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAc                                                                                                                >  3:179349/1‑149 (MQ=255)                                    cTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAcgc                                                                                                              >  4:232681/1‑149 (MQ=255)                                     tCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAcgcg                                                                                                             >  3:260309/1‑149 (MQ=255)                                      cgccgcCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCgg                                                                                                            >  1:279982/1‑149 (MQ=255)                                       gccgccCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAGACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGCGTGCGGGGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCgg                                                                                                            >  1:235913/1‑149 (MQ=255)                                        ccgcccgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGc                                                                                                          <  1:28997/149‑1 (MQ=255)                                        ccgcccgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCACCGGGGAACGCGGGc                                                                                                          >  3:89724/1‑149 (MQ=255)                                                aGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGAttt                                                                                                  <  4:32290/149‑1 (MQ=255)                                                 gcgcgcCAGCCGGCAAGGCCGCCATCAAGCAATGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTTCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTc                                                                                                 <  6:86345/149‑1 (MQ=255)                                                  cgcgcCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTcc                                                                                                <  8:191895/149‑1 (MQ=255)                                                     gcCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGc                                                                                             >  2:36064/1‑149 (MQ=255)                                                         gCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCa                                                                                         <  4:26848/149‑1 (MQ=255)                                                            ggCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGtt                                                                                      >  2:183131/1‑149 (MQ=255)                                                             gCAAGGCCGCCATCAAGCAAGGCA‑ACATTGAAGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGAGGTGTCCCGCCGCGGAACGCGGGCTAAGATTTCCAGCTCCAGttt                                                                                     >  4:30051/1‑149 (MQ=255)                                                              cAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCACCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTc                                                                                    <  4:89724/149‑1 (MQ=255)                                                              cAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGCCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTc                                                                                    <  5:278516/149‑1 (MQ=255)                                                                  gccgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACg                                                                                                                                          >  6:145931/1‑91 (MQ=255)                                                                  gccgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACg                                                                                                                                          <  5:145931/91‑1 (MQ=255)                                                                  gccgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACa                                                                                <  3:232681/149‑1 (MQ=255)                                                                  cgcgccATCAAGCAAGGCT‑ACATTGACTTAGCCCCAGGCCAGTAAACGCTGAGCGGCGTCTGCGGCGA‑GCCCCTCACCGGCATCGTAGAGTGTGAGTGCGGTGTCCCGCGGGGGAACGCGGGCGAAGATTTCCAGCTGAAGTTTCgg                                                                                  >  8:216217/3‑146 (MQ=255)                                                                    cgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAgg                                                                              >  1:54475/1‑149 (MQ=255)                                                                     gccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAggg                                                                             <  6:86639/149‑1 (MQ=255)                                                                        aTCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTa                                                                          >  8:252365/1‑149 (MQ=255)                                                                              caagGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTgg                                                                    <  1:183131/149‑1 (MQ=255)                                                                                 ggCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGGGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCcg                                                                 >  6:21361/1‑149 (MQ=255)                                                                                      aCATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGaa                                                             <  1:36064/149‑1 (MQ=255)                                                                                            aCGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCg                                                        >  7:54203/1‑148 (MQ=255)                                                                                                 gCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATgg                                                   <  3:74816/148‑1 (MQ=255)                                                                                                   ccAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGcc                                                >  2:312897/1‑149 (MQ=255)                                                                                                     aaGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTg                                              <  8:54203/149‑1 (MQ=255)                                                                                                      aGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCACCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTgg                                             <  5:309010/149‑1 (MQ=255)                                                                                                        gCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGc                                           <  2:279982/149‑1 (MQ=255)                                                                                                        gCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGc                                           <  7:252365/149‑1 (MQ=255)                                                                                                            gTAAACGCTGCGCGGCGACTGCGGCGACGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGa                                         >  1:14704/1‑148 (MQ=255)                                                                                                                    tGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGTGAACAGCGGATGGGCCTGGGCGAAAGGGtcttc                               >  4:187597/1‑149 (MQ=255)                                                                                                                      cgcgGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGtcttc                               <  1:164535/147‑1 (MQ=255)                                                                                                                           cgACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGt                        <  2:54475/149‑1 (MQ=255)                                                                                                                               ttcggcgA‑GCCCCTCACCGTCATCGTAGAAGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGt                    <  3:30051/147‑1 (MQ=255)                                                                                                                                  ggcgA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGTGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGa                 <  3:187597/149‑1 (MQ=255)                                                                                                                                   gcgA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGAt                <  8:225120/149‑1 (MQ=255)                                                                                                                                     gA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGACGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAATGGTCTTCTTCGCGTGAGTCGAt                >  2:304227/1‑147 (MQ=255)                                                                                                                                        gCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAgg            <  2:235913/149‑1 (MQ=255)                                                                                                                                           ccTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGcc         <  1:312897/149‑1 (MQ=255)                                                                                                                                           acTCACCGTCTACGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTa                       >  6:193638/2‑135 (MQ=255)                                                                                                                                           acTCACCGTCTACGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTa                       <  5:193638/134‑1 (MQ=255)                                                                                                                                                  gTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGctcct  >  7:101060/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AAAGCTTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCT  >  NC_002947/3649868‑3650159

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