Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,867,455	G→A	Q33* (CAG→TAG)	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,455	0	G	A	100.0%	136.5 / NA	44	Q33* (CAG→TAG)	rpoS	RNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/26);  total (18/26)

GCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTAC  >  NC_000913/2867222‑2867705                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           cccATAACGGGGGGCTATGTTTACCTCCAGTCGCTAGTGACCCTCGCTCAACCGGCGGCGAGAGGCGACCACTTCACGCAGTTCGCGATGCGCAAAATAATCTGCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACCTGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCtttt                                                                                                                                                                                                                                           <  2:158232/250‑1 (MQ=255)            gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAAcc                                                                                                                                                                                                                                <  1:194007/251‑1 (MQ=255)                aatgtttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCaa                                                                                                                                                                                                                            <  1:347194/247‑1 (MQ=255)                             aGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGGGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAGCACTTGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                                                                                                                                                                                                               <  1:436616/251‑1 (MQ=255)                               aCGCAAGTTACTCTCGATCAGCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaa                                                                                                                                                                                                             <  1:537677/251‑1 (MQ=255)                                      ttACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGc                                                                                                                                                                                                      >  1:507333/1‑251 (MQ=255)                                        aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                                                                                                                                                                    >  1:190514/1‑251 (MQ=255)                                               atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTAAATAACCAATCTCACCAAGGTAAAGCTGAGACGCGTCCAACACACGCTGTGAGGCTCCCTGCGCTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTAAACAACCTCAACTCCGTTCTCATCaacatcccaatcttaatt                                                                                                                                                                                             >  2:666852/1‑235 (MQ=255)                                                tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAttt                                                                                                                                                                                            >  1:611774/1‑251 (MQ=255)                                                 catcCGGCGGCGAGAGGCTCCGTCGCCACGCAGTGCGCGACGCGCAAAATAAACTTCTGCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa                                                                                                                                                                                           <  2:396272/251‑1 (MQ=255)                                                      ggcggcGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCCTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAGCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCa                                                                                                                                                                                      <  1:253636/251‑1 (MQ=255)                                                        cggcgAGAGGCGACAGCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACATTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                                                    <  1:413207/251‑1 (MQ=255)                                                          gcgAGAGGCGACATCTCCACGCCGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGGTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGaa                                                                                                                                                                                  <  2:546476/251‑1 (MQ=255)                                                           agagaGGCGACACCTCCACGCAGTGAGCGACGCCCAAAAAAATCTTCGTCTGCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAAc                                                                                                                                                                                 <  2:192211/250‑1 (MQ=255)                                                               aggagACATATCCACACAGTGCGCGACGCTCAAAACAAACTTCTTCTTCGTCCGTTAAAAGTGGTGAATAACCAATCTCACCCAGGTAAACCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACATTTCCTCTTCGGCCAAATCGGTATCACTTTGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTCCTCATCAAATTCCGCATCTTCATTTAAATCATGaatttgc                                                                                                                                                                             <  2:353854/247‑6 (MQ=255)                                                                   gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCg                                                                                                                                                                         <  1:648500/251‑1 (MQ=255)                                                                     aaactcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCATCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGTTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTATATCATGAGCTTTCAGCGTa                                                                                                                                                                       <  2:445063/248‑1 (MQ=255)                                                                                     gcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGt                                                                                                                                                       <  1:386986/251‑1 (MQ=255)                                                                                           cgcaaaataaacgacgtcgtcttccGTTAACAGTGCTGAATAACCAATGTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCGGTGTGGCTCCCTTCGATAACAGTTCCTCTTCGGGCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGCTCTCATCAAATTCCGCATTTTCATTTAAACCATGAACTTTCAGCGTATTCTCACTCATAAGGTGGCTcc                                                                                                                                                 <  2:300071/228‑1 (MQ=255)                                                                                                     aCTTCTTCTTCGGCCGTTAACAGTGGAGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTACGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCACATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGACTACTGACTCATCCGGTGGCTCCTACCCGTGAt                                                                                                                                       >  2:471135/1‑251 (MQ=255)                                                                                                       ttcttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATcc                                                                                                                                     >  1:395746/1‑251 (MQ=255)                                                                                                               cGGCCGTTAACAGGGGTGATTAACCAATCTAACCAAGGTAAAGCTGAGTCGAGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCACAAACCTCACCTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGAGTATTCTGACTCATAAGGCGGCACCTACCCTGGATCACTTGACgg                                                                                                                             >  2:89306/1‑251 (MQ=255)                                                                                                                 gCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGAAGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTTAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTAATCCCTTGACGGa                                                                                                                            >  1:538022/1‑250 (MQ=255)                                                                                                                   cGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAAc                                                                                                                          >  1:655161/1‑250 (MQ=255)                                                                                                                    gTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAt                                                                                                                        <  1:514825/251‑1 (MQ=255)                                                                                                                        aCAGTGGTGAATAACCAATCTCCCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCaa                                                                                                                    <  2:190514/251‑1 (MQ=255)                                                                                                                        aCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTACATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGCCGGAACATTCa                                                                                                                     >  2:114435/1‑250 (MQ=255)                                                                                                                                           cgcaacAAGGTATAGCGGAGTCGCGTCCAACACCCGTTTGGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCATATCTTTATCCCTGTGTTCCTATTCTACTAAGGCCTTTTCGTCAAAATCCTCAACTCCGTTCGCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTTAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGAGCCCTTGCCGGAACATTCAATCAAAAGCCTGGTTccgcc                                                                                                 <  2:517577/246‑1 (MQ=255)                                                                                                                                                            gtgtctcgTCCCACACTCGCTTTGTGGCCCCCGTCGAGAACAGTTCTTCTTCGGCCAAATCTTTATCACTGGGTTCATATTCTACTAAGGCCTTTTCGTCAAAACCCTCATCTCCTTTCTCACCAAATTCCGCCTTTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCTCTGCGGc                                                                                 <  2:577761/244‑1 (MQ=255)                                                                                                                                                                 gcgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATaa                                                                           <  1:597041/251‑1 (MQ=255)                                                                                                                                                                        acacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCt                                                                                                                                                   >  2:268103/1‑172 (MQ=255)                                                                                                                                                                        acacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCt                                                                                                                                                   <  1:268103/172‑1 (MQ=255)                                                                                                                                                                                           ttGCGATAACAGTTCCTCTTCGGCAACATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTAAACTCCGTTTTCATCAATTTCAGCATCTTAATTTAAATCATGAACTTTCAGCGTATTCTGACTCATACGGTGGCTCCTACCCGTGATCCATTGACGGAACATTCAATCAAAAGCATGGTTCCGCCGCTTTATCGCTGCGGCAAATAAATCAGCGGTCTTACGGATTTcccatt                                                 >  2:348390/2‑248 (MQ=255)                                                                                                                                                                                                aTAACATTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGAGTTCCCCTTGTAAc                                            <  2:507333/251‑1 (MQ=255)                                                                                                                                                                                                            actcagCCAAATCTTTATCACGGGGTTCCTATTCTACTAAGGCCTTTTCGTCCAAAAACTCAACACCGTTCTCAGCAAATTATGCATCTTCATTTAAATCATGAACTTTCAGCGAATTCTGACTCATAAGATGGCTCATACCCGTGATCCCATGACGGAACATTAAAGCAACAGCCTGGTTACGCCAATTGATAGCTGCGGCAACTAACGCAGCGGGCTTCAGGATTTCACCTTGTAACGACTATCaaacg                                >  2:426168/6‑249 (MQ=255)                                                                                                                                                                                                               cGGCCAAATAGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCACAAACCTCAACTCCGCTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGc                              >  2:324272/1‑250 (MQ=255)                                                                                                                                                                                                                    aaaTCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGCCGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgt                        >  2:416800/1‑251 (MQ=255)                                                                                                                                                                                                                     aaTCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTAc                                                                                                                                              <  2:50383/132‑1 (MQ=255)                                                                                                                                                                                                                     aaTCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTAc                                                                                                                                              >  1:50383/1‑132 (MQ=255)                                                                                                                                                                                                                        cGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTg                    >  1:135398/1‑251 (MQ=255)                                                                                                                                                                                                                            aTCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAAc                 <  1:114435/250‑1 (MQ=255)                                                                                                                                                                                                                              tccTGGGTTCCTATTCTACTAATGCTTTTTCGTCAAAAACCGCAACTCCGTTCCCATCAAATTCCGCATTTTCATTTAAATCATAAACTTTCAGCGTATTCTGCCTCATAAGGTGCCTCCTACCCTTGTTCCTTTGCCGGAACATTCAAGCAAAAGCCTGTGTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTg               <  2:653876/248‑1 (MQ=255)                                                                                                                                                                                                                                   ggTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAACTAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCg         >  2:478979/1‑251 (MQ=255)                                                                                                                                                                                                                                      tCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTGCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGTGTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGt       <  1:297630/250‑1 (MQ=255)                                                                                                                                                                                                                                          gTTCTACTAAGGCCTTTTCGTCAAAAACCTCACCTCCGTTCTCATCACATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCACGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTACAGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCtcc  >  2:328163/1‑249 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           GCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTAC  >  NC_000913/2867222‑2867705

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