Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,823,105	C→A	R236S (CGC→AGC)	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,105	0	C	A	100.0%	229.1 / NA	60	R236S (CGC→AGC)	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 A (26/34);  total (26/34)

TTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTG  >  NC_000913/3822959‑3823248                                                                                                                                                   |                                                                                                                                                ttGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTTGTagcg                                                                                                                                              <  2:804706/150‑1 (MQ=255) ttGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCAGCGTCAGTGGTagcg                                                                                                                                              <  1:276826/150‑1 (MQ=255)  tGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTagcg                                                                                                                                              >  1:316029/1‑149 (MQ=255)    aGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCgaga                                                                                                                                           >  1:805350/1‑150 (MQ=255)    aGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCgaga                                                                                                                                           >  1:794837/1‑150 (MQ=255)    aGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGaa                                                                                                                                          <  1:189785/151‑1 (MQ=255)     ggCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCCCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAg                                                                                                                                         <  1:500032/151‑1 (MQ=255)         cTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATc                                                                                                                                     >  2:668216/1‑151 (MQ=255)           gTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATc                                                                                                                                     <  2:903818/149‑1 (MQ=255)           gTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCtt                                                                                                                                   >  2:780188/1‑151 (MQ=255)             aTCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTa                                                                                                                                 <  1:264281/151‑1 (MQ=255)               cccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTAt                                                                                                                                <  2:151641/150‑1 (MQ=255)                ccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAGCCGTAAATAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTc                                                                                                                              <  2:729937/151‑1 (MQ=255)                     cGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGAttt                                                                                                                         >  1:520489/1‑151 (MQ=255)                      ggtATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGAttt                                                                                                                         >  2:589897/3‑150 (MQ=255)                       ttATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTAc                                                                                                                       >  1:906017/1‑151 (MQ=255)                         aTCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACt                                                                                                                      <  1:288043/150‑1 (MQ=255)                                aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAAt                                                                                                               <  1:638044/150‑1 (MQ=255)                                    aaaGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGc                                                                                                          <  1:563202/151‑1 (MQ=255)                                      agaagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTc                                                                                                        <  1:167274/151‑1 (MQ=255)                                           tggtgAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCaa                                                                                                      >  1:945239/1‑148 (MQ=255)                                            ggtgAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAaga                                                                                                   <  2:1643/150‑1 (MQ=255)                                               gAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagcag                                                                                               <  1:464869/151‑1 (MQ=255)                                                aaaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagcagc                                                                                              >  2:664263/1‑151 (MQ=255)                                                 aaGCCGCGCGCGGGGACCCTAAAGAGATGATCCAGAAGATGCTTTCTGAAATCGGAGGGCGTTTGCGGGAAGCGGGAATACCGTGCCGCGTCAGTGGAACCGAGACGCATCTTTATTAGATGTACTGCAAACTGGGGCCCAACGagcagc                                                                                              >  2:344778/1‑150 (MQ=255)                                                   gCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGTGCGTTTGCAGGAAGCGGGAATATCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATGGGATTTACTGCAAAAGGGTGCTCAAAGagctg                                                                                               <  2:164093/147‑3 (MQ=255)                                                   gCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGtt                                                                                           <  2:837918/151‑1 (MQ=255)                                                      gcgcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTc                                                                                        >  2:548098/1‑151 (MQ=255)                                                      gcgcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTc                                                                                        >  1:760117/1‑151 (MQ=255)                                                         cgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACt                                                                                     <  1:94974/151‑1 (MQ=255)                                                            ggCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGa                                                                                  <  1:679420/151‑1 (MQ=255)                                                             gCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGAt                                                                                 <  2:583245/151‑1 (MQ=255)                                                              cAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATc                                                                                <  2:165436/151‑1 (MQ=255)                                                                acggTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAt                                                                              <  1:202388/148‑1 (MQ=255)                                                                 ccGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAt                                                                              <  1:676821/150‑1 (MQ=255)                                                                 ccGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATg                                                                             >  2:3512/1‑151 (MQ=255)                                                                  cGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATgg                                                                            >  2:751155/1‑151 (MQ=255)                                                                   gTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATgg                                                                            <  1:632429/150‑1 (MQ=255)                                                                        gagaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCt                                                                      <  1:886881/151‑1 (MQ=255)                                                                           atgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACg                                                                   <  1:547663/151‑1 (MQ=255)                                                                            tgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTAc                                                                    >  2:516144/1‑149 (MQ=255)                                                                                    aagaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCcgcg                                                          <  2:676927/151‑1 (MQ=255)                                                                                       cttctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTg                                                        >  2:20316/2‑150 (MQ=255)                                                                                       aTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGa                                                       <  1:199380/151‑1 (MQ=255)                                                                                       aTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGa                                                       <  2:447288/151‑1 (MQ=255)                                                                                            ttctGAAATCGGAGGGCGTTTGCAGGAAGCGGGAAGACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTGCTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACGCGATCATGTACATCTACGCTTTCCGCGTGATc                                                     <  2:178950/148‑1 (MQ=255)                                                                                                   aTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGAtt                                           >  2:103246/1‑151 (MQ=255)                                                                                                     cGAAGGGCGTTTGCAGGAAGCGGGAATACCGCGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGAACATCTACGCTTTCGGCGTGATCGTCAATGATTc                                          >  1:719250/1‑150 (MQ=255)                                                                                                      gAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTg                                        >  2:726817/1‑151 (MQ=255)                                                                                                             gTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTg                                 >  1:640608/1‑151 (MQ=255)                                                                                                                tGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGtt                               >  1:334038/1‑150 (MQ=255)                                                                                                                 gCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTa                              <  1:140103/150‑1 (MQ=255)                                                                                                                  cAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATc                            <  1:379436/151‑1 (MQ=255)                                                                                                                  cAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATc                            >  1:314781/1‑151 (MQ=255)                                                                                                                    ggAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATcgc                          <  1:282845/151‑1 (MQ=255)                                                                                                                     gAAGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATcgcg                         <  1:742644/151‑1 (MQ=255)                                                                                                                       aGCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTg                       >  2:262476/1‑151 (MQ=255)                                                                                                                        gCGGGAATACCGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGc                      >  1:658599/1‑151 (MQ=255)                                                                                                                                  cGTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGAt            <  1:790363/151‑1 (MQ=255)                                                                                                                                   gTGCCGCGTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATg           <  1:759633/151‑1 (MQ=255)                                                                                                                                       cgcgTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGcac        >  2:336836/1‑150 (MQ=255)                                                                                                                                          gTCAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGcc    >  1:266347/1‑151 (MQ=255)                                                                                                                                            cAGTGGTAGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTg  <  1:623359/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                TTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTG  >  NC_000913/3822959‑3823248

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