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%	175.2 / NA	50	Q33* (CAG→TAG)	rpoS	RNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/24);  total (26/24)

CACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTC  >  NC_000913/2867215‑2867696                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  cACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGGGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGGGGGGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGGTCCCTGCGATAACAGTTCCTCTTCGGCCCAATCGTTATCACTGGGGTCCTATTctacta                                                                                                                                                                                                                                           >  1:371970/1‑249 (MQ=255)              aCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCaa                                                                                                                                                                                                                            <  2:271780/251‑1 (MQ=255)                       aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGCGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTACCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACCAACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTAATACGGCCTTTTCGTCACAAACCTCa                                                                                                                                                                                                                    >  2:358472/1‑250 (MQ=255)                           tttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAACTCGTTATCACTGGGTTCCTATTCTACTAAGGCCtttt                                                                                                                                                                                                                                  >  2:277832/1‑232 (MQ=255)                           tttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCtttt                                                                                                                                                                                                                                  <  1:277832/232‑1 (MQ=255)                               ccaccaGACGCAAGGTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAACTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtc                                                                                                                                                                                                           >  2:127460/1‑251 (MQ=255)                                         cAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGCCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCCTCAAAtt                                                                                                                                                                                                 >  2:321442/1‑251 (MQ=255)                                         cAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt                                                                                                                                                                                                 >  2:365784/1‑251 (MQ=255)                                             ttACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGc                                                                                                                                                                                             >  1:426171/1‑251 (MQ=255)                                                  ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCtt                                                                                                                                                                                        >  1:168721/1‑251 (MQ=255)                                                     gATCACCCGGCGGCCAGAGGCGAGACATCCACGTTGTAAGCCACGCGCAAAATAATGTTCTTCTTAGACTCCTAACAGTCGTGAACAACATCTCTCACCAAGGTAAAGCTGAGACGCGTCTAACAAACGACGTGAGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGCTCTCAATGGGCTCCTATTCTACTAAGGCCTTTCCGTCAAAACCCTCAACTCCGTTCTCATCACATTCCTCATCTTCAt                                                                                                                                                                                     >  2:427288/1‑251 (MQ=255)                                                             ggccgcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGTTAAAGCTGACTCTCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTCCCTCTTCGGCCAAATCTGTATCACTTGGTTCCTATTCTACTAAGGCCCTCTCGTCAAAAGCCTCAACGCCGGTCTCACCAAATTCCGCATCGTCATTTAAATCa                                                                                                                                                                             <  2:278124/247‑1 (MQ=255)                                                                 gcgAGAGGCGACATCTCCACGCAGCGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGa                                                                                                                                                                          >  1:98055/1‑250 (MQ=255)                                                                    agagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACtt                                                                                                                                                                      >  1:386503/1‑251 (MQ=255)                                                                          gACATCTCCACGCTGTGCGCGACGCCAAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATACCCAATCTCACCAAGGTAAAGCTGCGTCGCGTCCCACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTACGGCCTTTTCGTCAAAAACCTCAACTCCGTTATCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCg                                                                                                                                                                >  2:138962/1‑251 (MQ=255)                                                                              cctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTAtt                                                                                                                                                            <  2:83333/250‑1 (MQ=255)                                                                                        gggcgcgaatcgcAAAACAAACTTCTTCTTCGGCCGTTAGCAGTGGTGAATAACCAATGTCACCCAGGTAAAGCTTAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCGTCGGCCAAATCATTATCACTGGGTTCCTATTCTACGAAGGCCTTTTTGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTTACTCATa                                                                                                                                                  <  2:132965/241‑1 (MQ=255)                                                                                         tGCGCGCCGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAACGCTAAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTAGGCCAAATCGTTATCACTGGGTTACTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTATCATCAAATTCCGCATCTTAATTTAAATCATGAACTTTCAGCGTATTCTGACTCATa                                                                                                                                                  >  2:139429/1‑250 (MQ=255)                                                                                               aCGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGc                                                                                                                                           <  2:369166/251‑1 (MQ=255)                                                                                                   gcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCt                                                                                                                                       <  1:387838/251‑1 (MQ=255)                                                                                                           aaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACTCTGTTTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGTCCTATTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGa                                                                                                                               <  1:321442/251‑1 (MQ=255)                                                                                                                      cGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGAAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACgg                                                                                                                    >  1:63181/1‑251 (MQ=255)                                                                                                                                  gtggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                              <  1:238831/181‑1 (MQ=255)                                                                                                                                  gtggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                              >  2:238831/1‑181 (MQ=255)                                                                                                                                         aTAACCTATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCt                                                                                                 >  1:130067/1‑251 (MQ=255)                                                                                                                                         aTAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCt                                                                                                 <  2:32963/251‑1 (MQ=255)                                                                                                                                          tAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTg                                                                                                >  1:101471/1‑251 (MQ=255)                                                                                                                                             ccAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTTATCCCTTGACGGAACATTCAAGCAAAAGCCTGGtt                                                                                             <  1:283958/251‑1 (MQ=255)                                                                                                                                             ccAATCGCACCAAGGGAAAGCTGAGGCGCGTCCAACACACGCTGTTGGGCTCCCTGCTATAACAGTTCCTCGTCGGCCAAATCGTTATCACTGGGTTCCTAGTCTACTAAGGCCTCTTCGTCAAAGACCTCAACGCCGTTCTCATCAAATTCCGCATCTTCATTTAAACCATGAACTTTCAGCGTATTCTGACTCAAAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAGTCAAGCAAAAGCCTCGtt                                                                                             <  2:285302/251‑1 (MQ=255)                                                                                                                                                        aaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                                                                                   >  1:474603/1‑250 (MQ=255)                                                                                                                                                                     gTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaa                                                                      >  1:80701/1‑250 (MQ=255)                                                                                                                                                                     gTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAAAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCACGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaa                                                                      >  2:445752/1‑250 (MQ=255)                                                                                                                                                                       cgcgTCCAACACACGCTTTGTGGCTCCCTGCGATACCAGTGCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCTCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCTTGATCCCTTGACGGAACATTCAAGCAAAATCCTGGTTCCGCCGATTTATCGCTGCGGCAAATa                                                                   <  1:358472/251‑1 (MQ=255)                                                                                                                                                                         cgTCCAACACACGCTGTGTGGCTCCCTGCGCTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAc                                                                 >  2:222850/1‑251 (MQ=255)                                                                                                                                                                                  cacGCTGGGGGGCTCCCTGCGATAACAGTTCCTCTGCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGGGACCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGt                                                        <  2:426171/251‑1 (MQ=255)                                                                                                                                                                                        gtttgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTGCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCTCTGCGGCAAATAACGCAGCGGGTTTACgg                                                  <  2:404189/248‑1 (MQ=255)                                                                                                                                                                                        gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACg                                                   <  2:130067/250‑1 (MQ=255)                                                                                                                                                                                          gtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCGGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGAAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGa                                                 <  1:445752/250‑1 (MQ=255)                                                                                                                                                                                                 ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATACCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCt                                         >  1:293827/1‑251 (MQ=255)                                                                                                                                                                                                            aGTTCCCCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAAtt                              <  2:386503/251‑1 (MQ=255)                                                                                                                                                                                                             gTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAACGCCTTTTCGTCAAAAACCTCAACTCCGTTCTAATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCACTTGCCGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAAttt                             >  1:285816/1‑251 (MQ=255)                                                                                                                                                                                                                ccTCTGCGGACAAATCGTTATCCCCTGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACGCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGATCTTTCAGCGTATTCTGCCTCATAAGGTGGCTCCTACCCGTGATCCCTTGGCGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATATCGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCa                           <  2:11839/250‑1 (MQ=255)                                                                                                                                                                                                                      cGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAAc                                                                                                                 <  2:376263/158‑1 (MQ=255)                                                                                                                                                                                                                      cGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAAc                                                                                                                 >  1:376263/1‑158 (MQ=255)                                                                                                                                                                                                                           aaaTCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgt               <  1:365784/251‑1 (MQ=255)                                                                                                                                                                                                                               cgttatcccctggTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCCCATCAAATTCCGCATCTTCATTTAAACCATGAACTTTCAGCGGATTCTGACTCATAAGGTGGCTCCTCCCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGTTTCCGCCGATTTATCGCTTCGGCAAATAACGCAGCGGGTTTACGGAGTGACCCTTGTAACGAATTTCAAAATGCAAGCGTGTTg           <  2:159204/240‑1 (MQ=255)                                                                                                                                                                                                                                gTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTg           >  1:21067/1‑250 (MQ=255)                                                                                                                                                                                                                                  tATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGaa         <  2:80701/250‑1 (MQ=255)                                                                                                                                                                                                                                     cACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTg      <  2:21067/250‑1 (MQ=255)                                                                                                                                                                                                                                        tGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTc  >  1:15253/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  CACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTC  >  NC_000913/2867215‑2867696

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