Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,823,700	C→T	S434L (TCG→TTG)	spoT →	bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,823,700	0	C	T	100.0%	322.6 / NA	80	S434L (TCG→TTG)	spoT	bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (39/41);  total (39/41)

GGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTC  >  NC_000913/3823554‑3823844                                                                                                                                                   |                                                                                                                                                 ggATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGc                                                                                                                                                >  2:670447/1‑149 (MQ=255) ggATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAg                                                                                                                                              <  2:898214/151‑1 (MQ=255)     gagaTTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGc                                                                                                                                          <  2:947465/151‑1 (MQ=255)      agaTTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCt                                                                                                                                         >  1:381873/1‑151 (MQ=255)           tACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTAcc                                                                                                                                     >  2:311676/1‑150 (MQ=255)              gTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCg                                                                                                                                 >  2:295929/1‑151 (MQ=255)               ttttCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCCATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCgg                                                                                                                                >  2:471906/1‑151 (MQ=255)                 ttCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTc                                                                                                                              >  2:625454/1‑151 (MQ=255)                   cacaCCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCa                                                                                                                             <  2:957488/150‑1 (MQ=255)                    acacCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCaa                                                                                                                            >  1:601166/1‑150 (MQ=255)                       ccGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCg                                                                                                                        <  2:628342/151‑1 (MQ=255)                        cGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGt                                                                                                                       >  1:209167/1‑151 (MQ=255)                         ggAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGtt                                                                                                                      >  2:683511/1‑151 (MQ=255)                         ggAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGtt                                                                                                                      <  2:906589/151‑1 (MQ=255)                          gAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTg                                                                                                                     <  2:205871/151‑1 (MQ=255)                               gcgcATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATc                                                                                                                >  1:387162/1‑151 (MQ=255)                                 gcATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCAt                                                                                                              >  1:60040/1‑151 (MQ=255)                                   aTTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCAtt                                                                                                             <  2:362389/150‑1 (MQ=255)                                   aTTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTa                                                                                                            <  2:288481/151‑1 (MQ=255)                                      gTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCg                                                                                                         <  1:719173/151‑1 (MQ=255)                                         gAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTc                                                                                                      <  2:365397/151‑1 (MQ=255)                                           gCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCg                                                                                                    <  1:558936/151‑1 (MQ=255)                                              gcctgccGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCgg                                                                                                   <  1:382698/149‑1 (MQ=255)                                                 tgccGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGgcgc                                                                                               >  2:279347/1‑150 (MQ=255)                                                    cGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGc                                                                                           <  2:818046/151‑1 (MQ=255)                                                    cGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGc                                                                                           <  1:973300/151‑1 (MQ=255)                                                     ggTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGcc                                                                                          <  1:100429/151‑1 (MQ=255)                                                     ggTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGcc                                                                                          >  1:584401/1‑151 (MQ=255)                                                             gCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGcc                                                                                  >  2:125135/1‑151 (MQ=255)                                                               ccGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGc                                                                                >  2:369738/1‑151 (MQ=255)                                                                 gTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCtt                                                                              <  1:899732/151‑1 (MQ=255)                                                                   cGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTg                                                                             <  2:564456/150‑1 (MQ=255)                                                                    gACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGc                                                                           >  1:88261/1‑151 (MQ=255)                                                                     aCTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGc                                                                           >  2:842506/1‑150 (MQ=255)                                                                     aCTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCt                                                                          >  1:973560/1‑151 (MQ=255)                                                                          gCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAAc                                                                      >  2:75412/1‑150 (MQ=255)                                                                             tATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACttt                                                                   <  2:28059/150‑1 (MQ=255)                                                                                gCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCg                                                               >  2:81406/1‑151 (MQ=255)                                                                                 cAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGt                                                              >  2:869689/1‑151 (MQ=255)                                                                                   gTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGtt                                                             <  2:31241/150‑1 (MQ=255)                                                                                   gTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGtt                                                             >  2:482085/1‑150 (MQ=255)                                                                                   gTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGtt                                                             >  2:766007/1‑150 (MQ=255)                                                                                     gCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGc                                                          >  1:668606/1‑151 (MQ=255)                                                                                      cATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCt                                                         >  1:553683/1‑151 (MQ=255)                                                                                       aTACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTc                                                        <  2:550164/151‑1 (MQ=255)                                                                                         aCCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCg                                                       >  1:753047/1‑150 (MQ=255)                                                                                         aCCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGa                                                      <  1:965238/151‑1 (MQ=255)                                                                                          ccGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGa                                                      <  1:768800/150‑1 (MQ=255)                                                                                               atCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgc                                                <  1:540207/151‑1 (MQ=255)                                                                                               atCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcg                                                 <  1:441781/150‑1 (MQ=255)                                                                                                 cGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcg                                               <  2:718735/150‑1 (MQ=255)                                                                                                   gTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcgcgc                                            >  2:890787/1‑151 (MQ=255)                                                                                                    tcttGCCTGCGCGGGCGCAAGCGTTGACGGCCAGCCTTATCCGCTGTTGCAGCCGGTTACCAGCGGCCAAACCGTTGGCATCATTACCGCTCCGGGCGCTAGCCAGAATGCCGCCTGGCCGTACTTTGTCGATACCTCGCAAgcgcgcg                                             >  2:743925/4‑149 (MQ=255)                                                                                                       tGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCaa                                         >  1:68974/1‑150 (MQ=255)                                                                                                        gCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCaaa                                        >  2:799175/1‑150 (MQ=255)                                                                                                        gCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCaaa                                        <  1:180757/150‑1 (MQ=255)                                                                                                            gCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTc                                    <  2:799827/150‑1 (MQ=255)                                                                                                              gTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTc                                 >  2:489407/1‑151 (MQ=255)                                                                                                              gTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTc                                 >  1:448301/1‑151 (MQ=255)                                                                                                                 ggCGCACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGt                              <  2:381867/151‑1 (MQ=255)                                                                                                                  gcgcACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGtt                             <  1:764441/151‑1 (MQ=255)                                                                                                                    gcACGCGTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTg                            >  2:952597/1‑150 (MQ=255)                                                                                                                       cgcgTTGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTg                         >  2:560610/1‑150 (MQ=255)                                                                                                                           ttGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGaaaa                     <  2:682864/150‑1 (MQ=255)                                                                                                                            tGACCGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGaaaaa                    <  2:419166/150‑1 (MQ=255)                                                                                                                               ccGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTc                <  2:443497/151‑1 (MQ=255)                                                                                                                               ccGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTc                <  1:60372/151‑1 (MQ=255)                                                                                                                                cGCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCa               <  1:840129/151‑1 (MQ=255)                                                                                                                                 gCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCaa              <  2:752611/151‑1 (MQ=255)                                                                                                                                 gCCAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCaa              <  2:115949/151‑1 (MQ=255)                                                                                                                                  ccAGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCa               >  2:247697/1‑149 (MQ=255)                                                                                                                                    aGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCg           >  1:53963/1‑151 (MQ=255)                                                                                                                                    aGCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCg           <  2:238776/151‑1 (MQ=255)                                                                                                                                     gCCTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGt          <  2:88258/151‑1 (MQ=255)                                                                                                                                       cTTACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtg         >  2:910924/1‑150 (MQ=255)                                                                                                                                        ttACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtga        >  1:818668/1‑150 (MQ=255)                                                                                                                                         tACCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgat       <  1:629887/150‑1 (MQ=255)                                                                                                                                          aCCCGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgatga     <  1:167708/151‑1 (MQ=255)                                                                                                                                           cccGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTgatgat    <  1:830489/151‑1 (MQ=255)                                                                                                                                             cGCTGTTGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTc  >  1:137313/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTC  >  NC_000913/3823554‑3823844

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