Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	340,051	A→C	7.0%	intergenic (+308/‑114)	yahG → / → yahI	DUF1116 family protein/carbamate kinase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	340,051	0	A	C	7.0%	198.7 / 4.8	68	intergenic (+308/‑114)	yahG/yahI	DUF1116 family protein/carbamate kinase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (27/36);  new base C (2/3);  total (29/39)
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

CTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTG  >  NC_000913/339907‑340195                                                                                                                                                 |                                                                                                                                                 cTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGc                                                                                                                                              >  3:14979/1‑149 (MQ=255)  tGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGcc                                                                                                                                             >  1:88803/1‑149 (MQ=255)  tGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGcc                                                                                                                                             >  5:204886/1‑149 (MQ=255)  tGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGAATTTGTCACCATTGCCGGATGCGGCGTGAACGcc                                                                                                                                             >  7:156737/1‑149 (MQ=255)    attaGCGTATCGGATCAGGCATTTGCCCCCACTGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCtt                                                                                                                                           <  8:505/146‑1 (MQ=14)       aGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATAAGGCGTTCACGccgcatc                                                                                                                                        <  8:97534/149‑6 (MQ=14)          gTAGCGCATCAGGCATTAGCCACAACTGCCTTATCCCAACAGCACCCCTTAACCGGACTACAACTAGCTCGAGAAACTGTAGGCCTGATAAGCGTATCGAATCAGGAATTTATCAACATTGCCGGATGCTGCCAGAAATCATTATCCgg                                                                                                                                     >  4:46763/1‑149 (MQ=9)             gcgcATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                   <  4:163333/148‑1 (MQ=255)             gcgcATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                   <  8:151524/148‑1 (MQ=255)                 aTCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                    <  6:181072/143‑1 (MQ=255)                 aTCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                    >  5:181072/1‑143 (MQ=255)                 aTCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGa                                                                                                                              <  1:247004/149‑1 (MQ=255)                 aTCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGa                                                                                                                              <  1:140385/149‑1 (MQ=255)                              caccacTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAt                                                                                                                 >  6:300669/1‑149 (MQ=255)                               accacTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAtt                                                                                                                >  5:14699/1‑149 (MQ=255)                                 cacTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTg                                                                                                               >  4:130975/1‑148 (MQ=255)                                  acTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTa                                                                                                             <  4:44179/149‑1 (MQ=255)                                   cTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAg                                                                                                            <  7:259439/149‑1 (MQ=255)                                      ccGTATGCGGCGTGTACGCCTTATCCGGCCTACGAGTGTCGCTAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGc                                                                                                          <  5:310279/148‑1 (MQ=255)                                      ccGGATGCGGCGTTTACGCCTTATCCGGCCTACGAGTGGCGCTAGAATCTTTTGGCCTGTTAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGc                                                                                                          <  4:169131/148‑1 (MQ=14)                                         gATGCGTCGTGTACGCCTTCTCCGGCCTACGACTGGCGCTAGAATCTGTAGGCCTGCTAAGCGTAGCGCATCAGGCATTTGTCTCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGa                                                                                                      <  4:60912/149‑1 (MQ=11)                                             cggcgTGGACGCCCTAGCCGGCCTACGATGGGCGCGAGACTCTGTCGGCCTGATAAGCGGAGCGCATCAGGCATTTGTCCCCATTGCCGCCGGCGGCGTGAACGCCTTATCCGGCCTACAATTCTCGCGCCACTCGTAGGCCGGATaa                                                                                                   <  6:244884/148‑1 (MQ=2)                                                  gggACGCCTTAGCCGGCCTCCGAGTGGTGTTTGTATCTGTTGGCCGGCTCAGCGTAGCGCTTCAGGCTTTTGTCACCATTGCTGGATGCGGCGTGACCGCCTTATCGGGCCTAGGAGGGGGGCGTTATTTGTAGGCCTGATAAGCGTAg                                                                                             <  2:102440/148‑1 (MQ=9)                                                          ttATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGcagcag                                                                                     >  3:173176/1‑149 (MQ=255)                                                          ttATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                     >  6:120356/1‑149 (MQ=255)                                                              ccGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAtt                                                                                                                >  6:202303/1‑118 (MQ=255)                                                              ccGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAtt                                                                                                                <  5:202303/118‑1 (MQ=255)                                                                ggCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTc                                                                               >  2:296398/1‑149 (MQ=255)                                                                 gCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGCTGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCt                                                                              <  5:300234/149‑1 (MQ=255)                                                                  ccTACGAGTGGCGCGAGCATCTGTAGGCCGGATACGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAGGCGTAGCGCATCAGGCAGTCTg                                                                             <  3:4312/149‑1 (MQ=255)                                                                        aGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTg                                                                       <  3:130975/149‑1 (MQ=255)                                                                        aGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTg                                                                       <  6:14699/149‑1 (MQ=255)                                                                           ggCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTc                                                                    <  3:125954/149‑1 (MQ=255)                                                                            tcgctAGAATCTGTAGGCCTGATAAGAGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGATCGCCTTATCCGGCCTACGAGTGGCGCGATAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCa                                                                   <  2:117397/148‑1 (MQ=255)                                                                            gcgcttgaATCTGTTGGCCTGATAAGCGTAGCGCATCAGGCATTTTTCACCATTTCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGGGGGGCTTTATTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCa                                                                   <  8:326098/143‑1 (MQ=17)                                                                               ctttaaTCTGTATGCCTGATAAGCGTATCGCATCAGTAATTTTTCACCATTGCCGTATGCGTCGTGAACGCCTTATCCGGCATACTATTGGCGCGATAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATaa                                                                <  2:129105/145‑1 (MQ=17)                                                                                  gaATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAccc                                                             >  1:74982/1‑149 (MQ=255)                                                                                    aTCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCa                                                           <  2:88803/149‑1 (MQ=255)                                                                                       tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCa                                                        >  2:35665/1‑149 (MQ=255)                                                                                           ggCCTGATAAGCGTAGCGCATCAGGCTTTTGTCACCATTGCCGGATGCGGCGTGAACGCATTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCct                                                    <  5:171061/149‑1 (MQ=255)                                                                                           ggCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAc                                                               <  2:230126/138‑1 (MQ=255)                                                                                           ggCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAc                                                               >  1:230126/1‑138 (MQ=255)                                                                                                 aTAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAg                                              >  2:297760/1‑149 (MQ=255)                                                                                                     gCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCGCCCTCTGTAGCGGa                                          >  2:182299/1‑149 (MQ=255)                                                                                                        tAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACAt                                       <  6:142663/149‑1 (MQ=255)                                                                                                          gcgcATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaa                                     <  1:296398/149‑1 (MQ=255)                                                                                                             cATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaacaa                                  >  4:153826/1‑149 (MQ=255)                                                                                                              aTCAGGCATTTGTCTCCATTGCCGTATGCGGCGTTAAGTCCTTATCTGGCCTAAGTGTGGCGCGATAATTTTAGCCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaacaac                                 <  3:46763/149‑1 (MQ=25)                                                                                                              aTCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaacaa                                  >  7:34420/1‑148 (MQ=255)                                                                                                              aTCAGGCATTTGTCAACATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaacaa                                  >  7:45001/1‑148 (MQ=255)                                                                                                                 aGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACAt                               >  3:129373/1‑148 (MQ=255)                                                                                                                  ggCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATg                              >  7:166112/1‑148 (MQ=255)                                                                                                                   gCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGa                             >  4:301525/1‑148 (MQ=255)                                                                                                                    cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGaaa                           <  2:74982/149‑1 (MQ=255)                                                                                                                    cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGaaa                           <  8:156737/149‑1 (MQ=255)                                                                                                                    cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGaaa                           <  6:204886/149‑1 (MQ=255)                                                                                                                    cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGaaa                           <  4:291498/149‑1 (MQ=255)                                                                                                                             ccATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTtgtg                  <  5:120356/149‑1 (MQ=255)                                                                                                                              cATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTgg                 >  1:226756/1‑149 (MQ=255)                                                                                                                                ttGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATGACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTc               >  3:279370/1‑149 (MQ=255)                                                                                                                                  gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATaa                                                                <  1:232000/98‑1 (MQ=18)                                                                                                                                  gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATaa                                                                >  2:232000/1‑98 (MQ=18)                                                                                                                                  gCCGGATAAGGCGGTCACGCCGCATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGCTCAGCGGAGCGCATCAGGCAGTCTGGCGTTGGTCATaa                                                                <  1:236603/98‑1 (MQ=18)                                                                                                                                   ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGtt            <  1:35665/149‑1 (MQ=255)                                                                                                                                     ggATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGc          >  6:299599/1‑149 (MQ=255)                                                                                                                                       ataaggcgTTCACGCCGCATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCa        >  8:5783/5‑149 (MQ=37)                                                                                                                                        tGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGTGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCAt       <  1:277368/149‑1 (MQ=255)                                                                                                                                             cgTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATtggtg  <  4:14979/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 CTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTG  >  NC_000913/339907‑340195

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