Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,964,613	C→T	7.8%	D249N (GAC→AAC)	fleQ ←	transcriptional regulator FleQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,964,613	0	C	T	7.8%	185.0 / 4.1	64	D249N (GAC→AAC)	fleQ	transcriptional regulator FleQ
Reads supporting (aligned to +/- strand):  ref base C (27/32);  new base T (3/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.46e-01

GTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGAC  >  NC_002947/4964465‑4964757                                                                                                                                                     |                                                                                                                                                 gTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTc                                                                                                                                                  <  2:20446/149‑1 (MQ=255)   gCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGc                                                                                                                                                >  8:58373/1‑149 (MQ=255)    ccGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGcc                                                                                                                                               >  1:158862/1‑149 (MQ=255)           cGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTtcg                                                                                                                                        >  5:264274/1‑149 (MQ=255)            gATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTtc                                                                                                                                         <  8:179156/147‑1 (MQ=255)             aTCATGGTGTCGAGGTTCTTGTGCGGCGCAGCGATGATGCGCACGTCGATGCTTTGGGGCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTtcgtc                                                                                                                                      <  6:250005/149‑1 (MQ=255)                aTGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTAACCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGa                                                                                                                                    <  8:185466/148‑1 (MQ=255)                      tCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACa                                                                                                                              <  2:67654/148‑1 (MQ=255)                        gAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAggg                                                                                                                           <  3:120732/149‑1 (MQ=255)                        gAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAgg                                                                                                                            <  1:22217/148‑1 (MQ=255)                        gAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAgg                                                                                                                            <  8:206012/148‑1 (MQ=255)                           gTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTAc                                                                                                                        >  7:61873/1‑149 (MQ=255)                                 gtgCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTAACCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGt                                                                                                                  >  2:119612/1‑149 (MQ=255)                                  tgCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGtt                                                                                                                 <  6:305804/149‑1 (MQ=255)                                    cGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTgg                                                                                                               <  1:305346/149‑1 (MQ=255)                                    cGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTg                                                                                                                <  6:80202/148‑1 (MQ=255)                                     gTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGc                                                                                                              <  7:58373/149‑1 (MQ=255)                                         cAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGt                                                                                                          >  8:69984/1‑149 (MQ=255)                                          aGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGtt                                                                                                         <  3:130520/149‑1 (MQ=255)                                              atgatgCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGaa                                                                                                     <  5:50896/149‑1 (MQ=255)                                                 atgCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAAc                                                                                                   <  3:316294/148‑1 (MQ=255)                                                      cACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTTGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCgg                                                                                             <  4:10044/149‑1 (MQ=255)                                                        cGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCa                                                                                           <  2:158862/149‑1 (MQ=255)                                                           cGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACgg                                                                                        >  3:180096/1‑149 (MQ=255)                                                            gATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTTGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGc                                                                                       >  7:311767/1‑149 (MQ=255)                                                                      gTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGc                                                                             <  2:33474/149‑1 (MQ=255)                                                                           gttgCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCATACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCg                                                                         <  1:167352/148‑1 (MQ=255)                                                                           gttgCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGccc                                                                           <  4:180096/146‑1 (MQ=255)                                                                           gttgCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCgg                                                                        >  7:86477/1‑149 (MQ=255)                                                                             tgCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCAAGGCTGGTGATGGCCCCGGt                                                                       >  8:310276/1‑148 (MQ=255)                                                                                    cACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTAGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGgcgc                                                               >  2:152789/1‑149 (MQ=255)                                                                                     aCGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTTGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGcc                                                              <  8:311767/149‑1 (MQ=255)                                                                                        cgcTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCtt                                                           <  8:61873/149‑1 (MQ=255)                                                                                         gcTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTtc                                                          >  8:86519/1‑149 (MQ=255)                                                                                             gAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGa                                                                                 <  5:300163/122‑1 (MQ=255)                                                                                             gAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGa                                                                                 >  6:300163/1‑122 (MQ=255)                                                                                             gAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGt                                                      >  4:301081/1‑149 (MQ=255)                                                                                               aGGTACGCTCCTGAAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTgg                                                    >  5:199567/1‑149 (MQ=255)                                                                                                ggTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTGGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGc                                                   >  6:206945/1‑149 (MQ=255)                                                                                                 gTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTTGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGcc                                                  >  4:226734/1‑149 (MQ=255)                                                                                                    cGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGaa                                               <  1:152789/149‑1 (MQ=255)                                                                                                       tCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAg                                            >  2:276081/1‑149 (MQ=255)                                                                                                         cTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCgg                                                                                             <  6:301139/98‑1 (MQ=255)                                                                                                            cagcaCACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCGCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGc                                       >  7:146154/1‑149 (MQ=255)                                                                                                            cagcaCACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTc                                         <  1:119612/147‑1 (MQ=255)                                                                                                            cagcaCACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCg                                        <  7:86519/148‑1 (MQ=255)                                                                                                               cacaCGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCttt                                    <  6:264274/149‑1 (MQ=255)                                                                                                               cacaCGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGCCACGGCTGGTGATGGCCCCGGTAAAGGCGCACTTCTCGTGGCCGAACAGTTCGCtt                                     >  4:282335/1‑148 (MQ=255)                                                                                                                  acGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCca                                 <  3:301081/149‑1 (MQ=255)                                                                                                                     cagcagCTTGACCTGCATCGGCAACGGCATGTCGCCGATTCCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCtt                                     >  8:84361/1‑142 (MQ=255)                                                                                                                     cagcagCTTGACCTGCATCGGCAACGGCATGTCGCCGATTCCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCtt                                     <  7:84361/142‑1 (MQ=255)                                                                                                                      agcagcTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCcagcag                             >  4:176441/1‑149 (MQ=255)                                                                                                                              gACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCg                     >  2:283705/1‑149 (MQ=255)                                                                                                                               aCCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCgg                    >  6:146559/1‑149 (MQ=255)                                                                                                                                ccTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCggg                   <  7:69984/149‑1 (MQ=255)                                                                                                                                   gCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATc                >  5:5567/1‑149 (MQ=255)                                                                                                                                   gCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATc                >  4:28140/1‑149 (MQ=255)                                                                                                                                    cATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATcg               >  6:47112/1‑149 (MQ=255)                                                                                                                                      tCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATcgcg             <  3:176441/149‑1 (MQ=255)                                                                                                                                       cGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCgcgc            >  4:94602/1‑149 (MQ=255)                                                                                                                                         gCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCca          >  2:148290/1‑149 (MQ=255)                                                                                                                                           aaCGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCcaca        <  3:94602/149‑1 (MQ=255)                                                                                                                                              ggCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCTAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGtt     <  6:199567/149‑1 (MQ=255)                                                                                                                                                 aTGTTGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCACCAGCTCGGCCGGGATCGCGCCACAGTTGAc  >  4:180254/1‑149 (MQ=255)                                                                                                                                                 aTGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGAc  <  5:146559/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGCAGCACACGCAGCAGCTTGACCTGCATCGGCAACGGCATGTCGCCGATTTCGTCGAGGAACAGGGTACCGCCGTTGGCCAGTTCGAAACGCCCGGCACGGCTGGTGATGGCCCCGGTAAAGGCGCCCTTCTCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGAC  >  NC_002947/4964465‑4964757

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