Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,070,248	+GA	intergenic (+48/+94)	aroP‑I → / ← PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,070,246	1	.	G	89.9%	239.1 / 20.7	80	intergenic (+46/+96)	aroP‑I/PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (5/3);  new base G (33/38);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_002947	1,070,246	2	.	A	89.7%	235.6 / 21.3	80	intergenic (+46/+96)	aroP‑I/PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (5/3);  new base A (32/38);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAAC  >  NC_002947/1070104‑1070384                                                                                                                                                     ||                                                                                                                                           gccgccGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑tgt                                                                                                                                              >  6:410746/1‑143 (MQ=255) gccgccGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑Gagaac                                                                                                                                      >  4:25529/1‑145 (MQ=255)  ccgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑Gagatcg                                                                                                                                     <  2:209867/149‑6 (MQ=255)    gccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑Gagatcgag                                                                                                                                   <  5:12386/149‑8 (MQ=255)    gccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTGGTTCGGAAGGGGTATGACGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑Gagatcgag                                                                                                                                   >  6:456745/1‑142 (MQ=255)            gCATCTCGGTAATGCTGATCCCGGCATGGATCGACTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGATGAACGGGTTTTTTCGTAAGGGTTATGCCGTTAGGTTTTCGGCGCATT‑‑TGA‑‑GATCGAGCGCcgcg                                                                                                                              >  1:54497/1‑149 (MQ=255)            gCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCcgcgc                                                                                                                             >  6:426922/1‑149 (MQ=255)            gCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCcgcgc                                                                                                                             >  5:123624/1‑149 (MQ=255)            gCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCcgcgc                                                                                                                             >  6:255857/1‑149 (MQ=255)            gCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCcgcgc                                                                                                                             >  5:63071/1‑149 (MQ=255)            gCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGCTTTCATCGC‑‑C‑‑TGTGAGATCGAGcgccgccc                                                                                                                             >  2:286469/1‑147 (MQ=38)               tctcGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑Gagatcga                                                                                                                                    <  7:120383/137‑7 (MQ=255)               tctcGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑Gagatcga                                                                                                                                    >  8:120383/1‑131 (MQ=255)                  cGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGGGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGGTGGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGCGATCGAGCGCCGCgcgggcgg                                                                                                                       >  5:380624/1‑149 (MQ=255)                    gTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGgcggcg                                                                                                                     <  5:375405/149‑1 (MQ=255)                       attcTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                  <  2:389315/146‑1 (MQ=38)                       aTGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                  <  4:57326/149‑1 (MQ=39)                        tGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCg                                                                                                                 <  5:271455/149‑1 (MQ=255)                         gCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CATTGA‑‑GATCGAGCGCCGCGCGGgcggcgggcgg                                                                                                                >  2:326004/1‑144 (MQ=255)                         gCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                >  8:72546/1‑149 (MQ=255)                           tGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                              >  3:355887/1‑149 (MQ=255)                               cccGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATAGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                           >  1:243763/1‑148 (MQ=255)                                 cGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAACGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGCGGTAGGCCTTTAGGCTTTCATCGC‑‑C‑‑TGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAg                                                                                                        <  1:286469/149‑1 (MQ=25)                                        gATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATg                                                                                                  <  7:262214/148‑1 (MQ=255)                                         aTCG‑CTGAGCTGTGGGTGGCTTTCAAGCTGAAGACCGCTCGCCAGGCCAAGGGCTGATCGGGTATGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaa                                                                                                <  5:456745/149‑1 (MQ=38)                                          tCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                               <  7:72546/149‑1 (MQ=255)                                            g‑ctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                             <  3:25529/149‑1 (MQ=255)                                              ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGcc                                                                                            >  5:181965/1‑149 (MQ=255)                                                 tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGCTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaa                                                                                         <  5:410746/149‑1 (MQ=255)                                                  gctgtgGGTGGCTTTCAAGCTGAAGACGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                                        <  3:139657/149‑1 (MQ=255)                                                   ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                                        <  5:255857/148‑1 (MQ=255)                                                      tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                           >  7:168405/1‑126 (MQ=255)                                                      tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                           <  8:168405/126‑1 (MQ=255)                                                      tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                                                     <  3:116379/132‑1 (MQ=255)                                                      tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                                                     >  4:116379/1‑132 (MQ=255)                                                            gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGCATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                                    <  3:386893/127‑1 (MQ=37)                                                            gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGCATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                                    >  4:386893/1‑127 (MQ=37)                                                              tttCAAGCTGAAGAAGGCCCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                            >  4:399990/1‑149 (MQ=255)                                                                tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTgg                                                                           <  4:355887/148‑1 (MQ=255)                                                                tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTgg                                                                           <  1:326004/148‑1 (MQ=255)                                                                   aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                                       >  5:88353/1‑149 (MQ=255)                                                                    gCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGcc                                                                      <  5:426922/149‑1 (MQ=255)                                                                    gCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGcc                                                                      <  2:243763/149‑1 (MQ=255)                                                                       gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCCCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                                                    >  8:69407/1‑148 (MQ=255)                                                                           aagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgc                                                               >  1:397821/1‑149 (MQ=255)                                                                                  gCCAGGCCCAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGAGCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                          <  6:380624/147‑1 (MQ=38)                                                                                     aggcacaGGGCTGAGCGGGTGTGTTCGTACGGGGTAGGCCTTTCGGTTTGAGGCGC‑‑C‑‑TGTTCGCTCGAGCGCCGCGCGGGCGGCGCTCGATCGCACAGGCGATGAAAGCCTAAAGGCATGCGGCTGGCGGCCCTGACGCGTTCCCCTTc                                                     <  7:87453/143‑1 (MQ=25)                                                                                       gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCa                                                   >  6:407075/1‑149 (MQ=255)                                                                                        ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                                  >  6:224547/1‑149 (MQ=255)                                                                                          aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                                                >  3:152808/1‑149 (MQ=255)                                                                                             ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGCGCGCCCCGCGGGCGGCGCGCGGTCTGACAGGCGATGAAAGCCTAAAGGCCAGCGGGGTGCGGCCCGGCCGCGGTCCCCGTGCATGcccc                                             >  6:292381/1‑149 (MQ=255)                                                                                             ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                                              <  5:292381/148‑1 (MQ=255)                                                                                               cTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCt                                            >  7:75034/1‑148 (MQ=255)                                                                                                  aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAt                                        >  6:408698/1‑149 (MQ=255)                                                                                                   tCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                                       >  2:358483/1‑149 (MQ=255)                                                                                                    cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTAGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGCCGCGTTCCCCTTCCATGCCCCTCCAtgc                                      >  3:401457/1‑149 (MQ=38)                                                                                                     gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                     <  6:63071/149‑1 (MQ=255)                                                                                                     gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                     <  2:54497/149‑1 (MQ=255)                                                                                                     gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                     >  5:328478/1‑149 (MQ=255)                                                                                                        tttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGCCGCGTTCCCCTTCCATGCCCCTCAATGCTGcc                                  >  1:278691/1‑149 (MQ=255)                                                                                                           gTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAg                                                                                       <  7:346066/93‑1 (MQ=25)                                                                                                           gTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAg                                                                                       >  8:346066/1‑93 (MQ=25)                                                                                                           gTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc                                <  6:123624/148‑1 (MQ=255)                                                                                                              cGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGc                             <  4:152808/148‑1 (MQ=255)                                                                                                               gTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                            <  7:69407/148‑1 (MQ=255)                                                                                                               gTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATAGAGCGCCGCGCGGGCGGCGCGCGATCTGACAGGCGATGAAAGCCTAAAGGCATGCGTCTTGCGGCCCGGACGCGTTACCCATACATTCCCCTCAATGCAGCCTcacgc                            >  4:12699/1‑144 (MQ=25)                                                                                                                  aGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc                        <  6:88353/149‑1 (MQ=255)                                                                                                                    gggTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac                      <  1:358483/149‑1 (MQ=255)                                                                                                                        aTGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc                                   <  5:8090/132‑1 (MQ=255)                                                                                                                        aTGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc                                   >  6:8090/1‑132 (MQ=255)                                                                                                                        aTGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc                                   >  6:1815/1‑132 (MQ=255)                                                                                                                        aTGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgctgc                                   <  5:1815/132‑1 (MQ=255)                                                                                                                            cGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg              <  8:75034/149‑1 (MQ=255)                                                                                                                            cGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg              >  1:244019/1‑149 (MQ=255)                                                                                                                             gTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg             <  7:325616/149‑1 (MQ=255)                                                                                                                                aGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggtg           <  1:403185/148‑1 (MQ=255)                                                                                                                                  gTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg              <  6:252762/143‑1 (MQ=255)                                                                                                                                  gTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg              >  5:252762/1‑143 (MQ=255)                                                                                                                                   ttttCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa       <  3:399990/149‑1 (MQ=255)                                                                                                                                      tCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGGCTCCGCCAAGCACGCACAGTGGTGTACAGc     <  1:39137/148‑1 (MQ=255)                                                                                                                                        ggCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAac  <  5:408698/149‑1 (MQ=255)                                                                                                                                                     ||                                                                                                                                           GCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCG‑CTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAAC  >  NC_002947/1070104‑1070384

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