Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,479,012	A→T	9.5%	intergenic (+123/+22)	rraB → / ← yjgM	protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,479,012	0	A	T	9.5%	134.0 / 4.2	47	intergenic (+123/+22)	rraB/yjgM	protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (26/16);  new base T (2/3);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

TTGAATGCCGATCTGATCGATGCCCAGGTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCCCGGTTAATGGCG  >  NC_000913/4478725‑4479305                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                      tggagttcagacgtgtgctcttccgatcttGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCgg                                                                                                                                                                                                                                                                                          <  2:546313/273‑1 (MQ=255)       gCCGATCTGATCGTTGCCCAGGTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttc                                                                                                                                                                                                                                                                                    <  1:188992/301‑2 (MQ=255)       gCCGATCTGATCGATGTCCAGGTTGTACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACGTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCGGGCAttc                                                                                                                                                                                                                                                                                    <  2:276901/301‑2 (MQ=255)                              tGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggagatcggaagagcgtcgtgtagggaaaga                                                                                                                                                                                                                                                             >  1:545640/1‑272 (MQ=255)                                      tgattaCGCTGTCAGATAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGGATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCACACTTcac                                                                                                                                                                                                                                                     <  1:214644/301‑1 (MQ=255)                                                       aaaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                                                                                                                                                                                                                    >  2:298471/1‑301 (MQ=255)                                                                     gTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTAGCGCAACATCCACACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                                                                                                                                                                                                                      >  2:164453/1‑301 (MQ=255)                                                                                    cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGGCGATATGCTCaaa                                                                                                                                                                                                       >  1:287877/1‑301 (MQ=255)                                                                                        aagttggAAGATCCTAAGGGTGAAGATGGGGA‑GGTGGAGTGTTTGGCGGATGAGGAGGAGGACGAGGTCGGCCACTATGAAAGTGTGGCGAACAGGGCCAAGGTAGTTCCGGCAACTGGGGAAATTTGGAGCCCGGTTAAGGCATTCAGGGAACATCCGGCATAACAAAAACGTAGGTGGTTAAGGTGCGAAGGCCCGCGAGTGGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGccc                                                                                                                                                                                                   <  1:105547/295‑1 (MQ=255)                                                                                              gAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACCTCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCCCAAAGCCCAAATGc                                                                                                                                                                                             >  1:497307/1‑301 (MQ=255)                                                                                                 gATCCTTACGGCGAAGATGGTGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCAACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTGGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCa                                                                                                                                                                                          <  1:6794/301‑1 (MQ=255)                                                                                                          ggCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGATGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGAAATGCTCAAAGCCCAAAAGCACAAAAAACGCa                                                                                                                                                                                 >  1:393613/1‑301 (MQ=255)                                                                                                            cGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGAAATGCTCAAAGCCCAAATGCTCATAAAGCGCaaa                                                                                                                                                                               >  2:565901/1‑300 (MQ=255)                                                                                                            cGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAACGCATAGACGAAATGCTCACAGCCCAAATGCCAAAAAAGCGCAAt                                                                                                                                                                               >  1:122049/1‑301 (MQ=255)                                                                                                                                          ttGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGATGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCa                                                                                                                                                  <  2:394127/300‑1 (MQ=255)                                                                                                                                            gAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGACCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGa                                                                                                                                                >  2:553738/1‑301 (MQ=255)                                                                                                                                                             tCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAAATAGCAGCGTTTGAAAcc                                                                                                                               >  1:224256/1‑301 (MQ=255)                                                                                                                                                                     aaTTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCAtctctc                                                                                                                       >  2:110533/1‑301 (MQ=255)                                                                                                                                                                       atAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTcgc                                                                                                                     <  1:110409/300‑1 (MQ=255)                                                                                                                                                                        tAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTcgcg                                                                                                                    >  2:475210/1‑301 (MQ=255)                                                                                                                                                                         ctttcccttcacgacgctcttccgatctCCTTCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCgcgc                                                                                                                   <  1:409621/274‑1 (MQ=255)                                                                                                                                                                          aGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGcc                                                                                                                  <  2:224499/301‑1 (MQ=255)                                                                                                                                                                                      tcagtcgtgtgctcttccgatctcTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGcc                                                                                                      <  2:55078/278‑1 (MQ=255)                                                                                                                                                                                          cccAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTATCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCCCCATCGCCCttt                                                                                                  >  1:635895/1‑300 (MQ=255)                                                                                                                                                                                                 aGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAg                                                                                           <  2:497954/301‑1 (MQ=255)                                                                                                                                                                                                  gTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAga                                                                                          >  2:211512/1‑300 (MQ=255)                                                                                                                                                                                                     ccTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCgcgcagatcggaagagcacacgtctgaactcc                                                                                       >  2:410184/1‑273 (MQ=255)                                                                                                                                                                                                             cTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGATCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCACCATCGCCagaacagaagagcctcgtgtagg                                                                               >  1:55010/1‑279 (MQ=255)                                                                                                                                                                                                                ttGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGccc                                                                            <  2:122173/300‑1 (MQ=255)                                                                                                                                                                                                                                    aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGgaa                                                        >  2:604881/1‑301 (MQ=255)                                                                                                                                                                                                                                    aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTACTCCATCGCCCATAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGgac                                                        >  1:583606/1‑300 (MQ=255)                                                                                                                                                                                                                                    aGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATACTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGgaa                                                        >  2:630918/1‑301 (MQ=255)                                                                                                                                                                                                                                      gcatgcatGTCACATCCGGCATAAACAAAACGTACTTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGgaaga                                                      <  1:565187/301‑1 (MQ=255)                                                                                                                                                                                                                                                cacaTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGAAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTAGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCAAGCCCTTTGCCACGGACAGCCGACAGAAAATACAACt                                            >  1:160486/1‑301 (MQ=255)                                                                                                                                                                                                                                                   aTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACAACtttt                                         >  1:204064/1‑301 (MQ=255)                                                                                                                                                                                                                                                     ccgtttgctcttccgatctcgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGc                                       <  2:600709/282‑1 (MQ=255)                                                                                                                                                                                                                                                          aTAAACAAAACGGCCGTTTTTAATTTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATACTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTc                                  <  1:630041/301‑1 (MQ=255)                                                                                                                                                                                                                                                              aCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGTAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCAGACAGAAAAAACAACttttcaagtacacaa                              >  1:441645/1‑290 (MQ=255)                                                                                                                                                                                                                                                                     ttccgtttgttttaaaaaaaTGCCCGCGTGTTCAACTCGCGGTCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCg                       <  1:113329/283‑1 (MQ=255)                                                                                                                                                                                                                                                                        cgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGAAAGCCGGAAGAAAATAAAACTTTTGCAGATCGGaagaacgtcgtg                    >  1:599948/1‑291 (MQ=255)                                                                                                                                                                                                                                                                          ttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGactca                  >  2:553951/1‑301 (MQ=255)                                                                                                                                                                                                                                                                              ctctaaaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGagatcggaa              >  2:113453/10‑293 (MQ=255)                                                                                                                                                                                                                                                                                 tGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCCCgg           <  2:160637/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                   aaaaaaTGCCCGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATACCCGGAAGAACACCCAACTTTTGCAGTTCGAACATAcccacaccactccccttt         >  1:11543/5‑286 (MQ=255)                                                                                                                                                                                                                                                                                    tCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGCCGTTTCCCGCTAGCAGCGTTTGAAACCCATCTCTCGAGCCTGCTCCATCGCCATTAAGGCCAGATTTTTTGCCAGCCCTTTGCCGCGGATCGCCGGAAGAAAAGACACCTATTGCCGTTCACAAATAACCTACTCACCCCCGGttt        >  1:484433/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                     cGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATAACAGACTCACTCCCGGTTaa       >  2:332467/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                          gCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCCCGGTTAATGgcg  <  2:204307/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                      TTGAATGCCGATCTGATCGATGCCCAGGTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTC‑GATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCCCGGTTAATGGCG  >  NC_000913/4478725‑4479305

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑