Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,521,873	C→T	V50I (GTC→ATC)	yjhF ←	putative transporter YjhF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,521,873	0	C	T	96.5%	185.3 / NA	57	V50I (GTC→ATC)	yjhF	putative transporter YjhF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/28);  minor base G (0/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACGGGCAATAGCGAATAAAATAGGTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTAT  >  NC_000913/4521629‑4522118                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      tCACGGGCAATAGCGAATAAAATAGGTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGCGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAggg                                                                                                                                                                                                                                                 >  1:221036/1‑250 (MQ=255)                aaTAAAATAGGTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCa                                                                                                                                                                                                                                   >  1:454595/1‑250 (MQ=255)                  tAAAATAGGTACCAGAAGAATTAAACCGACTTCGAAACAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAAt                                                                                                                                                                                                                                >  1:78531/1‑251 (MQ=255)                    aaaTAGGTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATcc                                                                                                                                                                                                                              >  2:245910/1‑251 (MQ=255)                          gTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATAcc                                                                                                                                                                                                                         >  2:298088/1‑250 (MQ=255)                            aCCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACATACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCAc                                                                                                                                                                                                                       <  1:629024/250‑1 (MQ=255)                            aCCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATTCCAc                                                                                                                                                                                                                       <  2:112723/250‑1 (MQ=255)                               agaagaATTAAACCGACTGCGAAAAAAAGGGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTAGTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCTCCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCAATATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACGTTACTCAAATCCATACCACtggc                                                                                                                                                                                                                   <  2:328081/251‑1 (MQ=255)                                        aaaCCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCAATACCACTGGCGATGGCGAc                                                                                                                                                                                                          >  1:608062/1‑251 (MQ=255)                                          aCCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAg                                                                                                                                                                                                        <  1:309278/251‑1 (MQ=255)                                                  cGANAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCAGCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCgag                                                                                                                                                                                                <  2:608062/251‑1 (MQ=255)                                                  cGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATGCCACTGGCGATGGCGACAGCAATCgag                                                                                                                                                                                                <  1:245910/251‑1 (MQ=255)                                                   gAAAAAAAGTGCGATTCCGACAATAAACGCCGAACAGACCGCTGCCCAGTCAAGTTTATTTTTCCCGAAATAATGCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCAGCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCgaga                                                                                                                                                                                               <  2:403757/251‑1 (MQ=255)                                                   gAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCACTATGGCCCAGCGTACCGCCAAGTCCTGATTCGATTGAGGTGCCGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCg                                                                                                                                                                                                  >  2:146419/1‑248 (MQ=255)                                                       aaaaGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGAAATCGAGACaat                                                                                                                                                                                           >  2:619737/1‑251 (MQ=255)                                                         aaGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACaata                                                                                                                                                                                          >  1:579473/1‑250 (MQ=255)                                                                  cccGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAACGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGaaa                                                                                                                                                                                <  1:454054/250‑1 (MQ=255)                                                                   ccGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGAc                                                                                                                                                                                              <  1:212932/236‑1 (MQ=255)                                                                   ccGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGAc                                                                                                                                                                                              >  2:212932/1‑236 (MQ=255)                                                                             aCGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTAtt                                                                                                                                                                      >  2:226854/1‑250 (MQ=255)                                                                                 cGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCCCCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTaa                                                                                                                                                                    <  1:417675/248‑1 (MQ=255)                                                                                 cGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTaa                                                                                                                                                                    >  2:417675/1‑248 (MQ=255)                                                                                    aCAGACCACTGCCCAGTCCAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGAATTAAGCNTa                                                                                                                                                               >  2:594108/1‑250 (MQ=255)                                                                                      aGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAAtt                                                                                                                                                            <  1:525483/251‑1 (MQ=255)                                                                                       gACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAAttt                                                                                                                                                           >  1:575307/1‑251 (MQ=255)                                                                                       gACCACTGCCCAGTCAAGTTTATTTGTCCCGAAATAATTCATCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCATTCGAGACAATAATTAACGAAACAACCGTATTAAGCTTAAttt                                                                                                                                                           <  1:496489/251‑1 (MQ=255)                                                                                              gCCCAGGCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGGCCGGATTCGATTGAGGAGCCGACTTTACTCAAATCCATACCACTTGCGATGGCGACAGCAATCGATACAATAATTAACGCAACAAACGTATTAAGCTTAATTTTTATGGt                                                                                                                                                    >  2:13743/1‑251 (MQ=255)                                                                                                   gTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAAGCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTaaaa                                                                                                                                               <  1:102779/251‑1 (MQ=255)                                                                                                    tCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTaaa                                                                                                                                                >  1:242298/1‑249 (MQ=255)                                                                                                                ttCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAg                                                                                                                                   <  1:11487/250‑1 (MQ=255)                                                                                                                        ataataCAGCATCGTCAGCCCTATTCGCTGCGCACCACCCCCATCGGCCATCAAGCCACCGAGCCTAACCCCGGAGCCAAATAACAAACCCATAAGGCCCAGCGTTCCGCCAAGGCCGGATTCGAGTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACTTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTg                                                                                                                          <  2:384941/251‑1 (MQ=255)                                                                                                                        ataatTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTg                                                                                                                          <  1:226854/251‑1 (MQ=255)                                                                                                                             gcAGCATTGTCTGCGCTCTTCGCTGCGCACCTCCCGCCTCGTCTAGCAAGCCACCGAGCTTACCACCGAAGCCAAATATCAAACCAAGAGGGCCCCGCTTACCGCCAATTCCGGATTCGATTGAGTGGACGACTTTACTCAAATCCATACCTCGGGCTATGGCGATAGGAAGCGAGACAATAATTAGCGAAACAATCTGATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATTCCCGCcaa                                                                                                                      <  2:302416/249‑2 (MQ=255)                                                                                                                               aGCATGGTCAGCGCTATTCGCTGCGCACCACCAGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGCCCGGATTCGCTTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGTTGGTGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAAc                                                                                                                   >  1:549777/1‑251 (MQ=255)                                                                                                                                gCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACg                                                                                                                  >  1:410545/1‑251 (MQ=255)                                                                                                                                  aTGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGAt                                                                                                                >  2:215035/1‑251 (MQ=255)                                                                                                                                       cAGCGCTATTCGCTGCGCCCCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTa                                                                                                           <  1:155194/251‑1 (MQ=255)                                                                                                                                          cgcTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCAAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGCGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGATAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTg                                                                                                        >  1:474115/1‑251 (MQ=255)                                                                                                                                           gcTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTgg                                                                                                       >  1:560677/1‑251 (MQ=255)                                                                                                                                                tCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAATCAGGAGTAAAGCAATCCCTGCCACAACGATAATTATTGGCATa                                                                                                   <  1:564427/250‑1 (MQ=255)                                                                                                                                                 cgatgggaaccaccCGCAACGGCAGGAAAGGGACCGAGCATAACACCGAAGCAAAATAAAAAACCAATATGACCCAACGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACATTACTCAAATCCATACCACTGGAGAGGGCGAAAGAAATCGAAAAAATAATTAACGAAACAAACGTATTAAGCTTATTTTTTATGGTTAAAAGCAGGAGAAAAGCAGTCCCTGCCACAACGATAATTAGTGGCATAAc                                                                                                 <  2:485093/243‑1 (MQ=255)                                                                                                                                                   cTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGt                                                                                               >  1:447842/1‑251 (MQ=255)                                                                                                                                                     gcccaccacCCGCATCGGCCAGCAAGCAACCGAGCATAACTCCGAAGCCAAATATAAACCCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGAGTGAGGTGACGACTTTACTAAAATCAATACCACTGGCGATGGCGACAGCATTCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTc                                                                                              <  2:152384/247‑1 (MQ=255)                                                                                                                                                          ccaccCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCaa                                                                                                                                                                                                     >  1:659895/1‑142 (MQ=255)                                                                                                                                                          ccaccCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCaa                                                                                                                                                                                                     <  2:659895/142‑1 (MQ=255)                                                                                                                                                                 cATCGGCAAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGGCCGGATTCGATTGAGGTGACGACTTTACTCACATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAAAGATAATTAGTGGCATAACGTCCTCTCATAACAgg                                                                                 >  1:26063/1‑251 (MQ=255)                                                                                                                                                                    cGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGGGGCATAACGTCCTCTCATAACAGGGt                                                                               >  2:515134/1‑250 (MQ=255)                                                                                                                                                                          gCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGAcc                                                                        <  1:215035/251‑1 (MQ=255)                                                                                                                                                                             aGCGACCGAGCATAACACCAAAGCCAAATATAAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTAACGACTTTACACAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATAAAGCTTAATTTTTATGGTTAAAAGCAGGGGAAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCcat                                                                     <  2:287982/251‑1 (MQ=255)                                                                                                                                                                                gACCGAGCATAACCCCGAAGCCAAATAACAAACCAATAAGGCCCAGCGTACCGACAAGTCCGGATTCGATGGAGGTGACGATTTTCCTCAAATCCATACCACTGGCGATGGCGACAGCAAGCGAGACAATAATTAACGAAACAAACGTATTAAGCTAAATTGTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCcatcat                                                                  <  2:221036/251‑1 (MQ=255)                                                                                                                                                                                    gAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATc                                                              >  1:526395/1‑251 (MQ=255)                                                                                                                                                                                          aaCACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTg                                                        <  2:454595/251‑1 (MQ=255)                                                                                                                                                                                                   gCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATa                                                >  1:274666/1‑250 (MQ=255)                                                                                                                                                                                                            aaacccAAGATGGCCCAGCGTCCCGCCAATTCCGTGTTCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAatgatg                                      <  2:138280/247‑1 (MQ=255)                                                                                                                                                                                                                                           ccGGATTCGATTGATGTGACGACTTTTCTCAAATCCACACCACTGGCTATGGCGCCAGCAATCGAGACAATAATTAACGAAACAAACGTACTAAGCATAATTTTTATGGTTAAAAGCATGAGTCAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCCTCACAGGGTAGCGACCCATCATTACCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATg         <  1:120407/249‑1 (MQ=255)                                                                                                                                                                                                                                              gATTCGATTGAGGTGACGACTTTACTCAAATCCATACCACGGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACATATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGttttt    <  1:302145/251‑1 (MQ=255)                                                                                                                                                                                                                                                 tCGATTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTtttat  <  1:515134/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      TCACGGGCAATAGCGAATAAAATAGGTACCAGAAGAATTAAACCGACTTCGAAAAAAAGTGCGATACCGACAATAAACGCCGAACAGACCACTGCCCAGTCAAGTTTATTTTTCCCGAAATAATTCAGCATGGTCAGCGCTATTCGCTGCGCACCACCCGCATCGGCCAGCAAGCGACCGAGCATAACACCGAAGCCAAATATCAAACCAATATGGCCCAGCGTACCGCCAAGTCCGGATTCGACTGAGGTGACGACTTTACTCAAATCCATACCACTGGCGATGGCGACAGCAATCGAGACAATAATTAACGAAACAAACGTATTAAGCTTAATTTTTATGGTTAAAAGCAGGAGTAAAGCAATCCCTGCCACAACGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTAT  >  NC_000913/4521629‑4522118

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