Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,376,321	A→T	11.9%	intergenic (+30/+19)	kup → / ← PP_1201	potassium transport system protein Kup/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,376,321	0	A	T	11.9%	194.4 / 16.2	76	intergenic (+30/+19)	kup/PP_1201	potassium transport system protein Kup/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base T (5/4);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

CTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTGC  >  NC_002947/1376180‑1376468                                                                                                                                               |                                                                                                                                                    cTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGc                                                                                                                                              <  4:53515/149‑1 (MQ=255) cTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGc                                                                                                                                              <  2:279757/149‑1 (MQ=255)    gCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCt                                                                                                                                           <  2:56982/149‑1 (MQ=255)      ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGGACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                         <  8:208236/149‑1 (MQ=255)      ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                         <  5:171329/149‑1 (MQ=255)      ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTAGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                         <  8:77186/149‑1 (MQ=255)           aaCCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCtt                                                                                                                                    >  3:283633/1‑149 (MQ=255)                    gcgTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTCAGCCTAGTCACTTCt                                                                                                                           >  8:172370/1‑149 (MQ=255)                     cgTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGcc                                                                                                                                             >  7:50102/1‑130 (MQ=255)                     cgTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGcc                                                                                                                                             <  8:50102/130‑1 (MQ=255)                        tCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGAGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCt                                                                                                                        >  4:297194/1‑148 (MQ=255)                              tGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                         >  4:34942/1‑125 (MQ=255)                              tGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                         <  3:34942/125‑1 (MQ=255)                               gaagaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGGCAGTTCTGCTCAGCTACg                                                                                                                >  6:240091/1‑149 (MQ=255)                                aataaTGCCAACGGCAACTTGCTCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTAAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACg                                                                                                                <  8:132168/145‑1 (MQ=255)                                aagaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGc                                                                                                               >  8:238137/1‑149 (MQ=255)                                aagaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGc                                                                                                               >  3:65469/1‑149 (MQ=255)                                  gaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTg                                                                                                             <  5:194950/149‑1 (MQ=255)                                     tGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGAc                                                                                                          <  4:243536/149‑1 (MQ=255)                                       ccAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACtt                                                                                                        >  1:124294/1‑149 (MQ=255)                                       ccAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACtt                                                                                                        <  4:283633/149‑1 (MQ=255)                                        cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTg                                                                                                       >  5:275672/1‑149 (MQ=255)                                        cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCACGGCCCCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTAAGCTGGACTTg                                                                                                       >  4:175728/1‑149 (MQ=255)                                         aaCGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGc                                                                                                      >  4:210030/1‑149 (MQ=255)                                         aaCCGCAACTTGCGCGTCTTCAACCGGCCGTAGAACCGGGGGACCGAGCTAGGGA‑CCCAGGTCGAGATCTGACTGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGc                                                                                                      <  8:279920/149‑1 (MQ=255)                                              cAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAg                                                                                                                     <  3:56637/129‑1 (MQ=255)                                              cAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAg                                                                                                                     >  4:56637/1‑129 (MQ=255)                                                  ttGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGcgcg                                                                                             >  6:240442/1‑149 (MQ=255)                                                    gcgcTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGAGGCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGATCAGGGCCGCAGAGCGGCCGTGAGCCTTGTCAGTTGTGCTCAGCTACGCTGGACTTGCCTTGGCGCGt                                                                                            >  1:217982/1‑149 (MQ=255)                                                     cgcTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGttt                                                                                          <  2:217982/149‑1 (MQ=255)                                                        ttcttcAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGa                                                                                       <  6:275672/149‑1 (MQ=255)                                                              aaCCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCATTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGc                                                                                 >  5:271728/1‑149 (MQ=255)                                                              aaCCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCg                                                                                  <  7:172370/148‑1 (MQ=255)                                                              aaCCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGc                                                                                 <  3:210030/149‑1 (MQ=255)                                                               aCCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑GCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGTGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTTCCTTGGCGCGTTTCTATCTCGcc                                                                                <  1:140697/149‑1 (MQ=255)                                                               aCCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCg                                                                                        >  4:200583/1‑141 (MQ=255)                                                               aCCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCg                                                                                        <  3:200583/141‑1 (MQ=255)                                                                     cGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTGAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCaa                                                                          >  5:232986/1‑149 (MQ=255)                                                                      gTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAAc                                                                         <  4:65469/149‑1 (MQ=255)                                                                       ttGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAAcc                                                                        >  7:17979/1‑149 (MQ=255)                                                                       ttGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTAAGTTCTGCTCAGCAACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAAcc                                                                        >  8:118966/1‑149 (MQ=255)                                                                          aaCCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCAGAGCCTTGTCAGTTATGCTCAGCTACGCTGGACTTGACTTGGCGCGTTTCGATCTCGCCAATCAAACGCt                                                                     >  3:47906/1‑149 (MQ=255)                                                                           aCCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCtt                                                                    >  1:41750/1‑149 (MQ=255)                                                                           aCCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCtt                                                                    <  8:17979/149‑1 (MQ=255)                                                                             cGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGc                                                                                 >  2:206914/1‑134 (MQ=255)                                                                             cGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGc                                                                                 <  1:206914/134‑1 (MQ=255)                                                                               ggTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGc                                                                 >  6:118093/1‑148 (MQ=255)                                                                                gTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGa                                                               <  8:106453/149‑1 (MQ=255)                                                                                      gAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCg                                                         <  8:200753/149‑1 (MQ=255)                                                                                              gA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCACACGGCCGCAGAGCGGCCCAGAGCCGTGTGAGTTGTGCTAAGCAAAGCTGGAAATGGCATGGCGGGTTTAGATCTAGCCAATAAAACGCATGGGGAGTGGCGGGTAGtta                                                 >  4:153914/1‑148 (MQ=255)                                                                                               a‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTACCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtct                                                >  3:232569/1‑149 (MQ=255)                                                                                                 cccAGGTCGAGATCTGAATGGGGTGATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtctc                                               <  7:111548/149‑1 (MQ=255)                                                                                                    aGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGt                                                      >  6:184151/1‑139 (MQ=255)                                                                                                    aGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGt                                                      <  5:184151/139‑1 (MQ=255)                                                                                                    aGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgtc                                            <  5:240091/149‑1 (MQ=255)                                                                                                       tCGAGATCTGAATGGGGTTATACAGAAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                         >  5:69495/1‑53 (MQ=18)                                                                                                       tCGAGATCTGAATGGGGTTATACAGAAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                         <  6:69495/53‑1 (MQ=18)                                                                                                          agaTCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGAGTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAgtg                                      <  7:262355/149‑1 (MQ=255)                                                                                                           gaTCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCCATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAgtgg                                     >  1:157803/1‑149 (MQ=255)                                                                                                                 aaTGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGc                                  <  7:238137/146‑1 (MQ=255)                                                                                                                     gggTTATACAGCAGGCTCAGAGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCa                            <  3:297194/148‑1 (MQ=255)                                                                                                                      ggTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCAGCCAgg                          >  8:154015/1‑149 (MQ=255)                                                                                                                      ggTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGGTTCGATCTCGCCCATCAACCGCGTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCACGCCAgg                          >  4:25060/1‑149 (MQ=255)                                                                                                                        ttATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCa                        <  2:124294/149‑1 (MQ=255)                                                                                                                            aCAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCGGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGttt                    >  7:160502/1‑149 (MQ=255)                                                                                                                             gacaagGCTCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGttt                    >  5:143816/5‑148 (MQ=255)                                                                                                                              agcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGt                                                             <  5:97820/106‑1 (MQ=255)                                                                                                                              agcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGt                                                             >  6:97820/1‑106 (MQ=255)                                                                                                                                cagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTa                                                     <  6:227987/112‑1 (MQ=255)                                                                                                                                aagGCTCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAgg                >  2:273400/2‑149 (MQ=255)                                                                                                                                     tCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCGTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCt            >  4:285317/1‑148 (MQ=255)                                                                                                                                     tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAggtggtg           >  1:109752/1‑144 (MQ=255)                                                                                                                                     tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCt            >  7:283291/1‑148 (MQ=255)                                                                                                                                     tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTc           >  3:7382/1‑149 (MQ=255)                                                                                                                                          gtcGCAGAGAGGGCCTGTTCCTTGTCAGTTCTGCTCTGCTTGGGTGGAGTTGCGTTGGCGGGTTTCGATCTCGGGAAGCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCtc      <  3:175728/147‑2 (MQ=255)                                                                                                                                              cAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTGc  <  3:153914/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    CTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGA‑CCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTGC  >  NC_002947/1376180‑1376468

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑