Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	68,476	A→G	C525R (TGC→CGC)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	68,476	0	A	G	98.2%	174.6 / NA	55	C525R (TGC→CGC)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (33/21);  minor base T (0/1);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACGACGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGCACT  >  NC_000913/68232‑68721                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      tAACAACGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGtt                                                                                                                                                                                                                                                  <  2:258264/250‑1 (MQ=255)  aacaacgTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGc                                                                                                                                                                                                                                                >  2:48610/5‑251 (MQ=255)  aacaacgTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGc                                                                                                                                                                                                                                                >  2:309899/5‑251 (MQ=255)      acgTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAggg                                                                                                                                                                                                                                            >  1:138585/1‑251 (MQ=255)           cTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTtc                                                                                                                                                                                                                                       >  1:230850/1‑251 (MQ=255)             cGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAAACACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTtctc                                                                                                                                                                                                                                     >  1:299630/1‑251 (MQ=255)                   gTTGGGTTAGCTGACGCAGGGTTTCCGGGCCAGACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTGTCAAAGATCGTCATTATCGTGTCCGTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGGTGTTCGGCGCTCAGCGCCCATTGCTGATAGCGGCGATTAAGGTGTTCAAGGCGTTGGGCGTGTTGGCTGCGCGGTTGCAGGGTTTTCTCGACCGc                                                                                                                                                                                                                               <  1:345561/251‑1 (MQ=255)                               gACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCAttt                                                                                                                                                                                                                    <  1:175311/250‑1 (MQ=255)                                      ggTTTCCGGGCCTTGCAGAGGCTGGCTGCCAATGACAAACCCCACTTCATAATGATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGAGAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACTGCACTGGCCGTTTTTCgctgt                                                                                                                                                                                                            <  2:249840/251‑2 (MQ=255)                                            cGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGAt                                                                                                                                                                                                      >  1:109416/1‑251 (MQ=255)                                              ggCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATgg                                                                                                                                                                                                    >  1:52155/1‑251 (MQ=255)                                                       gATGCTGGCTGCCAATGACAAACCACACTTCATCATTCTCAAAAATCGTCATTACCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGACGGTGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGCTAGCGGCGATACAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCCGTGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGCTGGGATGTCTGc                                                                                                                                                                                           >  2:161264/1‑251 (MQ=255)                                                             ggctgCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCCCTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCAc                                                                                                                                                                                     >  1:221567/1‑251 (MQ=255)                                                               ctgCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACtt                                                                                                                                                                                   >  1:292636/1‑251 (MQ=255)                                                                   cAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTcgc                                                                                                                                                                               >  2:132578/1‑251 (MQ=255)                                                                    aaTGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTcgcg                                                                                                                                                                              >  1:176158/1‑251 (MQ=255)                                                                        aCAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTGCGCGGCGa                                                                                                                                                                          <  1:324751/251‑1 (MQ=255)                                                                         cAAACCACCCTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGAc                                                                                                                                                                         <  1:331528/251‑1 (MQ=255)                                                                          aaaCCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACg                                                                                                                                                                        >  2:66520/1‑251 (MQ=255)                                                                                     cATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAAt                                                                                                                                                             <  1:158543/251‑1 (MQ=255)                                                                                     cATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAAt                                                                                                                                                             >  1:96499/1‑251 (MQ=255)                                                                                            aTCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCAtttttt                                                                                                                                                                                                                 <  2:158244/192‑1 (MQ=255)                                                                                            aTCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCAtttttt                                                                                                                                                                                                                 >  1:158244/1‑192 (MQ=255)                                                                                                  aaTCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAgc                                                                                                                                                >  2:74514/1‑251 (MQ=255)                                                                                                         aTTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCCGGGCAGAAGTTGCACGATAGTGTTGTTCGGCGCTCACCGCACATTGCTGATAGCGGAGATAAACATGTTCAAAGCGTTGTGCCTGTTCGCAGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTACACTTTAGCGGCTACGGCAGCAAAAATCGCCTCACCGCGAGCACAGc                                                                                                                                         >  1:248095/1‑251 (MQ=255)                                                                                                                     ttATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAgg                                                                                                                             >  2:149422/1‑251 (MQ=255)                                                                                                                       ataGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGCTGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAACTCGCCGCACCGAGCGCACAGCCCTGGTCAGAGGc                                                                                                                            >  2:320477/1‑250 (MQ=255)                                                                                                                        taGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATACAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGTGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAACATCGCCGCACCGAGCGCACAGCACTGGTCAGCGGca                                                                                                                           >  2:257401/1‑250 (MQ=255)                                                                                                                              cGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTGGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAAttt                                                                                                                    <  2:307955/251‑1 (MQ=255)                                                                                                                                             aGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTc                                                                                                      >  2:365778/1‑250 (MQ=255)                                                                                                                                                  gtgCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACg                                                                                                <  1:82344/251‑1 (MQ=255)                                                                                                                                                             gAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGc                                                                                     <  1:132578/251‑1 (MQ=255)                                                                                                                                                               aGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGAGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCCACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGcac                                                                                   >  2:103472/1‑249 (MQ=255)                                                                                                                                                                gTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTg                                                                                  <  2:138585/251‑1 (MQ=255)                                                                                                                                                                    gAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATa                                                                              <  2:310934/251‑1 (MQ=255)                                                                                                                                                                             gttgttCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTgg                                                                     <  1:329007/251‑1 (MQ=255)                                                                                                                                                                                       gcTCATCGCCCATTGCTGTTAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCcgcg                                                           >  1:225305/1‑251 (MQ=255)                                                                                                                                                                                         tCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATCCAGCACGTAGCAGCAGGCCTGCATAATGACCTGGTTTTTCcgcgcg                                                         >  2:310735/1‑251 (MQ=255)                                                                                                                                                                                          cATCGCCCATTGCTGTTAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGa                                                        <  2:225305/251‑1 (MQ=255)                                                                                                                                                                                              gCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGCCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATgcc                                                    >  2:6933/1‑251 (MQ=255)                                                                                                                                                                                              ccccATTGCTGATAGCGGCGATTAATCTTTTCAAAGCGTTGTGCCTGTGCGCTGCTCGGTTTCAGGGTTTTCTCGACCGCTGTGGCCATGGTTGGCTTGGCTGATGGGATGTCTGCGGGCACTTTCGCGGCGTCGGCAGCAAAAATCGCCTCACCGTGCGCACATCACTGTTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCTCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATgcc                                                    <  1:176169/250‑1 (MQ=255)                                                                                                                                                                                                   ttGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGccc                                                >  2:151223/1‑250 (MQ=255)                                                                                                                                                                                                      cTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAgcg                                            >  1:11631/1‑249 (MQ=255)                                                                                                                                                                                                           aGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTGCCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACCTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCAt                                        <  2:11631/250‑1 (MQ=255)                                                                                                                                                                                                             cggcgATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCAc                                     <  1:151223/251‑1 (MQ=255)                                                                                                                                                                                                                 gATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGtta                                 <  1:309899/251‑1 (MQ=255)                                                                                                                                                                                                                         tGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGc                         >  1:147413/1‑251 (MQ=255)                                                                                                                                                                                                                          gTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCg                        >  1:52432/1‑251 (MQ=255)                                                                                                                                                                                                                               aaGCGTTGTGCCTGTTCGCGGCGGGGTTGCAGGGTTTCCTCTCCCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATccc                   <  2:269250/251‑1 (MQ=255)                                                                                                                                                                                                                                      gtgCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCg            >  1:206160/1‑251 (MQ=255)                                                                                                                                                                                                                                            gttCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTaaa       <  1:223024/248‑1 (MQ=255)                                                                                                                                                                                                                                             ttCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGc     >  2:366402/2‑251 (MQ=255)                                                                                                                                                                                                                                              tCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGc     <  2:292636/250‑1 (MQ=255)                                                                                                                                                                                                                                               cGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCAc   >  2:134751/1‑251 (MQ=255)                                                                                                                                                                                                                                                gCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCACt  >  1:141644/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      TAACGACGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTGCAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGCACT  >  NC_000913/68232‑68721

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