Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,295,013	Δ1 bp	9.2%	intergenic (+29/+41)	slyB → / ← nhaA‑I	outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,295,013	0	A	.	9.2%	199.8 / 14.1	65	intergenic (+29/+41)	slyB/nhaA‑I	outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base . (4/2);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGATGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGC  >  NC_002947/1294870‑1295155                                                                                                                                                |                                                                                                                                                aGATGATGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACAGTAGTGCAGTAATAAAAAACCCCGACATGgct                                                                                                                                            >  6:244978/1‑149 (MQ=255)     gatgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGctgctg                                                                                                                                        >  3:325003/1‑149 (MQ=255)      atgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCtgctgc                                                                                                                                       <  5:292372/149‑1 (MQ=255)              aTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCggg                                                                                                                              <  6:245248/149‑1 (MQ=255)               tGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGcc                                                                                                                                      <  7:185647/140‑1 (MQ=255)               tGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGcc                                                                                                                                      >  8:185647/1‑140 (MQ=255)                  gcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGttt                                                                                                                          >  8:71980/1‑149 (MQ=39)                     ccTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGtttttt                                                                                                                       <  7:71980/149‑1 (MQ=39)                      cTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGttgttt                                                                                                                       >  2:246215/1‑148 (MQ=255)                            gCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGAGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGAGCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCt                                                                                                                >  8:169135/1‑149 (MQ=255)                            gCAGGCCGTGCAGAAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCt                                                                                                                >  1:200017/1‑149 (MQ=255)                             cAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGTCTGTAGACGGTACCAGTCGCGTTTCGCACTGAACGTAGTGCAGTAATAAAAAACCCCGACATGTCTTCTGCCA‑GTCGGGGTTTTTTCATGCCTc                                                                                                               <  8:317105/149‑1 (MQ=255)                              aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCa                                                                                                              <  1:251860/149‑1 (MQ=255)                              aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCa                                                                                                              >  2:251860/1‑149 (MQ=255)                                gCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGc                                                                                                            >  2:49445/1‑149 (MQ=255)                                 ccGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGcc                                                                                                           <  5:244978/149‑1 (MQ=255)                                         gagaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAgg                                                                                                   >  1:265155/1‑149 (MQ=255)                                         gagaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAgg                                                                                                   <  8:305201/149‑1 (MQ=255)                                         gagaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTAGCACTGACCGTAGTGCAGTAATAAAAAACCCCCACATGGCTGCTGCCA‑GTCGGGGTTTTTTAATTCCTCAGCCAATTCAgg                                                                                                   >  2:317763/1‑149 (MQ=255)                                          agaACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGc                                                                                                  <  4:325003/149‑1 (MQ=255)                                               gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGttt                                                                                                                          <  1:134752/120‑1 (MQ=255)                                               gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGttt                                                                                                                          >  2:134752/1‑120 (MQ=255)                                                     ttCCGTTTTGGAGACCGACTGCCCATCATCACAGTATACGGTACAAGTCGCGATTCGCAATTACCGTAGTGCAGTAATAAAAAACCCCGAAATGAATGCTGACA‑GTCAGGGTTTTTTCATGCCTCAGCCAAGTCATGCAATGCTGGttt                                                                                       >  4:186770/1‑149 (MQ=255)                                                         gtgtTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCg                                                                                   >  4:144035/1‑149 (MQ=255)                                                          tgtTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCgg                                                                                  <  1:49445/149‑1 (MQ=255)                                                               gCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTgg                                                                              <  1:246215/148‑1 (MQ=255)                                                                     gAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGc                                                                       >  2:128291/1‑149 (MQ=255)                                                                        tGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGt                                                                    <  5:112381/149‑1 (MQ=255)                                                                             atcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACgg                                                               >  3:9355/1‑149 (MQ=255)                                                                              tcatGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGc                                                              >  2:183045/1‑149 (MQ=255)                                                                               catGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTg                                                                        >  4:25689/1‑138 (MQ=255)                                                                               catGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTg                                                                        <  3:25689/138‑1 (MQ=255)                                                                                  gACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTc                                                                                                                      >  8:21878/1‑89 (MQ=255)                                                                                  gACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTc                                                                                                                      <  7:21878/89‑1 (MQ=255)                                                                                  gACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAg                                                          >  6:229600/1‑149 (MQ=255)                                                                                   aCTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAg                                                          <  5:30/148‑1 (MQ=255)                                                                                       tAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCg                                                       <  2:265155/147‑1 (MQ=255)                                                                                         gACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCa                                                   >  7:130917/1‑149 (MQ=255)                                                                                          aCGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAc                                                  >  1:50512/1‑149 (MQ=255)                                                                                           cGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc                                                 <  2:50512/149‑1 (MQ=255)                                                                                                ccAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgcag                                             <  1:128291/148‑1 (MQ=255)                                                                                                 cAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgcagc                                            >  6:244248/1‑148 (MQ=255)                                                                                                       cgTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaac                                     >  7:300364/1‑149 (MQ=255)                                                                                                        gTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCaa                                                                                                         <  1:261279/80‑1 (MQ=255)                                                                                                        gTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCaa                                                                                                         >  2:261279/1‑80 (MQ=255)                                                                                                           tCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaa                                  <  5:244248/148‑1 (MQ=255)                                                                                                                   aCCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAg                          <  1:317763/148‑1 (MQ=255)                                                                                                                       tAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGc                     <  8:130917/149‑1 (MQ=255)                                                                                                                       tAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGc                     >  3:222001/1‑149 (MQ=255)                                                                                                                            cAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCAt                 <  3:144035/148‑1 (MQ=255)                                                                                                                             aGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGGTGGTTTTGCGGATGGCCATTGCGGTGACGGAGTATCCGATCACAGCAGCAAAGAACGAACCAGTAAGAATGCCCATGc               >  1:115902/1‑149 (MQ=255)                                                                                                                                  taaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                 >  8:110498/1‑126 (MQ=255)                                                                                                                                  taaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                 <  7:110498/126‑1 (MQ=255)                                                                                                                                  gaaaaaaCCCCGAC‑TGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTc           >  8:305233/2‑148 (MQ=38)                                                                                                                                   aaaaaaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCa         >  4:88924/1‑149 (MQ=255)                                                                                                                                    aaaaaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCAt        <  4:9355/149‑1 (MQ=255)                                                                                                                                     aaaaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTGTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATg       >  3:212316/1‑149 (MQ=255)                                                                                                                                       aaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAAt                       <  7:255708/131‑1 (MQ=255)                                                                                                                                       aaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAAt                       >  8:255708/1‑131 (MQ=255)                                                                                                                                         ccccGAC‑TGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCgg   >  5:79437/1‑149 (MQ=255)                                                                                                                                          cccGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCAt        <  5:210835/143‑1 (MQ=255)                                                                                                                                          cccGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCAt        >  6:210835/1‑143 (MQ=255)                                                                                                                                          cccGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc  >  8:311143/1‑149 (MQ=255)                                                                                                                                               cATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc  <  6:354217/144‑1 (MQ=255)                                                                                                                                               cATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGc  >  5:354217/1‑144 (MQ=255)                                                                                                                                                |                                                                                                                                                AGATGATGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGC  >  NC_002947/1294870‑1295155

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