Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,936,917	A→T	8.8%	intergenic (‑328/‑90)	flhA ← / → PP_4345	flagellar export pore protein/GntR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,936,917	0	A	T	8.8%	188.8 / 8.0	68	intergenic (‑328/‑90)	flhA/PP_4345	flagellar export pore protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base T (3/3);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

CCTCTACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCG  >  NC_002947/4936770‑4937061                                                                                                                                                     |                                                                                                                                                 ccTCTACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGaa                                                                                                                                                  >  2:285937/1‑148 (MQ=255)    ctACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTCGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAcat                                                                                                                                               >  8:32955/1‑148 (MQ=255)       cGGCAATGGCCAAAGGCTTAAGCTCTTTGAAAGTCAAAGATAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATg                                                                                                                                           <  7:277109/149‑1 (MQ=255)       cGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAcatcat                                                                                                                                            <  7:32955/148‑1 (MQ=255)        ggCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGc                                                                                                                                          <  6:258713/149‑1 (MQ=255)            aTGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGccttt                                                                                                                                      <  8:710/149‑4 (MQ=255)             ttgCCAAACGCATAATCCTTTTGAAATCCACACACATGGTTGCGGCAACAAAAAAAAAATCTTACCCAACTTATCATTCAAAACCTCACGGGTTCAAAGCCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAAc                                                                                                                                      <  1:285937/147‑1 (MQ=255)             tGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACt                                                                                                                                     <  1:209877/149‑1 (MQ=255)               gCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGcctttga                                                                                                                                    <  5:138390/148‑6 (MQ=255)                ccAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGcctttga                                                                                                                                    >  2:193356/1‑142 (MQ=255)                ccAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGcctttga                                                                                                                                    <  1:193356/147‑6 (MQ=255)                 cAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGc                                                                                                                                          <  3:301077/140‑1 (MQ=255)                 cAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTc                                                                                                                                 <  8:117516/149‑1 (MQ=255)                    aGGCTTAAGCTCTTGGAAACCCACAGACAGGGTCGCGCCACGACAAAAAAAATCTGAGCGACTTGATCAAGCAAAAGCTCACGCTTTCTATC‑CGGATCTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGAt                                                                                                                              <  8:217398/149‑1 (MQ=255)                     ggCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATg                                                                                                                             <  3:180241/149‑1 (MQ=255)                        ttAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGCGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATgg                                                                                                                            >  5:226600/1‑147 (MQ=255)                        ttAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATgg                                                                                                                            >  3:185486/1‑147 (MQ=255)                        ttAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTa                                                                                                                          <  4:185486/149‑1 (MQ=255)                            gCTCTTGGAAAGCCACAGAGAGTGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTGGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGctc                                                                                                                      <  8:154271/149‑1 (MQ=255)                            gCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGctc                                                                                                                      <  4:213649/149‑1 (MQ=255)                                ttGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCGAATGACTACTCAAAGGCATGAACATCATGCCAAGTGCTCGATGGTAGCTCTAGt                                                                                                                  >  5:192028/1‑149 (MQ=255)                                   gAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAg                                                                                                               >  4:265192/1‑149 (MQ=255)                                   gAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAg                                                                                                               >  1:73771/1‑149 (MQ=255)                                     aaGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTc                                                                                                             >  3:188314/1‑149 (MQ=255)                                       gCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCtg                                                                                                           <  6:226600/149‑1 (MQ=255)                                                  ggTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCg                                                                                                <  2:73771/149‑1 (MQ=255)                                                             gACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGt                                                                                                                  <  3:276754/120‑1 (MQ=255)                                                             gACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGt                                                                                                                  >  4:276754/1‑120 (MQ=255)                                                             gACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtttt                                                                                     <  5:142596/149‑1 (MQ=255)                                                             gACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtttt                                                                                     <  3:265192/149‑1 (MQ=255)                                                               aaaaaaaaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTATAGATCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTTAGTCTGTGATACCTTCGCCAGCGCTTTTGa                                                                                   <  5:374423/148‑1 (MQ=255)                                                                            aGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGa                                                                      <  8:289993/149‑1 (MQ=255)                                                                                aCTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCg                                                                  <  8:275525/149‑1 (MQ=255)                                                                                  ttGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGga                                                                <  5:137536/149‑1 (MQ=255)                                                                                    gATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGgag                                                               >  8:308968/1‑148 (MQ=255)                                                                                     aTCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtt                                                                                       >  7:75334/1‑123 (MQ=255)                                                                                     aTCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCtt                                                                                       <  8:75334/123‑1 (MQ=255)                                                                                        aTGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAt                                                          >  3:190592/1‑149 (MQ=255)                                                                                         tGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtg                                                         <  4:188314/149‑1 (MQ=255)                                                                                         tGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtg                                                         <  5:138548/149‑1 (MQ=255)                                                                                          gCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCgg                                                                 <  3:221626/140‑1 (MQ=255)                                                                                          gCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCgg                                                                 >  4:221626/1‑140 (MQ=255)                                                                                          gCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtgg                                                        >  6:279834/1‑149 (MQ=255)                                                                                            aaaaGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtggtg                                                      >  5:224191/1‑149 (MQ=255)                                                                                              acgCTCACGTCTTCTAGG‑CGGATGTAGAGACCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGGGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGc                                                    <  4:261555/147‑1 (MQ=255)                                                                                                    aCGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAg                                                                                                                   >  8:187728/1‑80 (MQ=255)                                                                                                    aCGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAg                                                                                                                   <  7:187728/80‑1 (MQ=255)                                                                                                    aCGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGa                                                                                   >  1:210692/1‑112 (MQ=255)                                                                                                    aCGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGa                                                                                   <  2:210692/112‑1 (MQ=255)                                                                                                    aCGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAActc                                              >  5:334534/1‑149 (MQ=255)                                                                                                    aCGGATTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAAct                                               <  8:309070/148‑1 (MQ=255)                                                                                                      ggTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAActctc                                            <  7:308968/149‑1 (MQ=255)                                                                                                       gTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAActctct                                           <  6:224191/149‑1 (MQ=255)                                                                                                            tATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGa                                        >  5:265818/1‑147 (MQ=255)                                                                                                               g‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCt                                   >  4:332278/1‑149 (MQ=255)                                                                                                                  ggATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCtt                                  >  5:372341/1‑148 (MQ=255)                                                                                                                   gATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTa                                >  6:86089/1‑149 (MQ=255)                                                                                                                        agagagCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGtt                           >  2:331416/1‑149 (MQ=255)                                                                                                                             gCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGa                      >  3:87956/1‑149 (MQ=255)                                                                                                                              ccAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGAc                     <  5:279834/149‑1 (MQ=255)                                                                                                                               cAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATAc                                                                                                     <  8:133744/68‑1 (MQ=255)                                                                                                                               cAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATAc                                                                                                     >  7:133744/1‑68 (MQ=255)                                                                                                                                         tCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCa          >  6:343178/1‑149 (MQ=255)                                                                                                                                         tCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCa          >  5:381633/1‑149 (MQ=255)                                                                                                                                         tCAAAGGCATGAAAATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTATCTCGATAGCTTTATCGTTCCTGACAAACATCGCCa          >  4:344887/1‑149 (MQ=255)                                                                                                                                           aaaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCAt         >  2:374590/1‑148 (MQ=255)                                                                                                                                            ttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATa        >  7:6085/3‑148 (MQ=255)                                                                                                                                              ggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATaaaa     >  3:340639/1‑149 (MQ=255)                                                                                                                                                 aTGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCAATTGACTACATCGCTAGACACGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTACTGACAAACATAGCCATAAAAcag  >  1:311414/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CCTCTACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATG‑CGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCG  >  NC_002947/4936770‑4937061

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