Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,867,454	T→A	Q33L (CAG→CTG)	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,454	0	T	A	98.0%	158.4 / NA	49	Q33L (CAG→CTG)	rpoS	RNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (22/26);  minor base G (1/0);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCC  >  NC_000913/2867219‑2867697                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCt                                                                                                                                                                                                                                      <  1:51258/251‑1 (MQ=255)          aCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCaa                                                                                                                                                                                                                             <  1:463483/251‑1 (MQ=255)              cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGCGCGGCGACATCTCCACGCAGTGCGAGACGCGCAAAATAACCTTCTTCTTCGGCCGTTAACAGTGGAGAATAACCAATCTCACCAAGGTAAAGCTGCGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAAc                                                                                                                                                                                                                         >  2:222256/1‑251 (MQ=255)                     tttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt                                                                                                                                                                                                                  >  2:157990/1‑251 (MQ=255)                       tttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAGGGTAAAGCTGAGTCGCGTCCACCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGTCCAAATCGTTATCACTGGGTTCCGGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                                                                                                                                                                                                                >  1:564932/1‑251 (MQ=255)                         tACCACCAGACGCAAGTCACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGGGCGCGACGCGCAAAATAAACTTCTTCTGCGGCCGTTAACAGTGGTGAATAACCAATCTCCCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTCCGGCCAAATCGTTATCTCTTGGTTCCAGTTCTGCTGAGGCCTTTTCGTCAAAAACTTCAACTCCg                                                                                                                                                                                                               <  2:582659/250‑1 (MQ=255)                              ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAATCTGAGTCGCGTCCAACGCACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTGCAGTTCTACTAAGTTCTTTTCGTCAAAAACCTCAACTCCGTTCtca                                                                                                                                                                                                         <  1:254536/251‑1 (MQ=255)                                     cAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAtt                                                                                                                                                                                                  >  1:470175/1‑251 (MQ=255)                                           aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                                                                                                                                                            <  1:377974/251‑1 (MQ=255)                                            ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATc                                                                                                                                                                                           >  2:532481/1‑251 (MQ=255)                                             tctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt                                                                                                                                                                                          >  1:30518/1‑251 (MQ=255)                                                    catcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa                                                                                                                                                                                   <  1:147243/251‑1 (MQ=255)                                                          gcggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCAt                                                                                                                                                                             <  1:488767/251‑1 (MQ=255)                                                             gcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGaa                                                                                                                                                                          <  1:532481/251‑1 (MQ=255)                                                                 gagGCGACATCTCCACGCAGTGCGCTACGCGCAAAATAAACCTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGTCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACttt                                                                                                                                                                      <  2:166661/251‑1 (MQ=255)                                                                          tctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGATGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTCAATCATGAACTTTCAGCGTAtt                                                                                                                                                             >  1:41550/1‑251 (MQ=255)                                                                             ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACtt                                                                                                                                                                       <  2:307001/238‑1 (MQ=255)                                                                             ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACtt                                                                                                                                                                       >  1:307001/1‑238 (MQ=255)                                                                              cACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGAGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGAGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATAGTTATCACTGGGATCAAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTAAAATCATGAACTTTCAGCGTATTCTg                                                                                                                                                          >  1:377658/1‑250 (MQ=255)                                                                              cACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTGCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCGTCATTTAAATCATGAACTTTCAGCGTATTCTGa                                                                                                                                                         <  2:485219/251‑1 (MQ=255)                                                                                   aGTGCGCGACGCGCAAACTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAAACCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCCTTTAAATCATGCACTTTCAGCGTATTCTGACTCAt                                                                                                                                                    >  2:87359/1‑251 (MQ=255)                                                                                         cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTg                                                                                                                                              >  2:76244/1‑251 (MQ=255)                                                                                             gcgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTc                                                                                                                                          >  1:149642/1‑251 (MQ=255)                                                                                               gcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCt                                                                                                                                        <  1:467134/251‑1 (MQ=255)                                                                                                  aaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCCGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGCCTCATAAGGTGGCTCCTAcc                                                                                                                                     <  2:45256/251‑1 (MQ=255)                                                                                                     gaaaCTTCTGCTTCGGCCGTTAACAGTGGTTAATAACCAATCTCTCCGAGGTAAAGCTGAGTCGCGTTCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAGCACTTGGTTCCAGTTCTACTAAGGCCTTTTCTTCAAAAACCTAAACGCCGTTCTCATCAAATTCCGCATCTTCAGTTAAATCTTGATCTTTCAGCTTATTCTGTCTCATAAGTTggtttccacccg                                                                                                                                   <  2:239942/249‑11 (MQ=255)                                                                                                      aaccttcttcttCGGCCGTAAACAGAGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGCTAACAGTTCCTCTTCGGCCAAATCGTTATCCCTGGGTTCCAGATCTACTAAGGCCTTTTCGTCAAACACCTCAACTCAGTTCTCATCAAATTCCGCATCTTCATTTAAATCCTGCACTTTCACCGCATCTTGACTCATAAGGTGGCTCCTTCCCGTg                                                                                                                                 >  2:485938/4‑251 (MQ=255)                                                                                                           tcttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAATGTAAAGCTGAGTCGCGTCCAACCCACGCTGTGTGGCTCCCTTCGATAACAGTTCCTCTTCGTCCAAATCGTTATCACTGGGTTCCAGTTCTTCTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAAGCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTCCCCGTGATcc                                                                                                                             <  2:531859/250‑1 (MQ=255)                                                                                                                  cGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGATACCTGCGATAACAGCTCCTCTTCGGCCAAATCGTTATCACTGGGCTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACgg                                                                                                                     >  1:549457/1‑251 (MQ=255)                                                                                                                    gCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGa                                                                                                                    >  2:195414/1‑250 (MQ=255)                                                                                                                         tAACAGTGGTGAATAACCAATCTCACCAAGGGAAAGCTGCGCCGCGTCCAACACACGCTGGTTGGCGCCCTGCGATAACAGTTCCTCTCCGGCCAAATCGTGATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAATTAAAGCATGAACAATCAGCGTATTCTTACTCATAAGGTGGCTCCTACCCGTGATCCCCTGACGGAACAtt                                                                                                               <  2:564932/250‑1 (MQ=255)                                                                                                                                ggtgAATAACCAATCGCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTGCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAGCCGCAACTCCGTTCGCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTGGACGGAACATTCAAGCaa                                                                                                        <  1:57739/250‑1 (MQ=255)                                                                                                                                       aaCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACTCACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATGCCGCATCTTCATTTAAATCATGAACTTGCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTg                                                                                                 <  1:331974/250‑1 (MQ=255)                                                                                                                                       aaCCAAACTCACCAAGGTAAAGCTGAGTCGCGTCCCACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTGCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAGTCCGCATCTTCATGTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCGTGACTGAACAGTCAAGTAAAAGCCTg                                                                                                 <  2:248955/250‑1 (MQ=255)                                                                                                                                           aaTCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                                                                            <  2:394243/251‑1 (MQ=255)                                                                                                                                                   ccgacgcccttccgatctcgTCCAACACCCGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTGCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCATCGTAGTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTGGACGGAACATTCAAGCAAAGGCCTGGTTCCGCCGAttt                                                                                    <  2:90366/233‑1 (MQ=255)                                                                                                                                                    aaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTa                                                                                   <  1:195414/251‑1 (MQ=255)                                                                                                                                                      ggTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGCGGCTCCTACCCGTGATCCCTTGTCGGAACATTCAAGCAAAAGCATGGTTCCGCCGAGTTATc                                                                                 >  1:179744/1‑251 (MQ=255)                                                                                                                                                      ggTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCGCAACTCCGTTCTCATCAAATGCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCTTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCTCCGATTTATc                                                                                 <  2:543877/251‑1 (MQ=255)                                                                                                                                                       gTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGACTGGTTCCGCCGATTTATCg                                                                                >  1:487416/1‑251 (MQ=255)                                                                                                                                                          aaGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTg                                                                             >  2:214862/1‑251 (MQ=255)                                                                                                                                                             cTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCgg                                                                          >  1:2342/1‑251 (MQ=255)                                                                                                                                                               gAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGTAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGc                                                                         >  1:72905/1‑250 (MQ=255)                                                                                                                                                                   cgcttCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATa                                                                    <  2:149642/247‑1 (MQ=255)                                                                                                                                                                     cgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                                                                                    >  1:90366/1‑233 (MQ=255)                                                                                                                                                                                                       cAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAAt                                >  1:344731/1‑251 (MQ=255)                                                                                                                                                                                                           tCCTCTTCGGCCAAATCGTTATAACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTAAAAAACCTCAACTCCGTTCTCATAAAATTCCGCATCTTCATTTAAATCATGAAATTTCAGCGTATTCTGACTCATAAGTTGGCTCCTACACGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGAATTTCCCCTTGTAACGAATTTc                             >  2:252728/1‑250 (MQ=255)                                                                                                                                                                                                                 tCGGCCAAATCATTATCACTGGGTTCCAGTGCTACTAATGCCTGTGAGCCAAAAACCTCAACTCCGTTCTCATCAAATTCCTCAGCTTCATTTAAATCATGAACTTTCAGCGTATTCGGACTCATAAGGTGGCTCCTTCCCGTGATCTCTGTACGGAACATTCAAGCAAAAGGCTGGTGCCGGCGATTTATGGCGTCGGCAAAGAACGCAGCGGGTTTACGGATTTCCCCTTGTAACTAAGTTCAAAAt                        <  2:77919/249‑1 (MQ=255)                                                                                                                                                                                                                                     gggTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTcc  <  1:250506/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    ATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCC  >  NC_000913/2867219‑2867697

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