Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,863,829	T→C	6.0%	T145T (ACT→ACC)	PP_1667 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,863,829	0	T	C	6.0%	194.2 / 2.1	66	T145T (ACT→ACC)	PP_1667	transporter
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base C (2/2);  total (32/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 = 9.09e-01

GCTTGAGCATGTCGCGCTGCCGTGGGTACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTGATGGGCCTGG  >  NC_002947/1863685‑1863977                                                                                                                                                  |                                                                                                                                                      gCTTGAGCATGTCGCGCTGCCGTGGGTACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGc                                                                                                                                                     <  4:310554/147‑1 (MQ=255) gCTTGAGCATGTCGCGCTGCCGTGGGTACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGc                                                                                                                                                     <  6:7593/147‑1 (MQ=255) gCTTGAGCATGTCGCGCTGCCGTGGGTACAGGGCCGCCTTGGACTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACAGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCt                                                                                                                                                   >  7:261749/1‑149 (MQ=255)     gAGCATGTCGCGCTGCCGTGGGTACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCaa                                                                                                                                                <  4:390925/148‑1 (MQ=255)      aGCATGTCGCGCTGCCGTGGGTACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGc                                                                                                                                              <  5:103769/149‑1 (MQ=255)      aGCATGTCGCGCTGCCGTGGGTACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGc                                                                                                                                              <  8:304837/149‑1 (MQ=255)        cATGTCGCGCTGCCGTGGGTACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCt                                                                                                                                            >  1:180067/1‑149 (MQ=255)              gcgcTGCCGTGGGTACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGtt                                                                                                                                      >  7:358210/1‑149 (MQ=255)               cgcTGCCGGGGGAACTGGGCGGCCTGGCCCTGGCCTATGGTTTCTGCAAA‑TGCCACACGATCACGCCCGCCATCGGCGCGCACATGGGCCACACCACCGATACCGTCGGAATGCTGCGGCCCCCGGCCACTCCCTAAAGCCGGGCGtgg                                                                                                                                     <  6:349422/149‑3 (MQ=255)                           tACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTg                                                                                                                         >  7:92926/1‑149 (MQ=255)                           tACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTg                                                                                                                         >  4:179369/1‑149 (MQ=255)                           tACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACCGCCTCAAGCCTGGCGTTGGTGGCCTGGTTg                                                                                                                         >  1:322157/1‑149 (MQ=255)                           tACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACCGCCTCAAGCCTGGCGTTGGTGGCCTGGTTg                                                                                                                         >  5:37199/1‑149 (MQ=255)                              aGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGc                                                                                                                       >  4:236842/1‑148 (MQ=255)                              aGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGc                                                                                                                       >  2:38524/1‑148 (MQ=255)                                ggCCGCCTTGGCCTGGCCGATGGTTTCTGGAAATTTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGCCCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGcc                                                                                                                      >  1:48253/1‑148 (MQ=255)                                ggCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCaa                                                                                                                    >  2:65366/1‑149 (MQ=255)                                 gccgccTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCaa                                                                                                                    >  5:140042/1‑148 (MQ=255)                                  ccgccTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGACTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAata                                                                                                                  <  8:92926/149‑1 (MQ=255)                                         ggcctggccGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGc                                                                                                            <  4:156885/148‑1 (MQ=255)                                          gcctggccGATGGTTTCTGGAAA‑TTCGACAAGAGCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCt                                                                                                           <  3:236842/148‑1 (MQ=255)                                            ctggccGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGAt                                                                                                        >  4:31241/1‑149 (MQ=255)                                               gccGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATccc                                                                                                     >  5:267273/1‑149 (MQ=255)                                               gccGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCATCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATccc                                                                                                     >  8:168512/1‑149 (MQ=255)                                                 cGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGTCTGGGGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑gg                                                                                                  <  8:261749/149‑1 (MQ=255)                                                 cGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑gg                                                                                                  <  6:871/149‑1 (MQ=255)                                                    tGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAg                                                                                               <  1:38524/149‑1 (MQ=255)                                                             gAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTc                                                                                       >  2:136488/1‑148 (MQ=255)                                                                     gACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATctgct                                                                               >  1:56477/1‑149 (MQ=255)                                                                      aCAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTg                                                                                             >  5:279346/1‑134 (MQ=255)                                                                      aCAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTg                                                                                             <  6:279346/134‑1 (MQ=255)                                                                        aaGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTgcg                                                                            >  5:7745/1‑149 (MQ=255)                                                                          gATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTgcgc                                                                           >  7:133730/1‑148 (MQ=255)                                                                          gATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTgcgc                                                                           >  1:127406/1‑148 (MQ=255)                                                                           aTCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCTAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGa                                                                         >  5:112751/1‑149 (MQ=255)                                                                              aaGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGAtt                                                                       <  6:267273/148‑1 (MQ=255)                                                                              aaGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTg                                                                      <  1:180771/149‑1 (MQ=255)                                                                              aaGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACCGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGAtt                                                                       <  6:37199/148‑1 (MQ=255)                                                                              aaGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACCGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGAtt                                                                       <  2:322157/148‑1 (MQ=255)                                                                               aGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTgg                                                                     <  8:358210/149‑1 (MQ=255)                                                                                ggCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTg                                                                                                             >  2:236596/1‑109 (MQ=255)                                                                                ggCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTg                                                                                                             <  1:236596/109‑1 (MQ=255)                                                                                ggCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTggg                                                                    <  1:136488/149‑1 (MQ=255)                                                                                ggCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTggg                                                                    <  4:122780/149‑1 (MQ=255)                                                                                   cgccATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACt                                                                 >  2:345066/1‑149 (MQ=255)                                                                                     ccATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCCAGGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACtt                                                                <  1:65366/149‑1 (MQ=255)                                                                                           gcgcgcACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGatgatg                                                         <  2:127406/149‑1 (MQ=255)                                                                                              cgcACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGc                                                      <  7:168512/149‑1 (MQ=255)                                                                                                 acaTGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAg                                                   <  4:77600/149‑1 (MQ=255)                                                                                                       ggCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGa                                                                   <  1:399525/127‑1 (MQ=255)                                                                                                       ggCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGa                                                                   >  2:399525/1‑127 (MQ=255)                                                                                                         cAGACCACCGATATCGTCGTCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCgg                                           <  1:314041/149‑1 (MQ=255)                                                                                                           gACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGcc                                         <  6:146562/149‑1 (MQ=255)                                                                                                              caccGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCctgc                                      >  1:303945/1‑149 (MQ=255)                                                                                                              caccGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCctgc                                      >  3:102418/1‑149 (MQ=255)                                                                                                               accGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCctgct                                     >  6:372429/1‑149 (MQ=255)                                                                                                                                       tCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAg              <  2:303945/148‑1 (MQ=255)                                                                                                                                        cccAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGt             >  3:378646/1‑148 (MQ=255)                                                                                                                                         ccAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTg            <  2:180067/148‑1 (MQ=255)                                                                                                                                         ccAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTGa           >  5:41278/1‑149 (MQ=255)                                                                                                                                          ccggCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATAGGGTGCTGATCCC‑GGTAGTCGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTGAt          <  8:106663/147‑1 (MQ=255)                                                                                                                                           aGGCCACTGACTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTGATg         >  4:376740/1‑149 (MQ=255)                                                                                                                                            ggCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTGATgg        <  1:345066/149‑1 (MQ=255)                                                                                                                                            ggCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTGATgg        <  5:372429/149‑1 (MQ=255)                                                                                                                                             gCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTGATggg       >  6:165295/1‑149 (MQ=255)                                                                                                                                             gCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTGATggg       <  6:140042/149‑1 (MQ=255)                                                                                                                                                  tGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTGATGGGCCTgg  <  6:41278/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      GCTTGAGCATGTCGCGCTGCCGTGGGTACAGGGCCGCCTTGGCCTGGCCGATGGTTTCTGGAAA‑TTCGACAAGATCAAGGCCGCCATCGGCGCGCACATGGGGCAGACCACCGATATCGTCGGCATGCTGCTGTCCCAGGCCACTGCCTCAAGCCTGGCGTTGGTGGCCTGGTTGGCCAATATGGTGCTGATCCC‑GGTAGTGGGCTTCTATCTGCTGCGCGATTGGGACTTGATGATGGGCAAGCTGCGCGGCCTGCTGCCGCGCCAGCGCGAGCCGCAAGTGATGGGCCTGG  >  NC_002947/1863685‑1863977

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