Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,899,249	T→G	T49P (ACC→CCC)	yieK ←	putative 6‑phosphogluconolactonase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,899,249	0	T	G	100.0%	248.5 / NA	63	T49P (ACC→CCC)	yieK	putative 6‑phosphogluconolactonase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (31/32);  total (31/32)

TTCTCTTCTTTGATCCCCGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGTGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCAT  >  NC_000913/3899104‑3899393                                                                                                                                                  |                                                                                                                                                 ttCTCTTCTTTGATCCCCGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGggggtg                                                                                                                                              >  1:642576/1‑150 (MQ=255)    tcttctTTGATCCCCGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACC‑TTAAACAGGGGGGTgggg                                                                                                                                           >  1:252559/1‑149 (MQ=255)     cttcttTGATCCCCGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAgg                                                                                                                                           <  1:405635/149‑1 (MQ=255)     cttcttTGATCCCCGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGt                                                                                                                                          >  2:375650/1‑150 (MQ=255)          ttGATCCCCGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTAtt                                                                                                                                       <  1:1017053/148‑1 (MQ=255)            gATCCCCGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATAc                                                                                                                                  <  1:865047/151‑1 (MQ=255)               ccccGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATg                                                                                                                               <  2:1033630/151‑1 (MQ=255)               ccccGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATg                                                                                                                               >  2:680295/1‑151 (MQ=255)                  cGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCt                                                                                                                            >  2:31149/1‑151 (MQ=255)                    cAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCttt                                                                                                                          >  1:1077808/1‑151 (MQ=255)                    cAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCttt                                                                                                                          <  1:372684/151‑1 (MQ=255)                     aGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTg                                                                                                                         >  2:1038373/1‑151 (MQ=255)                         gTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCCGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCg                                                                                                                     <  1:915448/151‑1 (MQ=255)                                 cAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACcg                                                                                                             <  1:379842/151‑1 (MQ=255)                                 cAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACcg                                                                                                             <  1:798465/151‑1 (MQ=255)                                  aGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACcgg                                                                                                            >  1:109993/1‑151 (MQ=255)                                  aGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACcgg                                                                                                            >  1:69228/1‑151 (MQ=255)                                     ttACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACcggcgg                                                                                                         >  1:881221/1‑151 (MQ=255)                                       aCGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACcggcgg                                                                                                         >  2:600119/1‑149 (MQ=255)                                       aCGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGt                                                                                                        >  2:468683/1‑150 (MQ=255)                                       aCGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTa                                                                                                       <  2:868698/151‑1 (MQ=255)                                          cAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAAtt                                                                                                    <  2:1017997/151‑1 (MQ=255)                                                  tAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTTCCAGGtt                                                                                            >  2:226463/1‑151 (MQ=255)                                                   aaTCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCTTACCAGGGCTTACCGTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGtt                                                                                            <  2:550489/150‑1 (MQ=255)                                                      cGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAac                                                                                         <  2:815999/150‑1 (MQ=255)                                                           cGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACacgacg                                                                                   <  1:1054642/151‑1 (MQ=255)                                                           cGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACacgac                                                                                    <  2:276399/150‑1 (MQ=255)                                                             ccTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGAcgcg                                                                                 <  1:679494/151‑1 (MQ=255)                                                             ccTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAAAACGAcgc                                                                                  >  1:1045596/1‑150 (MQ=255)                                                                 ctcCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGAAATTGGCAGGATAACacgacgacgctt                                                                             >  1:618821/1‑145 (MQ=255)                                                                   cccTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTc                                                                            <  1:117578/150‑1 (MQ=255)                                                                   cccTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCg                                                                           >  1:391493/1‑151 (MQ=255)                                                                    ccTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGa                                                                          <  1:123082/151‑1 (MQ=255)                                                                          tGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGAGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGGAATTGCCAGGGTAACACGACGGGTCTTCGACata                                                                      >  1:778179/1‑149 (MQ=255)                                                                           gCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACAtata                                                                    >  1:767099/1‑150 (MQ=255)                                                                             cgcgAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAAcc                                                                 <  1:79607/151‑1 (MQ=255)                                                                                    tGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCaa                                                               >  1:927816/1‑146 (MQ=255)                                                                                    tGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAAATGCCAGGGTAACACGACGCGTCTTTGACATATAAACaa                                                               >  1:837718/1‑146 (MQ=255)                                                                                      gAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACCCGTCTTCGACATATAAACAAGCa                                                            >  1:1009736/1‑147 (MQ=255)                                                                                         tttCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCCGGGTAAAACGACGCGTCTTCGACATATAACCAAACAGATgg                                                       >  2:873008/1‑149 (MQ=255)                                                                                           tcatcaAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGgtgt                                                    <  1:22442/150‑1 (MQ=255)                                                                                           tcatcaAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGGAATTGCCAGGTTAAAACGACGCGTCTTAGACATATAAACAAGAAGATGgtgt                                                    >  1:151334/1‑150 (MQ=255)                                                                                            catcaAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAAACAAGCAGATGGtgtg                                                   >  1:206726/1‑150 (MQ=255)                                                                                               caAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCcgc                                               <  2:529101/151‑1 (MQ=255)                                                                                                aaaaTTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGGTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGGGTGcc                                                 >  1:1061430/1‑148 (MQ=255)                                                                                                 aaaTTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCcgc                                               >  2:311691/1‑149 (MQ=255)                                                                                                        aGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGc                                       <  2:444779/150‑1 (MQ=255)                                                                                                           aaTAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCAt                                   <  1:285541/151‑1 (MQ=255)                                                                                                              aGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTc                                <  1:421911/151‑1 (MQ=255)                                                                                                                cAGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAAAACCACGCGTCTTCGACATATAACCAAGGAGATGGTGGGGCGCGACACGGCCCATTTTCt                              >  2:813053/1‑151 (MQ=255)                                                                                                                 aGTTATCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTg                             <  1:154316/151‑1 (MQ=255)                                                                                                                     aTCGTACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTaa                         <  2:873157/151‑1 (MQ=255)                                                                                                                         tACCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCtt                     >  1:424311/1‑151 (MQ=255)                                                                                                                          aCCAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTc                    <  2:90897/151‑1 (MQ=255)                                                                                                                           ccAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCg                   <  1:645067/151‑1 (MQ=255)                                                                                                                           ccAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCAACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCg                   <  2:336561/151‑1 (MQ=255)                                                                                                                           ccAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAAGTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGAGGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTc                    <  2:1025598/150‑1 (MQ=255)                                                                                                                            cAGGGCTTACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCgg                  >  1:637093/1‑151 (MQ=255)                                                                                                                                   tACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATg             <  1:522683/149‑1 (MQ=255)                                                                                                                                   tACCTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGa            <  1:623576/150‑1 (MQ=255)                                                                                                                                      cTTTAACCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGGAATTGCCCGGCTAACACAACCCCTGTTCGACAAATAAACAAAGAGATGGGGTGCCGCGACACGGCTCAATTCCTGGGAAATTTCGGAAATGGTTa         >  1:128563/1‑150 (MQ=255)                                                                                                                                          aaCCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTc    <  1:313551/151‑1 (MQ=255)                                                                                                                                           aCCAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTTCTGGTAATCTTCGGTAATGATTAAttt     >  1:1005940/1‑149 (MQ=255)                                                                                                                                            ccAGGGGGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGGTAAAACGACGCGGCTTCGACATATAAACAAACAGAAGGGGTGCCGCGACACGGCTCATTTCCTGGGAATCTTTGGGAATGATTAATTTTAt  >  2:394860/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TTCTCTTCTTTGATCCCCGCAGGGGTGAAAAACAGATTACGCAGATTGGTAATCGTTACGCCTTCTCCCTCTTTGCCGCGAAATGGAATTTCATCAAAATTATAGAAATAGCAGTTATCGTACCAGGGCTTACCTTTAACCAGGGTGGTGAGGTATTCATACATGCCTTTGGGCGTGCTACCGGCGGTAATTGCCAGGTTAACACGACGCGTCTTCGACATATAACCAAGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCAT  >  NC_000913/3899104‑3899393

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