Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,452,005	G→A	T316T (ACG→ACA)	ytfR →	putative sugar ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,452,005	0	G	A	100.0%	240.9 / NA	63	T316T (ACG→ACA)	ytfR	putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/27);  total (36/27)

TCAAAAATTACGGCAAAAAAGGAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACGGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGG  >  NC_000913/4451857‑4452150                                                                                                                                                     |                                                                                                                                                  tCAAAAATTACGGCAAAAAAGGAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGc                                                                                                                                                 >  2:953058/1‑151 (MQ=255)   aaaaaTTACGGCAAAAAAGGAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGt                                                                                                                                               <  2:440034/151‑1 (MQ=255)      aaTTACGGCAAAAAAGGAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTg                                                                                                                                             >  2:90979/1‑150 (MQ=255)         tACGGCAAAAAAGGAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATc                                                                                                                                          >  1:372744/1‑150 (MQ=255)            ggCAAAAAAGGAACGATCGCCCCGTTTGATCTGGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAg                                                                                                                                      <  2:543636/151‑1 (MQ=255)             gCAAAAAAGGAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAgg                                                                                                                                     <  1:90601/151‑1 (MQ=255)               aaaaaaGGAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCa                                                                                                                                   <  2:762392/151‑1 (MQ=255)                   aaGGAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAAcc                                                                                                                               <  2:458397/151‑1 (MQ=255)                    aGGAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCg                                                                                                                              >  2:264752/1‑151 (MQ=255)                     ggAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATGGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGc                                                                                                                             <  2:743565/151‑1 (MQ=255)                       aaCGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCa                                                                                                                            >  1:140325/1‑150 (MQ=255)                          gATCGCCCCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCaaaa                                                                                                                         <  1:899611/150‑1 (MQ=255)                          gATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAAc                                                                                                                        >  2:363572/1‑151 (MQ=255)                               cACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCg                                                                                                                   >  1:853510/1‑151 (MQ=255)                                aCCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCg                                                                                                                   >  2:366335/1‑150 (MQ=255)                                 ccGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGa                                                                                                                  >  2:708730/1‑150 (MQ=255)                                 ccGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGAt                                                                                                                 <  2:609674/151‑1 (MQ=255)                                 ccGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGAt                                                                                                                 >  2:433905/1‑151 (MQ=255)                                   gTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATc                                                                                                                <  1:127030/150‑1 (MQ=255)                                        aTCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCcac                                                                                                           >  1:945199/1‑150 (MQ=255)                                             gaagaaCCCCCCGGGGAGTTCGTCGGTGGGGCTGGATTCCTGGGTTCAGGACGTGCCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCATATCTgg                                                                                                     <  1:384939/151‑1 (MQ=255)                                                  aCGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCg                                                                                                <  1:675843/151‑1 (MQ=255)                                                   cGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTCCCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCg                                                                                                <  1:811397/150‑1 (MQ=255)                                                    gCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGt                                                                                              >  1:51745/1‑151 (MQ=255)                                                     ccccGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTCCCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGt                                                                                              <  1:26436/150‑1 (MQ=255)                                                     ccccGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGt                                                                                              >  2:316797/1‑150 (MQ=255)                                                       ccGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACt                                                                                           >  1:582424/1‑151 (MQ=255)                                                        cGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACt                                                                                           >  1:463846/1‑150 (MQ=255)                                                                tcgtcgGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTgg                                                                                  >  1:726268/1‑151 (MQ=255)                                                                    cgGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCtt                                                                               <  1:515490/150‑1 (MQ=255)                                                                      gTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTg                                                                            <  1:299810/151‑1 (MQ=255)                                                                          ggctggATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCg                                                                        >  1:4153/1‑151 (MQ=255)                                                                               gATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGa                                                                      >  1:494688/1‑148 (MQ=255)                                                                                               aggtcccGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAAcc                                                         <  1:502625/140‑1 (MQ=255)                                                                                                 gTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGcat                                                 >  2:781359/1‑151 (MQ=255)                                                                                                    ccGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGcatca                                               <  1:995912/150‑1 (MQ=255)                                                                                                     cGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCatcatc                                             >  1:187866/1‑151 (MQ=255)                                                                                                      gAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCg                                            >  2:649959/1‑151 (MQ=255)                                                                                                        aaCCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCt                                          >  1:371664/1‑151 (MQ=255)                                                                                                        aaCCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCt                                          >  2:365988/1‑151 (MQ=255)                                                                                                        aaCCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCt                                          <  1:778031/151‑1 (MQ=255)                                                                                                            gccgAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTgcc                                       <  2:976594/150‑1 (MQ=255)                                                                                                                aaGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTc                                  >  2:583012/1‑151 (MQ=255)                                                                                                                aaGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTc                                  <  2:141592/151‑1 (MQ=255)                                                                                                                     aTCTTCGGTATCAAACCTGCTGACACCAGCACAGCTTTGTTCAAAGCCAAACCGCAAAATCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTAGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTgc                                        >  1:790809/1‑140 (MQ=255)                                                                                                                         tCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGc                             <  1:811295/147‑1 (MQ=255)                                                                                                                         tCGGTATCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGc                             <  2:322891/147‑1 (MQ=255)                                                                                                                              aTCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAat                     >  2:975357/1‑150 (MQ=255)                                                                                                                              aTCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAat                     <  2:746018/150‑1 (MQ=255)                                                                                                                               tCAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAatat                   >  1:193284/1‑151 (MQ=255)                                                                                                                                cAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAatat                   <  1:699705/150‑1 (MQ=255)                                                                                                                                cAAACCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATatc                  >  1:272817/1‑151 (MQ=255)                                                                                                                                  aaCCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATatcat                >  2:422477/1‑151 (MQ=255)                                                                                                                                   aCCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATatcatc               >  2:335266/1‑151 (MQ=255)                                                                                                                                   aCCTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATatcatc               >  1:674996/1‑151 (MQ=255)                                                                                                                                    ccTGCTGACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATatcatc               <  2:606231/150‑1 (MQ=255)                                                                                                                                     ctgctgACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCt             >  1:177989/1‑151 (MQ=255)                                                                                                                                     ctgctgACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCt             >  1:564255/1‑151 (MQ=255)                                                                                                                                       gctgACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCg           <  1:823570/151‑1 (MQ=255)                                                                                                                                        ctgACAGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGc          >  1:241661/1‑151 (MQ=255)                                                                                                                                             aGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTcc     >  2:850247/1‑151 (MQ=255)                                                                                                                                             aGCGGCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTcc     >  2:527244/1‑151 (MQ=255)                                                                                                                                                ggCACAGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAgg  <  1:122376/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TCAAAAATTACGGCAAAAAAGGAACGATCGCACCGTTTGATCTCGAAGTACGCCCCGGCGAGATCGTCGGTCTGGCTGGATTGCTGGGATCAGGACGTACCGAAACCGCCGAAGTGATCTTCGGTATCAAACCTGCTGACAGCGGCACGGCGTTGATCAAAGGCAAACCGCAAAACCTGCGATCGCCACATCAGGCTTCGGTACTGGGCATTGGCTTCTGCCCGGAAGACAGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGG  >  NC_000913/4451857‑4452150

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