Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,964,846	A→G	23.8%	L171P (CTG→CCG)	fleQ ←	transcriptional regulator FleQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,964,846	0	A	G	23.8%	110.1 / 31.8	63	L171P (CTG→CCG)	fleQ	transcriptional regulator FleQ
Reads supporting (aligned to +/- strand):  ref base A (26/22);  new base G (8/7);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGAC  >  NC_002947/4964702‑4964991                                                                                                                                                  |                                                                                                                                                  tcGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCAc                                                                                                                                               >  6:6333/1‑149 (MQ=255)    tGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCt                                                                                                                                            <  1:189632/149‑1 (MQ=255)      gCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTg                                                                                                                                           <  5:6333/148‑1 (MQ=255)         gAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGt                                                                                                                                       >  5:49671/1‑149 (MQ=255)               ttCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCgg                                                                                                                                    >  5:27410/1‑146 (MQ=255)               ttCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTAt                                                                                                                                 >  3:136649/1‑149 (MQ=255)               ttCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGAAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTAt                                                                                                                                 >  4:31534/1‑149 (MQ=255)               ttCCCTTTCCAGCCGCTCGGCGGGGACCGCGCCACAGTTCCCCGCCACTAACGGCGCTTCGCGGCGCTTGTAAGGGTCGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCCCCCAGGAGCGCCACGCGGGCGTCgg                                                                                                                                    <  6:27410/146‑1 (MQ=255)                        cagcagCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACtt                                                                                                                        <  3:233510/149‑1 (MQ=255)                        cagcagCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACtt                                                                                                                        <  6:82926/149‑1 (MQ=255)                             gcTCGGCCGGGATCGCGCCACAGTTGACCGGCTCGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCTCGCTACCAATTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTtgctgc                                                                                                                   <  2:335032/149‑1 (MQ=255)                             gcTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTtgctg                                                                                                                    >  1:64174/1‑148 (MQ=255)                                cGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGcatc                                                                                                                <  3:12025/149‑1 (MQ=255)                                  gCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGcatcat                                                                                                              >  3:136998/1‑149 (MQ=255)                                    cGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCt                                                                                                            >  3:231204/1‑149 (MQ=255)                                          cgcgCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCCGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCa                                                                                                      >  4:169970/1‑149 (MQ=255)                                           gcgcCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCAc                                                                                                     <  8:15335/149‑1 (MQ=255)                                                    ttGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAAt                                                                                            >  7:189994/1‑149 (MQ=255)                                                      gACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTTAATg                                                                                           <  8:189994/148‑1 (MQ=255)                                                         cGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCAc                                                                                       >  8:20683/1‑149 (MQ=255)                                                           gCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGa                                                                                     >  6:236613/1‑149 (MQ=255)                                                           gCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGAGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGa                                                                                     >  2:289600/1‑149 (MQ=255)                                                            cACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCTGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGAc                                                                                    <  5:236613/149‑1 (MQ=255)                                                            cACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGAc                                                                                    <  6:49671/149‑1 (MQ=255)                                                              cGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTg                                                                                  >  8:29610/1‑149 (MQ=255)                                                                aaCGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGt                                                                                <  7:29610/149‑1 (MQ=255)                                                                aaCGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGt                                                                                <  6:308341/149‑1 (MQ=255)                                                                aaCGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGCGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGt                                                                                <  1:289600/149‑1 (MQ=255)                                                                      gcTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAAc                                                                          >  1:37407/1‑149 (MQ=255)                                                                        ttCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAAc                                                                          >  5:180725/1‑147 (MQ=255)                                                                          cgcgGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCCGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACCAGTCGTGCCattcag                                                                       <  3:169970/148‑6 (MQ=255)                                                                           gcggcgCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGc                                                                     >  2:255329/1‑149 (MQ=255)                                                                                gcTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCgg                                                                <  2:64174/149‑1 (MQ=255)                                                                                  ttGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGGTACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCgg                                                                >  1:250867/1‑148 (MQ=255)                                                                                      aaTGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATgg                                                                                          <  3:47620/117‑1 (MQ=255)                                                                                      aaTGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATgg                                                                                          >  4:47620/1‑117 (MQ=255)                                                                                            aGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAAt                                                                                            >  3:290604/1‑109 (MQ=255)                                                                                            aGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAAt                                                                                            <  4:290604/109‑1 (MQ=255)                                                                                                 aGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTcgcg                                               >  6:215071/1‑149 (MQ=255)                                                                                                    ttGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTg                                            <  3:31534/149‑1 (MQ=255)                                                                                                          gcTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGAcc                                      >  1:309142/1‑149 (MQ=255)                                                                                                              ccACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGc                                  >  3:300855/1‑149 (MQ=255)                                                                                                                aCTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAAc                                                             <  5:102247/120‑1 (MQ=255)                                                                                                                aCTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAAc                                                             >  6:102247/1‑120 (MQ=255)                                                                                                                aCTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTc                                >  7:116901/1‑149 (MQ=255)                                                                                                                    ccTTGCCGGTGCCGG‑ACTCACCCAGGATCGGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTcgcg                             <  4:231204/148‑1 (MQ=255)                                                                                                                     cTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTcgcgcg                           >  4:251063/1‑149 (MQ=255)                                                                                                                        gccggtgccgg‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGt                                                                                >  4:8387/1‑93 (MQ=255)                                                                                                                        gccggtgccgg‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGt                                                                                <  3:8387/93‑1 (MQ=255)                                                                                                                        gccggtgccgg‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGcc                         <  2:250867/148‑1 (MQ=255)                                                                                                                         ccggtgccgg‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTg                       >  3:109105/1‑149 (MQ=255)                                                                                                                            gtgccgg‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTc                    <  2:187603/149‑1 (MQ=255)                                                                                                                                 gg‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACa               >  4:297454/1‑149 (MQ=255)                                                                                                                                    aCTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACAtt             <  4:136998/149‑1 (MQ=255)                                                                                                                                        aCCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTcgc         <  2:37407/149‑1 (MQ=255)                                                                                                                                         cccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTcgcg        <  4:109105/149‑1 (MQ=255)                                                                                                                                          ccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACAttt            >  2:145989/1‑144 (MQ=255)                                                                                                                                          ccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACAttt            <  1:145989/144‑1 (MQ=255)                                                                                                                                          ccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCgg       >  4:51918/1‑149 (MQ=255)                                                                                                                                          ccAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCgg       >  1:265195/1‑149 (MQ=255)                                                                                                                                           cAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCgg       >  3:281597/1‑148 (MQ=255)                                                                                                                                             ggATCAGAACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGATCAGATTGGGTTCGCCCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCgg       >  6:271399/1‑146 (MQ=255)                                                                                                                                               aTCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGAc  <  3:73284/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TCGTGGCCGAACAGTTCGCTTTCCAGCAGCTCGGCCGGGATCGCGCCACAGTTGACCGGCACGAACGGCGCTTCGCGGCGCTTGGAATGGTAGTGCAGGTTGCGCGCTACCACTTCCTTGCCGGTGCCGG‑ACTCACCCAGGATCAGCACGCTGGCGTCGGTATCCGCCACTTGCTGCATCATCTGCCGCACGTGCTGAATGGCACGACTGGTGCCAACCAGGCTGCGGAACAGATTGGGTTCGCGCTGGCGACCCCGCTCGCGCGCCTGGTCATACATTTCGCGGTAGAC  >  NC_002947/4964702‑4964991

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