Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,176,701	A→G	6.1%	D461G (GAT→GGT)	PP_3676 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,176,701	0	A	G	6.1%	192.8 / 2.0	65	D461G (GAT→GGT)	PP_3676	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

ATACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGT  >  NC_002947/4176553‑4176849                                                                                                                                                     |                                                                                                                                                     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:31317/1‑149 (MQ=255)                                                        tAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGa                                                                                               >  6:62951/1‑149 (MQ=255)                                                           cATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATg                                                                                             >  7:275354/1‑148 (MQ=255)                                                             tCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcac                                                                                          <  1:343877/149‑1 (MQ=255)                                                                aCATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGc                                                                                            <  8:205142/144‑1 (MQ=255)                                                                aCATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGc                                                                                            >  7:205142/1‑144 (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                                                  <  7:29222/148‑1 (MQ=255)                                                                                                      tAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCaa                                                  <  5:236361/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                                               >  2:301207/1‑125 (MQ=255)                                                                                                                                cGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAaga                                               <  1:301207/125‑1 (MQ=255)                                                                                                                                      ctCAAATTTATCGCGGTAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGcc                   >  1:282147/2‑147 (MQ=255)                                                                                                                                      aTCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAg                                           >  3:125685/1‑123 (MQ=255)                                                                                                                                      aTCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAg                                           <  4:125685/123‑1 (MQ=255)                                                                                                                                       tCAAATTTATCACGGTAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGaaa                           <  2:282183/138‑1 (MQ=255)                                                                                                                                       tCAAATTTATCACGGTAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGaaa                           >  1:282183/1‑138 (MQ=255)                                                                                                                                            tttACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTa           <  2:239676/149‑1 (MQ=255)                                                                                                                                            tttACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTa           <  5:170115/149‑1 (MQ=255)                                                                                                                                                  gTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACGGGGAAATACATGCCAAGTTGTAGAGAtt     >  5:211865/1‑149 (MQ=255)                                                                                                                                                     aTAAATTTGAGGCCGCTTGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTg   <  6:235989/148‑1 (MQ=255)                                                                                                                                                     aTAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGt  <  2:30069/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGTTGTAGAGATTTGT  >  NC_002947/4176553‑4176849

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