Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,176,699	T→C	6.3%	R460R (CGT→CGC)	PP_3676 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,176,699	0	T	C	6.3%	193.4 / 2.5	66	R460R (CGT→CGC)	PP_3676	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base C (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

ATACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATT  >  NC_002947/4176553‑4176846                                                                                                                                                   |                                                                                                                                                    aTACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGa                                                                                                                                                   <  5:170538/149‑1 (MQ=255)     cAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATaaa                                                                                                                                               >  7:275840/1‑149 (MQ=255)      aaGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATaaa                                                                                                                                               >  8:280145/1‑148 (MQ=255)      aaGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATaaa                                                                                                                                               <  7:280145/148‑1 (MQ=255)          gcgcGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                          <  7:16059/149‑1 (MQ=255)          gcgcGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                          <  3:95681/149‑1 (MQ=255)          gcgcGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                          <  7:292384/149‑1 (MQ=255)          gcgcGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                          <  7:30139/149‑1 (MQ=255)            gcgGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatg                                                                                                                                        <  4:129267/149‑3 (MQ=255)                      cAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCAt                                                                                                                              >  8:174707/1‑149 (MQ=255)                            cGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCcg                                                                                                                        >  1:329312/1‑149 (MQ=255)                                aaTCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGa                                                                                                                    <  5:16711/149‑1 (MQ=255)                                aaTCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGa                                                                                                                    <  8:275840/149‑1 (MQ=255)                                   cAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTg                                                                                                                 >  3:128354/1‑149 (MQ=255)                                          gcctAAATTATTTCTAACATAGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCAccc                                                                                                          >  2:188917/4‑149 (MQ=255)                                             tAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTAt                                                                                                       >  3:61797/1‑149 (MQ=255)                                                        tAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGa                                                                                            >  6:62951/1‑149 (MQ=255)                                                        tAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGa                                                                                            >  6:31317/1‑149 (MQ=255)                                                           cATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATg                                                                                          >  7:275354/1‑148 (MQ=255)                                                             tCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcac                                                                                       <  1:343877/149‑1 (MQ=255)                                                                aCATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGc                                                                                         >  7:205142/1‑144 (MQ=255)                                                                aCATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGc                                                                                         <  8:205142/144‑1 (MQ=255)                                                                    cAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAggg                                                                                 >  5:239708/1‑148 (MQ=255)                                                                    cAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGATCACAATATTCAGTATTCGATGAACAAAggg                                                                                 >  5:235989/1‑148 (MQ=255)                                                                      aaGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCAt                                                                              <  5:62951/149‑1 (MQ=255)                                                                        gggAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACGCTATTCGGTATTCGATGCACACAGGGCATTc                                                                            <  1:188917/149‑1 (MQ=255)                                                                        gggAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTc                                                                            <  7:174707/149‑1 (MQ=255)                                                                         ggAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTc                                                                            >  4:215287/1‑148 (MQ=255)                                                                          gAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTa                                                                          <  8:275354/149‑1 (MQ=255)                                                                           taCGGCGCGTTTCTCAAATTTCATATGTTAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCGGCTGGCATCAGCCGCGGAGTGCTCACCCTATTCCGTTTTCGATGCAGACAGTGCATTCTa                                                                          <  4:61797/147‑1 (MQ=255)                                                                               gcgcGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAAtt                                                                     >  7:28148/1‑149 (MQ=255)                                                                                 gcgTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAAtttt                                                                   >  7:5903/1‑149 (MQ=255)                                                                                    tttCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTAc                                                                >  6:208565/1‑149 (MQ=255)                                                                                       ctcAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACgg                                                              >  6:236361/1‑148 (MQ=255)                                                                                       ctcAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACgg                                                              >  1:239676/1‑148 (MQ=255)                                                                                        tcAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAgg                                                                                  >  8:73580/1‑127 (MQ=255)                                                                                        tcAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAgg                                                                                  <  7:73580/127‑1 (MQ=255)                                                                                            aTTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACgg                                                         >  8:29222/1‑148 (MQ=255)                                                                                                  tAGGTAAAATTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCAtt                                                  <  2:329312/149‑1 (MQ=255)                                                                                                     gTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCa                                                >  1:30069/1‑148 (MQ=255)                                                                                                     gTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCa                                                >  8:26848/1‑148 (MQ=255)                                                                                                     gTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCa                                                >  2:279534/1‑148 (MQ=255)                                                                                                      tAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCaa                                               <  5:236361/148‑1 (MQ=255)                                                                                                      tAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCaa                                               <  7:29222/148‑1 (MQ=255)                                                                                                          cTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGAGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGt                                          <  2:55694/149‑1 (MQ=255)                                                                                                             aaTTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAg                                        <  8:28148/148‑1 (MQ=255)                                                                                                               ttcttcCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcac                                                                                       >  7:220636/1‑99 (MQ=255)                                                                                                               ttcttcCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcac                                                                                       <  8:220636/99‑1 (MQ=255)                                                                                                                  ttcCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACt                                  >  6:298013/1‑149 (MQ=255)                                                                                                                     cAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTa                                <  5:208565/148‑1 (MQ=255)                                                                                                                     cAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTAc                               >  3:116085/1‑149 (MQ=255)                                                                                                                       aGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTc                                                 >  8:286107/1‑129 (MQ=255)                                                                                                                         ggTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCTCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTggg                           <  6:239708/149‑1 (MQ=255)                                                                                                                             ttACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGaaa                        <  8:5903/148‑1 (MQ=255)                                                                                                                              tACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAAt                       <  4:128354/148‑1 (MQ=255)                                                                                                                               aCGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATa                      >  6:170115/1‑148 (MQ=255)                                                                                                                                cGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAaga                                            <  1:301207/125‑1 (MQ=255)                                                                                                                                cGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAaga                                            >  2:301207/1‑125 (MQ=255)                                                                                                                                      ctCAAATTTATCGCGGTAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGcc                >  1:282147/2‑147 (MQ=255)                                                                                                                                      aTCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAg                                        <  4:125685/123‑1 (MQ=255)                                                                                                                                      aTCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAg                                        >  3:125685/1‑123 (MQ=255)                                                                                                                                       tCAAATTTATCACGGTAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGaaa                        >  1:282183/1‑138 (MQ=255)                                                                                                                                       tCAAATTTATCACGGTAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGaaa                        <  2:282183/138‑1 (MQ=255)                                                                                                                                            tttACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTa        <  2:239676/149‑1 (MQ=255)                                                                                                                                            tttACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTa        <  5:170115/149‑1 (MQ=255)                                                                                                                                                  gTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACGGGGAAATACATGCCAAGTTGTAGAGAtt  >  5:211865/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    ATACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATT  >  NC_002947/4176553‑4176846

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