Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,585,149	Δ1 bp	5.3%	coding (1338/2280 nt)	ydeP ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,585,149	0	G	.	5.3%	260.8 / 7.4	76	coding (1338/2280 nt)	ydeP	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (39/32);  new base . (2/2);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

AGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCG  >  NC_000913/1585001‑1585297                                                                                                                                                     |                                                                                                                                                      agCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTg                                                                                                                                                       >  8:267725/1‑149 (MQ=255) agCGGTTAAGCTTAGTGCCTACGTGTACCGCCAAATCTAATTGCGTTTACGGTACAGCGCTAGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTg                                                                                                                                                       <  6:330053/149‑1 (MQ=255)     cttACGCTTAGTGGCTACGTGTACAGCCAAATCAACTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCaa                                                                                                                                                   <  2:363082/148‑1 (MQ=255)      ttAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAAt                                                                                                                                                  <  2:206081/149‑1 (MQ=255)      ttAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAAtt                                                                                                                                                 <  1:212747/149‑1 (MQ=255)            gtactGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCGCGCGTCCCGATCTGGCCTTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cag                                                                                                                                            <  1:275478/144‑1 (MQ=255)              aGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cagc                                                                                                                                           <  2:56398/147‑1 (MQ=255)               gTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATg                                                                                                                                        >  4:177290/1‑149 (MQ=255)                      cGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAgg                                                                                                                                  <  4:227332/148‑1 (MQ=255)                      cGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAg                                                                                                                                   <  1:64941/147‑1 (MQ=255)                                cAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATg                                                                                                                            >  4:216608/1‑144 (MQ=255)                                cAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATg                                                                                                                            <  3:216608/144‑1 (MQ=255)                                cAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTg                                                                                                                       <  3:117472/149‑1 (MQ=255)                                      tAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCAGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTg                                                                                                                  >  8:391063/1‑148 (MQ=255)                                         ttGCGTTAACGGTACAGCGCTCGCTTCCCGCTCAGGCATTCCCAGCGCAAAAATGCCCCCCATTCAGATGAATGCTCGAGCCTGCCCCGTACATATTGCATGCATGCTGGCAATTG‑CAGCATCTCCAGGTCCATGAGGTGGTGTGAAGc                                                                                                              >  2:93384/1‑149 (MQ=255)                                         ttGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CCGCATGTCCAGGTGCAGGAGGTGGGGTGAAGc                                                                                                              >  7:336898/1‑149 (MQ=255)                                             gTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCa                                                                                                            >  7:234303/1‑147 (MQ=255)                                             gTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                                                                          <  2:260024/149‑1 (MQ=255)                                                aaCGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑ca                                                                                                                                             <  7:382148/111‑1 (MQ=255)                                                aaCGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑ca                                                                                                                                             >  8:382148/1‑111 (MQ=255)                                                aaCGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcg                                                                                                       >  5:339341/1‑149 (MQ=255)                                                aaCGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATgcaa                                                                                                                                     <  2:96558/119‑4 (MQ=255)                                                aaCGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATgcaa                                                                                                                                     >  1:96558/1‑116 (MQ=255)                                                     tACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                                                                                                  >  6:324133/1‑149 (MQ=255)                                                       cAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAcc                                                                                                 <  7:267725/148‑1 (MQ=255)                                                          cgctcgcgTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCa                                                                                               <  1:93384/147‑1 (MQ=255)                                                           gctcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGAc                                                                                            <  8:328233/149‑1 (MQ=255)                                                             tcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGAc                                                                                            >  6:301919/1‑147 (MQ=255)                                                                           ggCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTggg                                                                                                                     <  4:279225/108‑1 (MQ=255)                                                                           ggCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTggg                                                                                                                     >  3:279225/1‑108 (MQ=255)                                                                            gCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTg                                                                             <  3:378846/147‑1 (MQ=255)                                                                            gCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCa                                                                           <  3:177290/149‑1 (MQ=255)                                                                              aTTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAgg                                                                                                                                  >  6:264179/1‑92 (MQ=255)                                                                              aTTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAgg                                                                                                                                  <  5:264179/92‑1 (MQ=255)                                                                              aTTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGa                                                                         >  3:13532/1‑149 (MQ=255)                                                                              aTTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGa                                                                         >  2:261306/1‑149 (MQ=255)                                                                                tGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgggg                                                                                                                    >  3:304746/1‑104 (MQ=255)                                                                                tGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgggg                                                                                                                    <  4:304746/104‑1 (MQ=255)                                                                                  ccAGCGCAAAATTGCCGCCGATGCTGATCAATGCTGGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCTTGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCTGACGGt                                                                     <  8:234303/149‑1 (MQ=255)                                                                                  ccAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGt                                                                     >  4:86318/1‑149 (MQ=255)                                                                                        cAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCg                                                               >  5:302482/1‑149 (MQ=255)                                                                                                  ccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGa                                                     >  7:389581/1‑149 (MQ=255)                                                                                                  ccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGa                                                     >  8:257573/1‑149 (MQ=255)                                                                                                  ccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGa                                                     <  3:86318/149‑1 (MQ=255)                                                                                                  ccccATGCAGATCAATGCTCGAGCCTGCCCCGTACAAATTCCTTTCAAGCTCCCAATAG‑CAGCATGACCAGCTGAATTAGGTGGCTTTAAGCCATAGCTCTCACCCACACCAGCCACAAACTCAGCATACGTTTTCTCGATGATACCGa                                                     >  7:206082/1‑149 (MQ=255)                                                                                                   cccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACg                                                                       >  1:271680/1‑130 (MQ=255)                                                                                                   cccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACg                                                                       <  2:271680/130‑1 (MQ=255)                                                                                                       tGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGtt                                                <  7:257573/149‑1 (MQ=255)                                                                                                        gCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTc                                               >  7:363703/1‑149 (MQ=255)                                                                                                        gCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATAG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGACAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTc                                               >  5:380335/1‑149 (MQ=255)                                                                                                            aTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTc                                           >  5:400456/1‑149 (MQ=255)                                                                                                              cAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGc                                         >  3:263301/1‑149 (MQ=255)                                                                                                                aTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGccc                                       >  8:196907/1‑149 (MQ=255)                                                                                                                 tGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGcc                                        <  4:207565/147‑1 (MQ=255)                                                                                                                   cTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTg                                     <  1:261306/148‑1 (MQ=255)                                                                                                                    tCGAGTCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTa                                   >  3:336422/1‑149 (MQ=255)                                                                                                                    tCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCg                                              >  8:153169/1‑138 (MQ=255)                                                                                                                         ccTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTa                              <  4:13532/149‑1 (MQ=255)                                                                                                                             ccccGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                            <  4:336422/147‑1 (MQ=255)                                                                                                                             ccccGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                            <  6:302482/147‑1 (MQ=255)                                                                                                                             ccccGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                            <  4:263301/147‑1 (MQ=255)                                                                                                                                cGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTc                       >  6:130542/1‑149 (MQ=255)                                                                                                                                 gTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTa                              >  7:133366/1‑141 (MQ=255)                                                                                                                                 gTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTa                              <  8:133366/141‑1 (MQ=255)                                                                                                                                 gTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTc                       >  8:302868/1‑148 (MQ=255)                                                                                                                                       atTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAg                >  7:280044/1‑149 (MQ=255)                                                                                                                                       atTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAg                >  8:388493/1‑149 (MQ=255)                                                                                                                                        ttgcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGt               >  2:40376/1‑149 (MQ=255)                                                                                                                                        ttgcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGt               >  5:328996/1‑149 (MQ=255)                                                                                                                                        ttgcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGt               <  5:301919/149‑1 (MQ=255)                                                                                                                                           tggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTggg            >  7:39914/5‑149 (MQ=255)                                                                                                                                            ttgcATGGTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGACAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCAGGTAGTggg            >  1:119427/1‑148 (MQ=255)                                                                                                                                              gcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCgg                                                              >  4:167536/1‑96 (MQ=255)                                                                                                                                              gcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCgg                                                              <  3:167536/96‑1 (MQ=255)                                                                                                                                                   cTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATAcc    >  4:21475/1‑149 (MQ=255)                                                                                                                                                   cTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGGGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATAc     <  7:391063/148‑1 (MQ=255)                                                                                                                                                    tGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGCTAccc   >  1:318974/1‑149 (MQ=255)                                                                                                                                                    tGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATAccc   >  4:311865/1‑149 (MQ=255)                                                                                                                                                     ggCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCg  >  2:26149/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCG  >  NC_000913/1585001‑1585297

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