Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,152,625	G→C	5.6%	R69P (CGA→CCA)	PP_5738 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,152,625	0	G	C	5.6%	232.5 / 4.1	72	R69P (CGA→CCA)	PP_5738	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base C (2/2);  total (38/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 = 8.21e-01

GATTGGCCAGCTTGCCTGGCCAAGTCCGCGTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGACCTCGCGCCAGAC  >  NC_002947/6152481‑6152773                                                                                                                                                 |                                                                                                                                                     gATTGGCCAGCTTGGCTGGCCAAGTCCGCGTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTc                                                                                                                                                  >  7:290958/1‑149 (MQ=255) gATTGGCCAGCTTGCCTGGCCAAGTCCGCGTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTc                                                                                                                                                  <  8:229612/149‑1 (MQ=255)       ccAGCTTGCCTGGCCAAGTCCGCGTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACTCGGTCGACatc                                                                                                                                            <  3:324218/149‑3 (MQ=255)              gCCTGGCCAAGTCCGCGTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGtt                                                                                                                                       <  5:125403/147‑1 (MQ=255)                  ggCCAAGTCCGCGTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGc                                                                                                                                 >  5:195462/1‑149 (MQ=255)                      aaGTCCGCGTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGtt                                                                                                                             <  4:205758/149‑1 (MQ=255)                        gTCCGCGTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGc                                                                                                                           >  5:322750/1‑149 (MQ=255)                           cgcgTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTGCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGt                                                                                                                         <  5:65076/148‑1 (MQ=255)                           cgcgTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGt                                                                                                                         >  8:177532/1‑148 (MQ=255)                            gcgTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCa                                                                                                                       <  3:144943/149‑1 (MQ=255)                                       tACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCCGCACTCaaaaa                                                                                                            >  8:168642/1‑149 (MQ=255)                                       tACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCaaaa                                                                                                             >  7:227277/1‑148 (MQ=255)                                        aCAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCaaaaaa                                                                                                           >  5:29546/1‑149 (MQ=255)                                        aCAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCaaaaa                                                                                                            >  1:122793/1‑148 (MQ=255)                                         cagcagAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCaaaaaa                                                                                                           <  4:76997/148‑1 (MQ=255)                                            cagAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAaagaag                                                                                                       <  2:57722/149‑1 (MQ=255)                                             agagagCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAgg                                                                                                      <  2:307451/149‑1 (MQ=255)                                                gagCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAggcg                                                                                                    <  8:227277/148‑1 (MQ=255)                                                    tGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCaa                                                                                                >  7:144813/1‑148 (MQ=255)                                                    tGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAAt                                                                                               >  3:319201/1‑149 (MQ=255)                                                       cACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAAtctc                                                                                            <  6:29546/149‑1 (MQ=255)                                                       cACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAAtctc                                                                                            <  2:107554/149‑1 (MQ=255)                                                        aCGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCtcg                                                                                           >  8:143321/1‑149 (MQ=255)                                                             tcgtTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCtcgtcgtc                                                                                      <  8:241753/149‑1 (MQ=255)                                                                 ttCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTCATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTTGTCGAct                                                                                  >  5:168634/1‑149 (MQ=255)                                                                   cTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGActct                                                                                >  6:755/1‑149 (MQ=255)                                                                    tGAGCGTGTAGGCCAAGTCACGCAAGACGATATGAGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGActct                                                                                >  1:226468/1‑148 (MQ=255)                                                                      aGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCAc                                                                             >  8:167859/1‑149 (MQ=255)                                                                            tAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGtt                                                                       <  6:195462/149‑1 (MQ=255)                                                                              ggCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTa                                                                      <  8:290958/148‑1 (MQ=255)                                                                               gCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTAc                                                                     >  6:43374/1‑148 (MQ=255)                                                                                ccAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGa                                                                   >  4:244290/1‑149 (MQ=255)                                                                                  aaGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATc                                                                 <  5:755/149‑1 (MQ=255)                                                                                       aCGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGaaa                                                             >  8:230169/1‑148 (MQ=255)                                                                                       aCGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGaaa                                                             >  6:231545/1‑148 (MQ=255)                                                                                          cAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACAtt                                                         <  7:143321/149‑1 (MQ=255)                                                                                              aCGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTa                                                                      >  7:15755/1‑132 (MQ=255)                                                                                              aCGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTa                                                                      <  8:15755/132‑1 (MQ=255)                                                                                                    tGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAgg                                               <  5:43374/149‑1 (MQ=255)                                                                                                       gagagaAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGCAGTTGATTGCGGTTGCGTCAGCACGCAAAAAAGAAGGCGGCACTCTCGTCGTCGACTCTCCCCCGGTTACGAACGAAACATTCGGTTACAGGGt                                             >  6:88177/1‑148 (MQ=255)                                                                                                             aGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCa                                                 >  6:349963/1‑138 (MQ=255)                                                                                                             aGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCa                                                 <  5:349963/138‑1 (MQ=255)                                                                                                              gTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCAt                                     <  4:319201/149‑1 (MQ=255)                                                                                                                 ctctgCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCAttt                                   <  5:88177/144‑1 (MQ=255)                                                                                                                 cGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCAttt                                   <  2:122793/148‑1 (MQ=255)                                                                                                                 cGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTg                                  <  3:244290/149‑1 (MQ=255)                                                                                                                       tctGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAg                            >  8:312917/1‑149 (MQ=255)                                                                                                                       tctGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAg                            >  5:124349/1‑149 (MQ=255)                                                                                                                              aCGCGCCAAGATGGCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGtt                     <  2:226468/149‑1 (MQ=255)                                                                                                                               cgcgCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTg                    >  1:240364/1‑149 (MQ=255)                                                                                                                                gcgcCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGt                                                     <  4:7058/115‑1 (MQ=255)                                                                                                                                gcgcCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGt                                                     >  3:7058/1‑115 (MQ=255)                                                                                                                                gcgcCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCCAAACATTCGGTTTCAGGGTGACGCCa                                      >  4:172481/1‑130 (MQ=255)                                                                                                                                    ctccctGTCGACTCGGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGa               >  7:104507/6‑149 (MQ=255)                                                                                                                                    cAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGaaa                                                             <  6:35840/103‑1 (MQ=255)                                                                                                                                    cAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGaaa                                                             >  5:35840/1‑103 (MQ=255)                                                                                                                                    cAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGa               <  2:240364/149‑1 (MQ=255)                                                                                                                                     aaGATGTCGCCCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGTGTGACGCCATTTGCAAAAGCGCAGTTGCATGAc              >  3:18761/1‑149 (MQ=255)                                                                                                                                     aaGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCTGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGAc              <  7:167859/149‑1 (MQ=255)                                                                                                                                       cctGTCGACTCGGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAgc                           <  3:83411/132‑1 (MQ=255)                                                                                                                                       cctGTCGACTCGGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAgc                           >  4:83411/3‑134 (MQ=255)                                                                                                                                        aTGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGACCTc           <  7:177532/149‑1 (MQ=255)                                                                                                                                           tCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCAt                 >  7:203356/1‑140 (MQ=255)                                                                                                                                           tCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGACCTcgc         <  6:322750/148‑1 (MQ=255)                                                                                                                                           gccaCCGAGTCTACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCAt                 <  8:203356/137‑1 (MQ=255)                                                                                                                                             gACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGACCTCgcgc       >  3:189193/1‑148 (MQ=255)                                                                                                                                             gACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGACCTCgcgc       >  2:466/1‑148 (MQ=255)                                                                                                                                             gACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGACCTCgcgc       >  2:212835/1‑148 (MQ=255)                                                                                                                                             gACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGACCTCgcgc       >  2:190707/1‑148 (MQ=255)                                                                                                                                              aCCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGCCTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGACCTCgcgc       >  4:116570/1‑147 (MQ=255)                                                                                                                                               ccGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCAttt                                   <  6:136858/118‑1 (MQ=255)                                                                                                                                               ccGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCAttt                                   >  5:136858/1‑118 (MQ=255)                                                                                                                                                 gAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGACCTCGCGCCAGAc  <  5:231545/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GATTGGCCAGCTTGCCTGGCCAAGTCCGCGTATTTACCTACAGCAGAGAGCTGGCACGCGTCGTTTCTGAGCGTGTAGGCCAAGTCACGCAAGACGATATGCGAGAGAAGTCCGCAGCTCTGATAACGCGCCAAGATGTCGACCGAGTCGACAGGGAGTTGATTGCGGTTGCGTCAGCACTCAAAAAAGAAGGCGGCAATCTCGTCGTCGACTCTCACCCGGTTACGATCGAAACATTCGGTTTCAGGGTGACGCCATTTGCAAAAGCGCAGTTGCATGACCTCGCGCCAGAC  >  NC_002947/6152481‑6152773

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