Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,146,137	A→G	6.1%	S153S (AGT→AGC)	yceF ←	m(7)GTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,146,137	0	A	G	6.1%	192.5 / 2.7	65	S153S (AGT→AGC)	yceF	m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base G (2/2);  total (33/32)
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

CACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTC  >  NC_000913/1145993‑1146282                                                                                                                                                 |                                                                                                                                                  cacaTCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCActct                                                                                                                                               >  1:130115/1‑149 (MQ=255) cacaTCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCActct                                                                                                                                               >  7:212409/1‑149 (MQ=255) cacaTCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCActct                                                                                                                                               <  7:336709/149‑1 (MQ=255)  acaTCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCtt                                                                                                                                              <  1:66966/149‑1 (MQ=255)     tCCGGCGCAAAAAACTATCCCATCAGCGGGCTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTaaa                                                                                                                                           <  8:212409/149‑1 (MQ=255)       cGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGc                                                                                                                                         <  8:163690/149‑1 (MQ=255)            cAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTAccc                                                                                                                                    <  2:213560/149‑1 (MQ=255)             aaaaaaCTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCcg                                                                                                                                   >  3:222256/1‑149 (MQ=255)                 aaCTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCa                                                                                                                               <  8:105192/149‑1 (MQ=255)                    tATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTg                                                                                                                            >  2:266798/1‑149 (MQ=255)                           tCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAAAGCCAAATCCATCACTCTTAAAGCTACCCGCGCAGTGCAGGGGa                                                                                                                     >  8:148823/1‑149 (MQ=255)                              gCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGc                                                                                                                  >  1:316849/1‑149 (MQ=255)                                   tttttCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCttt                                                                                                             <  5:199861/149‑1 (MQ=255)                                       tcccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACg                                                                                                          >  4:120949/1‑148 (MQ=255)                                        cccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCa                                                                                                        >  7:36466/1‑149 (MQ=255)                                         ccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCAc                                                                                                       >  5:121647/1‑149 (MQ=255)                                         ccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCAc                                                                                                       >  3:41685/1‑149 (MQ=255)                                                gACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTa                                                                                                >  1:105179/1‑149 (MQ=255)                                                  cGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTAACGCACGTAATTATc                                                                                              <  6:142545/149‑1 (MQ=255)                                                            gACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCGCTCGTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCAGcc                                                                                    >  3:224123/1‑149 (MQ=255)                                                            gACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCACCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTAGCTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGcc                                                                                    >  3:219464/1‑149 (MQ=255)                                                             acacAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCt                                                                                   >  3:181305/1‑149 (MQ=255)                                                               acaGCGCGATTAGCGGTAAACCCACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCGCGTAATTATCAATCTCCGCCTCg                                                                                 <  2:220954/149‑1 (MQ=255)                                                               acaGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAgggg                                                                                                                      <  7:276348/112‑1 (MQ=255)                                                               acaGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAgggg                                                                                                                      >  8:276348/1‑112 (MQ=255)                                                                   cgcgATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCa                                                                             <  4:181305/149‑1 (MQ=255)                                                                      gATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGt                                                                          >  2:44849/1‑149 (MQ=255)                                                                           gCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACg                                                                      >  5:117708/1‑148 (MQ=255)                                                                           gCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACg                                                                      <  3:120949/148‑1 (MQ=255)                                                                              gTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGTCGGaaa                                                                  <  6:121647/149‑1 (MQ=255)                                                                               tAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAAt                                                                 <  4:41685/149‑1 (MQ=255)                                                                               tAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAAt                                                                 <  2:105179/149‑1 (MQ=255)                                                                                   ccaaccaGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACGCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGAc                                                             <  7:148823/149‑1 (MQ=255)                                                                                   ccaaccaGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGAc                                                             <  8:36466/149‑1 (MQ=255)                                                                                   ccaaccaGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGAc                                                             >  1:27916/1‑149 (MQ=255)                                                                                   ccaaccaGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGAc                                                             <  2:316849/149‑1 (MQ=255)                                                                                   ccaaccaGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCAAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGAc                                                             <  2:52869/149‑1 (MQ=255)                                                                                    caaccaGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACg                                                            >  2:116940/1‑149 (MQ=255)                                                                                         aGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCaaa                                                       <  4:222256/149‑1 (MQ=255)                                                                                         aGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCaa                                                        >  8:186482/1‑148 (MQ=255)                                                                                          gCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGct                                                                                                                 >  1:319932/1‑90 (MQ=255)                                                                                          gCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGct                                                                                                                 <  2:319932/90‑1 (MQ=255)                                                                                          gCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCaa                                                        >  2:60296/1‑147 (MQ=255)                                                                                                  ggATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCAc                                              <  2:130115/149‑1 (MQ=255)                                                                                                    aTCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACtt                                            >  1:110102/1‑149 (MQ=255)                                                                                                      cACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTCCCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGCCTCCACTTCt                                          >  3:111761/1‑149 (MQ=255)                                                                                                       aCGCCCCTCTAAAAGCTCAAACAGCTTAATGCCAAATCCTTCGCTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCt                                          <  4:219464/148‑1 (MQ=255)                                                                                                       aCGCCACTCTAAACGCTCAAACATCGTAATGACAAATCCTTCGCTCTTAAAGCTACCAGCGCAGTGCAGTGGATGCTCTTTACGCACGTAATTATCAATCTCCGCATCGCTCAGGTGACTGAAATGGACGTCAAAAGGCTCCACTTCt                                          <  4:224123/148‑1 (MQ=255)                                                                                                         gCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGtt                                       <  1:266798/149‑1 (MQ=255)                                                                                                            ccTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAAt                                                                 <  6:79276/120‑1 (MQ=255)                                                                                                            ccTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAAt                                                                 >  5:79276/1‑120 (MQ=255)                                                                                                             ctctAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACg                                                                      <  2:89527/114‑1 (MQ=255)                                                                                                             ctctAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACg                                                                      >  1:89527/1‑114 (MQ=255)                                                                                                               ctAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGAc                                                             >  4:168854/1‑121 (MQ=255)                                                                                                               ctAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGAc                                                             <  3:168854/121‑1 (MQ=255)                                                                                                                tAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAg                                  <  1:199437/147‑1 (MQ=255)                                                                                                                 aaaCGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAg                                  >  3:243852/1‑146 (MQ=255)                                                                                                                 aaaCGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAg                                  <  4:243852/146‑1 (MQ=255)                                                                                                                  aaCGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAgg                                                    <  4:276201/127‑1 (MQ=255)                                                                                                                  aaCGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAgg                                                    >  3:276201/1‑127 (MQ=255)                                                                                                                       tCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACACGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCAtt                         >  7:102879/1‑149 (MQ=255)                                                                                                                                   aTGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTaa              <  8:102879/148‑1 (MQ=255)                                                                                                                                    tGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTaaa             >  5:164269/1‑148 (MQ=255)                                                                                                                                       cAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAaacaac         >  3:215077/1‑149 (MQ=255)                                                                                                                                              ttCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGccagtc  >  8:15969/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTC  >  NC_000913/1145993‑1146282

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