Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,387,978	C→T	A465T (GCC→ACC)	mscM ←	miniconductance mechanosensitive channel MscM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,387,978	0	C	T	100.0%	125.2 / NA	38	A465T (GCC→ACC)	mscM	miniconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/16);  total (22/16)

CGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCG  >  NC_000913/4387733‑4388221                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    cGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGc                                                                                                                                                                                                                                                >  2:5046/1‑251 (MQ=255)       gcgcAAGCCGTAACCCAGCGTCCTCCACAGCACAGGGAACGGTGACGCGACGAGCATCTACCAGAAAAGAGTGCGCACCGTCAGCCAGAAGTGATCCTGAGTCACGTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGTCGCGAGTAAATTCTGCAGCCGGCCAGAATCAACGCTCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCAGCACACTGGTTTTGTCCAact                                                                                                                                                                                                                                          <  2:161489/251‑1 (MQ=255)          cAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGACAGCATCATCACACTGGTTTTGCCCAactgac                                                                                                                                                                                                                                       >  1:137222/1‑251 (MQ=255)              ccGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAActgactga                                                                                                                                                                                                                                    >  2:242669/1‑250 (MQ=255)                        gCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGAGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGc                                                                                                                                                                                                                         <  2:172641/251‑1 (MQ=255)                                   aGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGGGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTGCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAgacg                                                                                                                                                                                                              <  2:83437/251‑1 (MQ=255)                                                gTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTgg                                                                                                                                                                                                 >  1:38702/1‑251 (MQ=255)                                                           cgaATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTcc                                                                                                                                                                                        >  2:67113/2‑249 (MQ=255)                                                                cGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGc                                                                                                                                                                                  <  1:242062/250‑1 (MQ=255)                                                                              tGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATc                                                                                                                                                                   >  1:33138/1‑251 (MQ=255)                                                                                         aGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTc                                                                                                                                                        >  2:77685/1‑251 (MQ=255)                                                                                          gCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTc                                                                                                                                                        >  2:233323/1‑250 (MQ=255)                                                                                           ccAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCCCGTc                                                                                                                                                        >  2:143316/1‑249 (MQ=255)                                                                                             aGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAACTAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTTTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGGGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGCTTTCCAGCGGCCAGGCGATGGTAATCGGGCGCACGTccg                                                                                                                                                      >  2:172191/1‑247 (MQ=255)                                                                                                         gAGTCACGTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAg                                                                                                                                         <  1:22777/250‑1 (MQ=255)                                                                                                             cACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGc                                                                                                                                    >  2:49919/1‑251 (MQ=255)                                                                                                                  tGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGc                                                                                                                               >  2:148589/1‑251 (MQ=255)                                                                                                                   gCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCACCCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGc                                                                                                                               <  1:118832/250‑1 (MQ=255)                                                                                                                          ggcgccgcCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGAGCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTTATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTc                                                                                                                       <  2:37861/249‑1 (MQ=255)                                                                                                                           tCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTTCGTTGCTTc                                                                                                                       <  2:37847/250‑1 (MQ=255)                                                                                                                                 ccgAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCAcc                                                                                                                >  1:65983/1‑251 (MQ=255)                                                                                                                                   gAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACctc                                                                                                              <  1:70041/251‑1 (MQ=255)                                                                                                                                                      tGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGt                                                                                           >  1:104400/1‑251 (MQ=255)                                                                                                                                                           tAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGcc                                                                                       <  1:112356/250‑1 (MQ=255)                                                                                                                                                           tAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGcc                                                                                       >  2:112356/1‑250 (MQ=255)                                                                                                                                                                   cgAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGgaga                                                                              <  1:233323/251‑1 (MQ=255)                                                                                                                                                                    gAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGAc                                                                             >  2:85193/1‑251 (MQ=255)                                                                                                                                                                               gCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTgg                                                                  <  1:5046/251‑1 (MQ=255)                                                                                                                                                                                  aCCGACCAGAATCAACGCTCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCATATCCTGGGCGATTTCCAGCGGCCAGGCTATGGTCATCTGGCTCACGTCAGAGGTCCAGAACAGATAGCGTTGCGTTGCTTCGTTCGCCTCTTTCAGCGCATCCTGCAGTTGCCCGTTTGAGACTTTCAGCTTGGTCa                                                               <  1:144481/251‑1 (MQ=255)                                                                                                                                                                                        cAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTcc                                                          >  1:222378/1‑250 (MQ=255)                                                                                                                                                                                              cAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAAtagt                                                    >  1:220091/1‑250 (MQ=255)                                                                                                                                                                                                     ccAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTc                                            >  1:228340/1‑251 (MQ=255)                                                                                                                                                                                                     ccAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTc                                            <  2:75111/251‑1 (MQ=255)                                                                                                                                                                                                     ccAAACAGCGGCAAAATTGTCTCTTTACTGGTCACCATCATCACACTGGTTTTGCCCAACTGACTGAAGGAGTCCAGCGAAATGAGACGACGCAGATCCTGAGCGATTTCCAGCGGCCAGGCGATGGTCATCTAGCGCACGTACGAGGTCAAGACCAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGAATCCTCCAGTTGCCCGTTGGAGACTTTACGCTTGGTCAGCTACAATAGAAGCGCGTc                                            >  1:228330/1‑251 (MQ=255)                                                                                                                                                                                                        aaCAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGGCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGATACTTTCAGCTTGGTCAGTTCCAATAGTAGCgttgcgcc                                         >  2:234569/1‑245 (MQ=255)                                                                                                                                                                                                                             ttGCTGGTCAGCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTc                      <  1:162340/249‑1 (MQ=255)                                                                                                                                                                                                                                       gCATCATCACACTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGc          >  2:41002/1‑251 (MQ=255)                                                                                                                                                                                                                                               acacTGGTTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTgcg  <  1:148589/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    CGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAACGCGCCAAACAGCGGCAAAATTGTCTCTTTGCTGGTCAGCATCATCACACTGGCTTTGCCCAACTGACTGAAGGTGTCCAGCGAAATGAGACGACGCAGATCCTGGGCGATTTCCAGCGGCCAGGCGATGGTCATCGGGCGCACGTCAGAGGTCCAGAACAGATAGCGGTGCGTTGCTTCGTTCACCTCTTTCAGCGCATCCTCCAGTTGCCCGTTGGAGACTTTCAGCTTGGTCAGTTCCAATAGTAGCGTGTCGCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCG  >  NC_000913/4387733‑4388221

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