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	95.8%	131.4 / NA	49	Q33* (CAG→TAG)	rpoS	RNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/26);  minor base T (0/2);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTG  >  NC_000913/2867208‑2867701                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTc                                                                                                                                                                                                                                                     <  1:131418/251‑1 (MQ=255)   cAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTcta                                                                                                                                                                                                                                                   <  1:294793/251‑1 (MQ=255)         ccGATTGCCATATCGGCGGGCAATTGTTACTCCCTGACTCAAGTTAGTCTCGACCAACCGGCGGCGAGAGGCTGCATCCCCACGCTGTGCGCGACGCGCCAAATAAACTTCTTCTGCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGGTCCTCTTCGGCCAAAGCGTTATCACTGGGTGCCTATTCTTCTAAgg                                                                                                                                                                                                                                             <  2:527198/250‑1 (MQ=255)               gCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCtttt                                                                                                                                                                                                                                       <  1:269706/251‑1 (MQ=255)                ccATAACGGCTGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGCGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTc                                                                                                                                                                                                                                      >  1:232485/1‑251 (MQ=255)                       ggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGCTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCCGCCAAATCGTTATCACTGGGTTCCTATTCTACTAACGCCTTTTCGTCaaa                                                                                                                                                                                                                                >  2:463339/1‑250 (MQ=255)                                     acccccataCGCAATTTACTCCCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGGTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAGCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt                                                                                                                                                                                                                 <  2:408410/243‑1 (MQ=255)                                                   gTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg                                                                                                                                                                                                   <  1:144248/251‑1 (MQ=255)                                                           cGATCATCCGGCGGCGAGAGGCGACATCTCAACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTc                                                                                                                                                                                            <  1:207286/250‑1 (MQ=255)                                                            gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCGAACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCCCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGAGTTCCTATTCTACTAAGGCCTTTTCGTCAAAACCCTCAACTCCGTTCTCATCAAATTCCACCTCTTCAt                                                                                                                                                                                          >  2:157369/1‑251 (MQ=255)                                                               catcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAttt                                                                                                                                                                                        <  1:212000/250‑1 (MQ=255)                                                                atcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa                                                                                                                                                                                       >  1:300776/1‑250 (MQ=255)                                                                                       cccACGCAGTGCGCGACGCTCACAATAAACTTCTTCTCCGGCCTTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCCCCCTGCGATAACAGTTCCTCCTCGTCCAAATCGTTATCACTGGGTTCCTATCCTCCTAAGGCCTCATCGTCAAAAACCTCAACTCCTTTCTCATCAAATTCCCCATCTCCATCCAAATCATGAACTTTCAGCGTTTTcg                                                                                                                                                               <  2:355702/250‑2 (MQ=255)                                                                                        ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGACTAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTATTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTCAATCATGAACTTTCAGCGTATTCTg                                                                                                                                                              >  1:366534/1‑251 (MQ=255)                                                                                        ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTg                                                                                                                                                              <  1:367896/251‑1 (MQ=255)                                                                                                     gACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTgg                                                                                                                                                 >  1:264916/1‑251 (MQ=255)                                                                                                                 aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                                                                                                                     >  1:209961/1‑251 (MQ=255)                                                                                                                        ttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCtt                                                                                                                              >  1:212123/1‑251 (MQ=255)                                                                                                                            tcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGTCTCCTACCCGTGATCCCTTGa                                                                                                                            <  2:457452/249‑1 (MQ=255)                                                                                                                                ccGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAAc                                                                                                                      >  1:308609/1‑251 (MQ=255)                                                                                                                                 cGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAAc                                                                                                                      >  1:311731/1‑250 (MQ=255)                                                                                                                                     aaCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTAGTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCa                                                                                                                 <  2:258049/251‑1 (MQ=255)                                                                                                                                                  aaCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTgg                                                                                                    <  1:95097/251‑1 (MQ=255)                                                                                                                                                      caTCTCCCCAAGGTAAAGTCGAGTCGCGTCCAACCCACGCTGTGTGGCTCCCTGCGATAACAGTTCCGCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGGGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                                                                                <  2:487090/250‑1 (MQ=255)                                                                                                                                                         ctcACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACgg                                                                                                                         >  2:483801/1‑223 (MQ=255)                                                                                                                                                         ctcACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACgg                                                                                                                         <  1:483801/223‑1 (MQ=255)                                                                                                                                                                   tAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGc                                                                                   >  1:214740/1‑251 (MQ=255)                                                                                                                                                                        cTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCgg                                                                              >  1:23781/1‑251 (MQ=255)                                                                                                                                                                                  tccgacacacGCTGTTTGTCTCCCTGCGATAACAGTTCCTCTTCAGCCACATCGTTATCCCTGGGTTCCTCTTCTACTACGGCCTTTTCGTCAACAACCTCAACTCCGTTCTCCTCACATTCCGCATCTTCATTTAAATCATGAACTTTCATCGTATTCTGACTCATCAGCTGGCTCATACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGTTTTATCGCTGCGGCAAATACCGc                                                                    >  2:80826/5‑251 (MQ=255)                                                                                                                                                                                  tCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGc                                                                    >  1:498091/1‑251 (MQ=255)                                                                                                                                                                                               gtgtgGCTCCCTGCGATAACTGTCCCTCTTCGGCCAAATCGTTAGCACTGGTGTCCTTTTCTACGAAGGCCTTTTCGTCAAAAACTTCAACTCCGTTCTCATCAAAGTCCGCATCTTCATGTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTTATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCTCTGCGGCAAAGAACGCAGCGGGTTTACg                                                        <  2:311731/250‑1 (MQ=255)                                                                                                                                                                                                   ggCTCCCTGCGATACCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTTTTCTTCTAAGGCCTTTTCGTCAAAATCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGATCTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAGGCCTGGTTCCGCCGATTTATCGCTTCGGCAAATAACGCAGCGGGTTTACGGAtt                                                    <  1:379786/250‑1 (MQ=255)                                                                                                                                                                                                      ccccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTccc                                                <  1:254491/250‑1 (MQ=255)                                                                                                                                                                                                           gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCTCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCtt                                             <  2:8369/249‑1 (MQ=255)                                                                                                                                                                                                               aaaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTTATAATGCCTTGTCGCTAAAAACCTCAACTCCGTTCTCATCAAATCCCGCATCTTCATTTAAATCATGATCTTTCAGCGTATTCTGACTCATAAGGGGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGTAAATAACGCAGCTGGTTTACGGATTTCCCCTTGTAACg                                       <  2:487529/250‑1 (MQ=255)                                                                                                                                                                                                                  cagtgcccctgctgCCAAATCGTTATCTATGGTTTCCTATTCTACTAAGGCCTTTTCTTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCATCGTATTCTGACTCATAAGTTGGCTCCTACCCGTGATCCCTTGACGTAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCATTGGAACGAAt                                    <  1:254146/238‑1 (MQ=255)                                                                                                                                                                                                                  cAGTTCCTCTTCGGCCCAATCGTTATCCCTGGGTCCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATATGGTGGCTCCTACCCGTGATCTCTTGTCGGAACATTCAAGCAAAAGCCTGGTTCCGCTGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAAt                                    <  1:261261/251‑1 (MQ=255)                                                                                                                                                                                                                  cAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAAt                                    <  1:487875/251‑1 (MQ=255)                                                                                                                                                                                                                   aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAAtt                                   >  1:233400/1‑251 (MQ=255)                                                                                                                                                                                                                      tccccgtCGGCCAAATCGTTATCACTTGGTATCTATTCTACTAAGGCCTTTTCTACAAAAACCTCAACTCCGATCTCATCAAATGCCGCATCGTACTTTAAATCATTAACTTTCAGCGTATTCTGACTCAAAAGGAGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTTACGAATTTc                                 <  2:358664/244‑1 (MQ=255)                                                                                                                                                                                                                          ctTCGGTCAAATAGTTATCACTGGGTTCCTATTCTACTAACGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGCTCCCTTGCCGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGAAGCGGGTTTACGGATTTCCCCTGATAACGAATTTCACAAt                            >  2:479637/1‑251 (MQ=255)                                                                                                                                                                                                                             cGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGc                          >  1:474821/1‑250 (MQ=255)                                                                                                                                                                                                                                        ttaccACTGGGTTCCTATTTCACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGGGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAAGGCAAGCGTGTTGa               <  2:308609/246‑1 (MQ=255)                                                                                                                                                                                                                                        ttATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGa               >  1:426925/1‑250 (MQ=255)                                                                                                                                                                                                                                         tATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAAc             >  1:405185/1‑251 (MQ=255)                                                                                                                                                                                                                                                gggTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGtt        <  1:236861/249‑1 (MQ=255)                                                                                                                                                                                                                                                gggTGCCTATTCTACTAAGGCCTGTTCGTCAAAAACCTCAACCCCTTTCTCATCAAACCCCTCATCTTCATTCAAATGTTGCACTGTTAGCGTATGCTGACTCATAAGGTGGCGCCTACCCGTGATCCCGTTACGGAACATGCAAGCAAAAGCCTGGTTCCGCCGTTCTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAGCGAAGTTCAAAATGCAAGCGTGTTGAACTGGtt        <  2:532923/249‑1 (MQ=255)                                                                                                                                                                                                                                                  gTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCgg    <  1:43379/251‑1 (MQ=255)                                                                                                                                                                                                                                                    tCCTATTCTACTAATGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCACTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTCTCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGACTTTCAAAATGCAAGCGTGTTGACCTGTTTCCGGTg  >  2:468373/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTG  >  NC_000913/2867208‑2867701

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: ‑