Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,629,476	T→C	G349G (GGA→GGG)	ettA ←	energy‑dependent translational throttle protein EttA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,629,476	0	T	C	100.0%	142.3 / NA	45	G349G (GGA→GGG)	ettA	energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/16);  total (29/16)

AAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCGGACCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGTTCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACCTTGACGTACCCATTCCAGCTCTTTCTCAAT  >  NC_000913/4629230‑4629718                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   aaagCGGCCAACGAAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGGTTTGCTGTTATCAATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGGGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCtt                                                                                                                                                                                                                                                >  1:58307/1‑251 (MQ=255)      ggCCAACGTAGGCGCGGCTTGGCATCTCGGAGTTGCCGATCTTCATAATATACAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTACAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGGTTTACAAGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCttaccgg                                                                                                                                                                                                                                           >  1:7805/1‑246 (MQ=255)                 gcgcGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCt                                                                                                                                                                                                                                <  2:432432/251‑1 (MQ=255)                       gttggcgtatcGGTGTTTCCGAGCTTCAAAACATCCGCCCCGCCGGAATCTTCTGCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAGCTGGTCAACCGATGCCAGTTTCACCGTTTCGCCTAAAGGGATGGTGCCGCTGTCCGGCTGTTCCTGCCCAGAGATCATACGGAACAGGGTCGATGTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGtc                                                                                                                                                                                                                          <  2:225078/242‑1 (MQ=255)                         tGGCATCTCGGTGTTGCCGATCTTCTTAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTGATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGtc                                                                                                                                                                                                                          <  2:336928/249‑1 (MQ=255)                         tGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGtc                                                                                                                                                                                                                          <  1:379551/249‑1 (MQ=255)                            cATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGtcatca                                                                                                                                                                                                                      >  1:38723/1‑250 (MQ=255)                                  ggTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATcagca                                                                                                                                                                                                               <  1:120550/251‑1 (MQ=255)                                     gtTGCCGATCTTCATAATAGCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGAGGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGGCGATTTACCCTCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATcagcagg                                                                                                                                                                                                             <  1:453668/250‑2 (MQ=255)                                                                   ggAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTg                                                                                                                                                                               >  1:33791/1‑250 (MQ=255)                                                                                                tATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGTTTGATGCCTTCCAGCACTTTATCGCACAGAcgcgtt                                                                                                                                                 >  1:480680/1‑249 (MQ=255)                                                                                                  tCCATTGAGTCACGGAACTGCTCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGCTGATCCCGACGATTGCCCCTTTCGGGATCGCGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTCGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGCCGCGGTcc                                                                                                                                               >  2:45226/1‑251 (MQ=255)                                                                                                   ccATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTcc                                                                                                                                               >  1:505314/1‑250 (MQ=255)                                                                                                    cATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTcc                                                                                                                                               >  1:92282/1‑249 (MQ=255)                                                                                                        gAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTgg                                                                                                                                         >  1:216284/1‑251 (MQ=255)                                                                                                               ggAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGGTGCTGACTTCCAGCACTTTCACCCCCCGACGCGGTCCAGGGGGACTAAAc                                                                                                                                  >  2:348897/1‑251 (MQ=255)                                                                                                                 aaCTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAAc                                                                                                                                  >  1:377956/1‑249 (MQ=255)                                                                                                                      aTCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGt                                                                                                                           >  1:500771/1‑251 (MQ=255)                                                                                                                           caggtgCCAGTTACACCGGTTCACCTAAAGTGTAGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTGCGGACCGATGATCCCGACGACTGCCCCTTGCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCCCCGTAGGATTTACGCAGATTGCTGACTTCCAGCATTTTATCGCCCATACGCGGTCCAGGTGGAATAAACAGTTCGTTGGtt                                                                                                                      <  2:161831/246‑1 (MQ=255)                                                                                                                              gggcaattttacccgttgcaccAAAAGGGATGGTGCCGATGTCCGGCTGCTCCTGACCAGATATCATTCGGAACAGGGACGATTTACCCGCGCCGTGCGGACCGAGGATCCCGACGAGTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCTGTCCAGGTGGAATAAACAGTTCGTTGGttgcg                                                                                                                   <  2:189775/239‑4 (MQ=255)                                                                                                                                          ccgtttcaccTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTg                                                                                                        >  1:72305/1‑250 (MQ=255)                                                                                                                                            gtttcaccTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTc                                                                                                                             <  1:114245/227‑1 (MQ=255)                                                                                                                                            gtttcaccTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTc                                                                                                                             >  2:114245/1‑227 (MQ=255)                                                                                                                                             tttcaccTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTT‑ATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATAtt                                                                                                   >  1:207051/1‑251 (MQ=255)                                                                                                                                                   cgaacgTGATGGTGCCGCTGTCCGGCTGTTACTGACCAGAGATCATACGGAACAGGGTCGGTTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCCATAAGCAGGCGATCACCGTAGGATTTACGCCGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGAGGTCCAGGTGGAATCAACAGCTCGTTTGTTTCGTTAAGTTTATGGTATTCGGTg                                                                                              >  2:261432/6‑251 (MQ=255)                                                                                                                                                       aGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGt                                                                                          >  1:183884/1‑251 (MQ=255)                                                                                                                                                           aTGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGTTCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCGAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTc                                                                                        >  2:64163/1‑249 (MQ=255)                                                                                                                                                                      gTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCaaa                                                                             <  1:398133/249‑1 (MQ=255)                                                                                                                                                                           gCTGGTCCTGACCAGAGATCAGACGGAACAGGGTCGATTTACCCGCGCCGTGCGGACCGATGATCCCGACGAGTGCCCCGTGCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGTTTGCTGACTTCCAGCACTGTATCGCCCAGACGCGGGCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCGCTTCAAAGCGTGcc                                                                      <  2:428842/251‑1 (MQ=255)                                                                                                                                                                             tGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGcc                                                                      >  1:410274/1‑249 (MQ=255)                                                                                                                                                                             tGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTAACCGCGCCGTTCGGACAGATGATCCCGACGATTGCCCCTTTCGGGCACGAGAAGCTCAGGTCATCAATCAGCAGGCGCTCACCGTAGGCTTTACGCAGGTCGCTGCCTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTAGGTTTATGGTATTCGGTGATGTTCAGCTCTTCACAGCGGGcccg                                                                    >  2:504792/1‑249 (MQ=255)                                                                                                                                                                                        agagATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCACAGCGTGCCAGACGTGCTTTAc                                                         >  2:10383/1‑251 (MQ=255)                                                                                                                                                                                        agagATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTAc                                                         >  2:232769/1‑251 (MQ=255)                                                                                                                                                                                           gaTCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCtt                                                      >  2:97315/1‑251 (MQ=255)                                                                                                                                                                                           gaTCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCtt                                                      >  2:137900/1‑251 (MQ=255)                                                                                                                                                                                              cATACGGACCAGGGTCGATTTACCCGCGCAGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGAGTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTGCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCGTTc                                                    <  2:38723/250‑1 (MQ=255)                                                                                                                                                                                               aTACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAGGCGTGCCAGACGTGCTTTACCTTTCg                                                   >  2:348287/1‑250 (MQ=255)                                                                                                                                                                                                tACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCCGGATCGAGCAGCTCAGGTCATCAATCAGCAGGCGATCAACGTAGAATTTACACAGGTTACTGACTTCCAGCACTTTATCGCCCAGACGCGGACCAGGTGGACTAAAAAGATCGTTGGATTCGTTCCGTTTCTGATATTCGGTGATGTTCAGCTCTTCAAAGCGTGCACGACGTGCTTAACCTTTCGAc                                                 >  1:58133/1‑251 (MQ=255)                                                                                                                                                                                                     aaCAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTggcg                                            >  1:294896/1‑248 (MQ=255)                                                                                                                                                                                                                     ccGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAAGCAGCAGGCGATCCCCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGtgccctga                            <  2:72305/251‑8 (MQ=255)                                                                                                                                                                                                                       gcgccGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACg                          <  1:121775/248‑1 (MQ=255)                                                                                                                                                                                                                                     gattatCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGc             <  1:52775/246‑1 (MQ=255)                                                                                                                                                                                                                                     gatgatCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCCATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGCATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCCGCCGTGCTTTACCTTTCGCCTGGCGGCCTTTAGTTCCCTGACGTACCCCCTCCAGcc            >  2:52775/1‑250 (MQ=255)                                                                                                                                                                                                                                            ccGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTtctc     >  1:420446/1‑251 (MQ=255)                                                                                                                                                                                                                                               acgaTTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGATCACCGTAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGCTCTTCAAAGCGTGCCAGACGTGCTTTACCTTTCGACTGGCGGCCTTTAGTTCCCTGACGTACCCACTCCAGCTCTTtctcgat  <  1:137900/251‑4 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   AAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCGGACCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGATCGCCATAGGATTTACGCAGGTTGCTGACTTCCAGCACTTTATCGCCCAGACGCGGTCCAGGTGGAATAAACAGTTCGTTGGTTTCGTTACGTTTCTGATATTCGGTGCTGTTCAGTTCTTCAAAGCGCGCCAGACGTGCTTTACCTTTCGACTGACGGCCTTTAGTACCTTGACGTACCCATTCCAGCTCTTTCTCAAT  >  NC_000913/4629230‑4629718

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