Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	340,058	T→C	6.2%	intergenic (+315/‑107)	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,058	0	T	C	6.2%	185.3 / 3.1	64	intergenic (+315/‑107)	yahG/yahI	DUF1116 family protein/carbamate kinase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (24/36);  new base C (2/2);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

ATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTG  >  NC_000913/339910‑340195                                                                                                                                                     |                                                                                                                                          attaGCGTATCGGATCAGGCATTTGCCCCCACTGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCtt                                                                                                                                           <  8:505/146‑1 (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:140385/149‑1 (MQ=255)              aTCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGa                                                                                                                              <  1:247004/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                                                                       <  6:14699/149‑1 (MQ=255)                                                                     aGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTg                                                                       <  3:130975/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                           <  4:291498/149‑1 (MQ=255)                                                                                                                 cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGaaa                           <  6:204886/149‑1 (MQ=255)                                                                                                                 cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGaaa                           <  2:74982/149‑1 (MQ=255)                                                                                                                 cATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGaaa                           <  8:156737/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                                                                >  2:232000/1‑98 (MQ=18)                                                                                                                               gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATaa                                                                <  1:232000/98‑1 (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)                                                                                                                                               aCGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACAt                               >  6:240406/1‑115 (MQ=255)                                                                                                                                               aCGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACAt                               <  5:240406/115‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                          ATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTG  >  NC_000913/339910‑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: ‑