Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,867,213	A→C	G113G (GGT→GGG)	rpoS ←	RNA polymerase sigma factor RpoS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,867,213	0	A	C	100.0%	125.7 / NA	39	G113G (GGT→GGG)	rpoS	RNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (21/18);  total (21/18)

CGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTG‑GCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCT  >  NC_000913/2866967‑2867459                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         cGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCAc                                                                                                                                                                                                                                                      <  1:361795/251‑1 (MQ=255)    aCTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGAt                                                                                                                                                                                                                                                   >  2:291444/1‑251 (MQ=255)     cTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACgctt                                                                                                                                                                                                                                                  >  2:237145/1‑248 (MQ=255)                 tCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAAcggcg                                                                                                                                                                                                                                      >  2:62535/1‑251 (MQ=255)                           gggACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAAttttt                                                                                                                                                                                                                            >  1:127125/1‑251 (MQ=255)                                           ggTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt                                                                                                                                                                                                            <  1:25660/251‑1 (MQ=255)                                              tCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCCGACGCAAGTTAc                                                                                                                                                                                                         >  1:258336/1‑251 (MQ=255)                                               cGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATATGAAGGGGTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGTGGGCGGATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTGTTACCACCAGaagaaagggact                                                                                                                                                                                                        <  1:306329/251‑12 (MQ=255)                                                          gTTCAGCTCCTTTACGATGTGAAGCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCTTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATcc                                                                                                                                                                                             <  2:189172/251‑1 (MQ=255)                                                                      tACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCCATTTTTACCACCCGGCGCAAGTTACTCTCGCTCATCCggcggccggacg                                                                                                                                                                                 >  1:306866/1‑245 (MQ=255)                                                                        cGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCg                                                                                                                                                                               <  1:286460/251‑1 (MQ=255)                                                                              gAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACAtc                                                                                                                                                                          <  2:337769/250‑1 (MQ=255)                                                                                  cGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAAATTTTTCCACCCGACGCAAGGTACTCCCGATCCTCCCGCGGCGAGAGGCGGCAACTCCAc                                                                                                                                                                     >  1:347442/1‑251 (MQ=255)                                                                                                  cGGGTTTGGTTCATAATCGCCCGTTCAAGCGTCTGGCGAATCCACCAGGTGGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTCCCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACAtc                                                                                                                                                                          <  1:182816/230‑1 (MQ=255)                                                                                                     gTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAAATTTTTCCACCCGACGCAAGTTCATCTCGATCATCCGGCGGCGGGAGGGCACATCTCCCCGCAGTGGGCGACGGG‑Cac                                                                                                                                                  >  1:278946/1‑249 (MQ=255)                                                                                                         ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTTCCACCCGACGCCAGGTACTCTCGCTCCACCGGCGGCGGGAGGCGGCCTCTcccccgatcggaagagcgtccgggtgg                                                                                                                                              >  2:71423/1‑226 (MQ=255)                                                                                                         ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCAc                                                                                                                                                                     <  1:71423/228‑1 (MQ=255)                                                                                                                aaTCGCCCGTTCAATAGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGGAATTTTTAACACCAGACGCAAGGTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCACGACGCG‑CAAAATAAActtctt                                                                                                                                      >  1:280575/1‑251 (MQ=255)                                                                                                                       cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTCTTCGGc                                                                                                                               <  1:19346/251‑1 (MQ=255)                                                                                                                                    ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAAATTTTTCCCCCCGACGCAAGATAATCTCCCTTCTCCCGCGGGGGGGGGCGACATCTCCACGCCGTGCGCCACGCG‑CCAAATAAAATTCTTTTTCGTCCCTTAACaagggg                                                                                                                  >  2:80647/1‑246 (MQ=255)                                                                                                                                           ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCCGACGCAAGTTACTCTCGATCATCcggcg                                                                                                                                                                                         >  1:330887/1‑174 (MQ=255)                                                                                                                                           ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcggcg                                                                                                                                                                                         <  2:330887/174‑1 (MQ=255)                                                                                                                                                            ttGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCTAGTTTCTCTCTATCATCCGGCGGAGAGGGGCGACATCTCCACGCCGAGCGCGACGCG‑CACAATAAACTTCTTCTTCGGCCCTTAACAGAGGTGAATAAACAATCTCACCAAGGTa                                                                                           >  1:52305/1‑250 (MQ=255)                                                                                                                                                             tGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGACGCCAGTTAATCTCGATACTCCCGCGGCGGGAGGCGACACCTCCCCCCCGGGCGCGGCCCGGCAAAA‑AAACCTCTTCTTCGGCCGCTTACAGTGGTGAATAACCacacccacacaggtaa                                                                                          >  2:240606/1‑235 (MQ=255)                                                                                                                                                                   gCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGGGCGCGACGCG‑CAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGGATAACCAATCTCACCAAGGTAAAACTg                                                                                     >  1:108964/1‑249 (MQ=255)                                                                                                                                                                            aCGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCa                                                                          <  2:306866/251‑1 (MQ=255)                                                                                                                                                                              gTTCCGGGTCAAAATTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCAACGATTGCCATAACGGCGGGCAATTTTTTCCCCCAGAAGCAAGTTACTCTTCAAAATCCCGCGGCGAGAGGGGACATCTCCACGCAGGTCGCGACGCG‑CCAAATAAAATTCTTTTTTGGCCGGTAACAGTGGGGCATAACCAATCTCCCCAAGGTAAAGCTGAGTTGCGTCCAac                                                                        >  1:373393/1‑251 (MQ=255)                                                                                                                                                                                                ctACCGCGCGGAGCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTG‑GCt                                                        <  2:328419/248‑1 (MQ=255)                                                                                                                                                                                                                 cccAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCAACCGACGCCAGGTACTCTCGATCATCCGGCGGCGAGAGGCGACAACTCCACGCAGTGCGCGACGCG‑CCAAATAAACTTCTTCTTTGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGGCCCACACACGCTGTGTG‑GGCCCCTGCGATAAAAGTTActc                                    >  1:219817/1‑251 (MQ=255)                                                                                                                                                                                                                     ggatgaccacatCGATAAGGTCCAGCAGCGCCAGCACACGATAGCCATAACGGTGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCAGCGCCGAGAGGCTACCTCCCCCCGCAGTGCGCGCCGCG‑CAAAATAAACTTCTTCCTCGGCCGTTAACAGTGCTGAATAACCAATCTCCCCAAGATAAAGCTGAGTCAAGTCCAACACTCGCGCTGTG‑GCTCCCTGCGATAAAAGTCCCTCTTca                                <  2:223110/240‑2 (MQ=255)                                                                                                                                                                                                                     ggTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTTTCACCCGACGCCAGTTACTCTCGTTCCTCCGGCGGCGGGAGGGGACATCTCCACGCCGGGCGCGACCCG‑GAAAATAACATTTTTCTTCCGCCGTTAACAGTGGGGGATAAACAAACTCCCCCAAGAAAAGCTTAGTCGCGTCCAAAACCCGCTCT‑TGTGCTCCCTGCCATCACAGTTtctcctct                                >  2:310614/1‑247 (MQ=255)                                                                                                                                                                                                                              ctTCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTCTGCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTG‑GCTCCCTGCGATAACAGTTCCTCTGCGGCCAAAt                         <  2:219817/249‑1 (MQ=255)                                                                                                                                                                                                                                tCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTG‑GCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCg                       <  1:41473/249‑1 (MQ=255)                                                                                                                                                                                                                                tCGATAAGGTCCAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTG‑GCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCg                       <  1:62535/249‑1 (MQ=255)                                                                                                                                                                                                                                       ggTCCATCAACGTCAGCCCACGATTGCCATAACGGCCGGCAATTTTTCCCCCCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTTGTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGGCCAACACACGCTTTGTG‑GCTCCCTGCGATAACAGTCCCTCTTCGGCCAAATTGTTATCACTg              <  2:174440/251‑1 (MQ=255)                                                                                                                                                                                                                                          ccAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGGGCGCGACGCG‑CAAAATAAAATTCTTCTTCGGGCGTTAACAGTGGGGAATAAACAAACTCAACAAGGTAAAGCTGAGTCGCGTCCCACACACGCTGTGTG‑GGCCCCTGCGAAAACAGTTCCTCTTCGGCCAAATCGTTATCACTgggg           >  2:26373/1‑250 (MQ=255)                                                                                                                                                                                                                                          ccAGCAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAAACGCAAGTTACTCTTGATCATCCGGCGGCGAGAGGCGACATCTCCCCGCAGTGGGCGACGCG‑CAAAATAAAATTCTTCTTCGGCCGTTAAAAGTGGTGAATAACCAATCTCACCAAGGGAAAGCTGAGTCGCGTCCAACAAAAGCTGTGTG‑GGTCCCTGCGATAACAGGTCCTCTTCGGCCAAATCGGTATCACTGGGt           >  1:271468/1‑251 (MQ=255)                                                                                                                                                                                                                                              cAACGCCAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTG‑GCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCt       <  2:355671/251‑1 (MQ=255)                                                                                                                                                                                                                                                   ccAGCCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGGAAAGCTGAGTCGCGTCCAACACACGCTGTGTG‑GCTCCCTGCGATAAAAGTTCCTCTTCGGCCAAATCGTTATCCCTGGGTTCCTGTTct  >  1:86176/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         CGCACTTGGTTCATGGTCCAGCTTATGGGACAACTCACGTGCGGTTCGCAGGTAAACGTTCAGCTCCTTTACGATGTGAATCGGCAAACGAATAGTACGGGTTTGGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCG‑CAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTG‑GCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCT  >  NC_000913/2866967‑2867459

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