Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,371,696	C→T	D224D (GAC→GAT)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,371,696	0	C	T	100.0%	160.5 / NA	47	D224D (GAC→GAT)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/28);  total (19/28)

ATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGACAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATC  >  NC_000913/4371454‑4371939                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    aTTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAAt                                                                                                                                                                                                                                             >  1:74680/1‑251 (MQ=255)   tGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTGCATCCTGCTGGCTGATAAGAAAAtc                                                                                                                                                                                                                                            <  2:256552/250‑1 (MQ=255)       cAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTTATCGCTGAAGCGAGGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGGCTCCTTACTTCAGCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTcc                                                                                                                                                                                                                                         <  2:310697/249‑1 (MQ=255)                      tcCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAAtgct                                                                                                                                                                                                                        >  2:27268/1‑251 (MQ=255)                       ccGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAAtgct                                                                                                                                                                                                                        <  1:258745/250‑1 (MQ=255)                       ccGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAAtg                                                                                                                                                                                                                          <  1:108957/248‑1 (MQ=255)                                    cgaAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAACTCTCCAACATCCGCGAAATGCTGCCGGTTCTGGa                                                                                                                                                                                                            >  2:223999/1‑249 (MQ=255)                                       aaCCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTg                                                                                                                                                                                                       >  1:73953/1‑251 (MQ=255)                                                      gATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGc                                                                                                                                                                                         >  1:132457/1‑250 (MQ=255)                                                              aaGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCgctg                                                                                                                                                                                <  2:129835/251‑1 (MQ=255)                                                                gCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCgctgct                                                                                                                                                                              >  2:334391/1‑251 (MQ=255)                                                                   aTGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGTAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGACATGCTGCCgg                                                                                                                                                                                                                   >  2:202751/1‑211 (MQ=255)                                                                   aTGGACAAAGGCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGACTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCGTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCgg                                                                                                                                                                                                                   <  1:202751/211‑1 (MQ=255)                                                                      gACAAAGTCGGGAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTCGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGCTCTGGAAGCTGTTGCCAACGCAGCAAAACCGCTGCTGatcat                                                                                                                                                                        >  2:326895/1‑251 (MQ=255)                                                                       aCAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGatcatc                                                                                                                                                                       <  2:270312/251‑1 (MQ=255)                                                                              cGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGaa                                                                                                                                                                 <  1:211055/250‑1 (MQ=255)                                                                                   aagaagGAGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCGGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGAGCATCGCTGAAGATGTa                                                                                                                                                           <  2:151645/251‑1 (MQ=255)                                                                                      aagGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGGGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGTAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCTCGAAATGCTGCCGGTTCTTGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGaa                                                                                                                                                        <  1:378332/251‑1 (MQ=255)                                                                                       agGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAg                                                                                                                                                       <  2:367051/251‑1 (MQ=255)                                                                                          cGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCg                                                                                                                                                    >  2:314165/1‑251 (MQ=255)                                                                                             tATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGaa                                                                                                                                                  <  1:314165/250‑1 (MQ=255)                                                                                              tgcACCGTGGAAGACGGTACCGGTCTTCAGGACCAACTGGACGTGGTTGAGGGTATGCAGTTCGACCGTGGCTCCCTGTCTCCTTTCTTCATCAGCAAGCCGGAAACTGCCGCAGTAGAACTGGAAAGCCCGTTCAGCCTGCTGGCTTATAAGAAAATCTCCAACATCCGCTAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGATGATGTATAAGGCGAAg                                                                                                                                                 <  2:201753/248‑1 (MQ=255)                                                                                                 aCCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCt                                                                                                                                             >  2:388200/1‑251 (MQ=255)                                                                                                                aCCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGgttg                                                                                                                               <  2:157144/250‑1 (MQ=255)                                                                                                                           ggACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGc                                                                                                                   <  2:195070/251‑1 (MQ=255)                                                                                                                                   tGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATc                                                                                                           >  2:395328/1‑251 (MQ=255)                                                                                                                                    gtgcGTGGTTGAAGGTATGGAGTTCGACCGTGGCTACCTGTCTCCTTACTTCAGCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGACCGTTCATCCTGCTGGCTTATAAGAAAAGCTCCAACAGCCGCTAAAGGCCGCCGTTTCTGGAAGCTGTTGCCAAAGCAGGTAAATCGCTGCTGATCATCGCTGATGATGTAGAAGGCGAAGCGCTGGCAACTCTTGTTGTTAACACCATGCGTGGCAt                                                                                                            <  1:225405/246‑1 (MQ=255)                                                                                                                                            ttGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTc                                                                                                  >  1:184993/1‑251 (MQ=255)                                                                                                                                             tGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTc                                                                                                  <  2:117429/250‑1 (MQ=255)                                                                                                                                             tGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCg                                                                                                 >  1:187242/1‑251 (MQ=255)                                                                                                                                              gAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGc                                                                                                <  1:182221/251‑1 (MQ=255)                                                                                                                                                     ggCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCTTTCATCCGGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCgttt                                                                                         <  1:250440/250‑4 (MQ=255)                                                                                                                                                      gCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTa                                                                                        <  2:132457/251‑1 (MQ=255)                                                                                                                                                              aCCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCg                                                                                >  2:358305/1‑251 (MQ=255)                                                                                                                                                                  tGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGCAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGCTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGCTGTTAACACCATGCGTGGCATCGTGCAAGTCGCTGCGGTTAAAGCACCGGGCt                                                                            >  1:192794/1‑251 (MQ=255)                                                                                                                                                                           gTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGAtc                                                                   <  2:317734/251‑1 (MQ=255)                                                                                                                                                                           gTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGAt                                                                    <  2:61457/250‑1 (MQ=255)                                                                                                                                                                                      tcatcaACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGc                                                         <  1:358305/250‑1 (MQ=255)                                                                                                                                                                                       catcaACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCt                                                        <  1:388200/250‑1 (MQ=255)                                                                                                                                                                                                     aaaCTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAAGCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATc                                         <  2:26019/251‑1 (MQ=255)                                                                                                                                                                                                               aGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCTGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTg                                >  2:130296/1‑250 (MQ=255)                                                                                                                                                                                                               aGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTAtg                                                     <  2:159068/229‑1 (MQ=255)                                                                                                                                                                                                               aGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTAtg                                                     >  1:159068/1‑229 (MQ=255)                                                                                                                                                                                                                         taaaGCCCGTTCATCCTGCTGGCTGGTACGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTACAAGGCGAAGCGCTGGCAACTCTGGTTGTTCACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCGGCAGAATATCGCAACCCTGACGGGCGGTAc                     >  2:85455/2‑251 (MQ=255)                                                                                                                                                                                                                               ccGNTAATCCTGCTGGCTGATAATAAAAGCTCCAACATCCGCGAAATTCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGCTCATCGATGAAGATGTAGAAGGTGAAGCGCTGGCAACTCTGGTTGTTAACGCCATGCGTGGCATCGTGAAAGTCGCGGCGGTTAAAGCAACGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGAt               <  2:375727/251‑1 (MQ=255)                                                                                                                                                                                                                                            tggctgATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAACGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCCTGATCCCTGAag       >  1:176106/1‑246 (MQ=255)                                                                                                                                                                                                                                             ggctgATAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATc  <  1:198146/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    ATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGACAAGAAAATCTCCAACATCCGCGAAATGCTGCCGGTTCTGGAAGCTGTTGCCAAAGCAGGCAAACCGCTGCTGATCATCGCTGAAGATGTAGAAGGCGAAGCGCTGGCAACTCTGGTTGTTAACACCATGCGTGGCATCGTGAAAGTCGCTGCGGTTAAAGCACCGGGCTTCGGCGATCGTCGTAAAGCTATGCTGCAGGATATCGCAACCCTGACTGGCGGTACCGTGATCTCTGAAGAGATC  >  NC_000913/4371454‑4371939

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