Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,306,037	A→G	G187G (GGT→GGC)	yjcS ←	alkyl sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,306,037	0	A	G	98.7%	289.3 / ‑6.0	78	G187G (GGT→GGC)	yjcS	alkyl sulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (37/40);  total (37/41)
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

TAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCC  >  NC_000913/4305903‑4306168                                                                                                                                       |                                                                                                                                    tAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATa                                                                                                                                >  2:898746/1‑140 (MQ=255)   gagCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGt                                                                                                                              >  2:97722/1‑140 (MQ=255)     gcgcACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTgg                                                                                                                            >  2:816730/1‑140 (MQ=255)      cgcgcgGCGGCTATTGTTGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGt                                                                                                                           <  1:1046609/140‑1 (MQ=255)      cgcCCGGCGGCTCATGATTTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGt                                                                                                                           <  1:651753/140‑1 (MQ=255)      cgcACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGt                                                                                                                           >  1:1021957/1‑140 (MQ=255)      cgcACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGt                                                                                                                           >  1:600243/1‑140 (MQ=255)        cACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCg                                                                                                                         >  2:563046/1‑140 (MQ=255)         aCGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGt                                                                                                                           >  2:714831/1‑137 (MQ=255)         aCGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGt                                                                                                                           <  1:714831/137‑1 (MQ=255)         aCGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCgg                                                                                                                        >  1:943268/1‑140 (MQ=255)          cggcggCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGgt                                                                                                                       <  1:25669/140‑1 (MQ=255)          cggcggCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGgt                                                                                                                       <  2:647682/140‑1 (MQ=255)           ggcggcTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGt                                                                                                                           <  1:956742/135‑1 (MQ=255)           ggcggcTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGt                                                                                                                           >  2:956742/1‑135 (MQ=255)             cggcTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGtgtg                                                                                                                    >  1:1024445/1‑140 (MQ=255)                cTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCt                                                                                                                 <  1:736373/140‑1 (MQ=255)                 tCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCtg                                                                                                                <  1:816730/140‑1 (MQ=255)                 tCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCtg                                                                                                                <  1:898746/140‑1 (MQ=255)                 tCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCtg                                                                                                                <  1:191017/140‑1 (MQ=255)                  cATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCtgt                                                                                                               <  2:638603/140‑1 (MQ=255)                  cATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCtgt                                                                                                               <  1:563046/140‑1 (MQ=255)                  cATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCtgt                                                                                                               <  1:271648/140‑1 (MQ=255)                  cATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCtgt                                                                                                               <  2:1004395/140‑1 (MQ=255)                     gatgTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAg                                                                                                            >  2:969190/1‑140 (MQ=255)                       gttgaCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAgtg                                                                                                          <  1:815617/136‑1 (MQ=255)                          tACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAg                                                                                                       >  1:622407/1‑140 (MQ=255)                           aCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGa                                                                                                      >  2:251405/1‑140 (MQ=255)                            cccGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCtg                                                                                                                >  2:89607/1‑129 (MQ=255)                            cccGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCtg                                                                                                                <  1:89607/129‑1 (MQ=255)                            cccGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGAt                                                                                                     >  2:7340/1‑140 (MQ=255)                               gCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTTGGTGTGCCTGTGAGTGTAGATAAc                                                                                                  >  2:965415/1‑140 (MQ=255)                               gCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAAc                                                                                                  >  2:902562/1‑140 (MQ=255)                                ccAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACg                                                                                                 >  2:322318/1‑140 (MQ=255)                                 cAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACgg                                                                                                >  1:274165/1‑140 (MQ=255)                                 cAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACgg                                                                                                >  1:587782/1‑140 (MQ=255)                                  aGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGc                                                                                               <  1:7340/140‑1 (MQ=255)                                      cGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGcaaca                                                                                           >  1:572447/1‑140 (MQ=255)                                                 aTGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTg                                                                                >  1:30981/1‑140 (MQ=255)                                                     cTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGa                                                                            <  1:902562/140‑1 (MQ=255)                                                     cTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGa                                                                            <  1:244995/140‑1 (MQ=255)                                                      ttCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGAc                                                                           <  2:1024445/140‑1 (MQ=255)                                                      acagTCCATAAAGCCTGCAGGGGCAACCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGAc                                                                           <  1:743126/137‑1 (MQ=255)                                                       tCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACg                                                                          <  2:457424/140‑1 (MQ=255)                                                         gTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGAt                                                                        <  2:401402/140‑1 (MQ=255)                                                            cATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCt                                                                     <  2:162547/140‑1 (MQ=255)                                                              tAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGa                                                                            >  1:61057/1‑131 (MQ=255)                                                              tAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGa                                                                            <  2:61057/131‑1 (MQ=255)                                                                 aGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAg                                                                <  2:274165/140‑1 (MQ=255)                                                                 aGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAg                                                                <  2:587782/140‑1 (MQ=255)                                                                   ccTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTa                                                              >  2:92022/1‑140 (MQ=255)                                                                   ccTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTa                                                              >  2:487355/1‑140 (MQ=255)                                                                     tacAGGGGCGATCACCTGAACTAAGCCGGGTGTGACATCGGGCTCAGGGATGAGGGCTTTCACGCCACCGACGCGGCCGGCGAGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGAGGCTGGAAgctgt                                                            >  2:743126/3‑136 (MQ=255)                                                                     tGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACg                                                                          >  2:587859/1‑126 (MQ=255)                                                                     tGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACg                                                                          <  1:587859/126‑1 (MQ=255)                                                                               aTCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGc                                                  >  1:295255/1‑140 (MQ=255)                                                                                   ccTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCt                                              <  2:532081/140‑1 (MQ=255)                                                                                    cTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCtt                                             >  1:314837/1‑140 (MQ=255)                                                                                                 gATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAggcgg                                >  2:507886/1‑140 (MQ=255)                                                                                                    ttAACATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGt                             <  1:969190/140‑1 (MQ=255)                                                                                                        cATCGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCAcc                         >  1:363665/1‑140 (MQ=255)                                                                                                           cGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGc                      <  2:600243/140‑1 (MQ=255)                                                                                                           cGGCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGc                      <  1:507886/140‑1 (MQ=255)                                                                                                            ggCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCg                     <  2:622407/140‑1 (MQ=255)                                                                                                            ggCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCg                     >  1:615149/1‑140 (MQ=255)                                                                                                             gCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCgg                    >  2:57398/1‑140 (MQ=255)                                                                                                             gCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCgg                    >  1:838331/1‑140 (MQ=255)                                                                                                             gCTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCgg                    <  1:433340/140‑1 (MQ=255)                                                                                                              cTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCggg                   <  2:807241/140‑1 (MQ=255)                                                                                                              cTTCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCggg                   <  2:363665/140‑1 (MQ=255)                                                                                                                tCAGAGATAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTc                 <  2:943268/140‑1 (MQ=255)                                                                                                                   gagaTAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGAt              <  2:838331/140‑1 (MQ=255)                                                                                                                   gagaTAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGAGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGAt              <  1:487355/140‑1 (MQ=255)                                                                                                                     gaTAATGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATaa            <  1:322318/140‑1 (MQ=255)                                                                                                                        aaTGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGGCGATAACaa         >  2:112172/1‑140 (MQ=255)                                                                                                                         aTGCCTTTCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaat        <  2:515360/140‑1 (MQ=255)                                                                                                                              tttCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGc   >  1:267571/1‑140 (MQ=255)                                                                                                                              tttCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGc   >  1:307145/1‑140 (MQ=255)                                                                                                                               ttCACGCCGCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGcc  >  1:547623/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                    TAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCC  >  NC_000913/4305903‑4306168

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑