Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	465,775	G→C	7.3%	A513G (GCC→GGC)	ybaE ←	putative protein YbaE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	465,775	0	G	C	7.3%	222.7 / 5.5	67	A513G (GCC→GGC)	ybaE	putative protein YbaE
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base C (3/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTATATTCCGGTGCTTCGCCAATT  >  NC_000913/465631‑465924                                                                                                                                                 |                                                                                                                                                      aGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTa                                                                                                                                                  <  2:1053990/150‑1 (MQ=255)    cAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAAt                                                                                                                                               <  1:1955369/150‑1 (MQ=255)            ggTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTttc                                                                                                                                       >  1:339285/1‑150 (MQ=255)            ggTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTttc                                                                                                                                       >  1:204845/1‑150 (MQ=255)            ggTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTttc                                                                                                                                       >  1:1422253/1‑150 (MQ=255)               aaaTTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTttcttc                                                                                                                                    <  2:917369/150‑1 (MQ=255)               aaaTTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTttcttc                                                                                                                                    >  2:1493872/1‑150 (MQ=255)                aaTTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTtca                                                                                                                                   <  1:546698/150‑1 (MQ=255)                   tCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGAGAGAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTACATCGGTTTTCTtcatcc                                                                                                                                >  2:1537744/1‑149 (MQ=255)                    cAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAg                                                                                                                               >  1:1720042/1‑150 (MQ=255)                    cAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAg                                                                                                                               >  1:213993/1‑150 (MQ=255)                     aaaCCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAgg                                                                                                                              >  1:365434/1‑150 (MQ=255)                     aaaCCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGCTTTTCTTCATCAgg                                                                                                                              >  1:1536772/1‑150 (MQ=255)                      aaCCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGc                                                                                                                             >  1:1463274/1‑150 (MQ=255)                        ccAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCat                                                                                                                           <  1:755682/150‑1 (MQ=255)                           gCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCattat                                                                                                                        >  2:1272945/1‑150 (MQ=255)                                 cgGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCAc                                                                                                                  <  2:211873/150‑1 (MQ=255)                                 cgGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCAc                                                                                                                  <  1:609033/150‑1 (MQ=255)                                       cAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATc                                                                                                            >  1:1106803/1‑150 (MQ=255)                                        aGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTGAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATcc                                                                                                           <  1:1933788/150‑1 (MQ=255)                                        aGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTAATTGCACGGCATcc                                                                                                           >  1:1807722/1‑150 (MQ=255)                                          tCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAg                                                                                                         <  1:1944077/150‑1 (MQ=255)                                          tCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAg                                                                                                         <  2:1035225/150‑1 (MQ=255)                                              acacCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTc                                                                                                     >  2:1204506/1‑150 (MQ=255)                                              acacCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTc                                                                                                     <  2:1720042/150‑1 (MQ=255)                                               cacCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCg                                                                                                    <  2:1272940/150‑1 (MQ=255)                                                   gTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTg                                                                                                >  1:1193106/1‑150 (MQ=255)                                                      cACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAg                                                                                             >  2:1405261/1‑150 (MQ=255)                                                      cACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAg                                                                                             <  2:1463274/150‑1 (MQ=255)                                                      cACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAg                                                                                             >  1:1467746/1‑150 (MQ=255)                                                      cACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAg                                                                                             <  2:1999021/150‑1 (MQ=255)                                                      cACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCCATTGTAg                                                                                             >  2:1460625/1‑150 (MQ=255)                                                       aCGCCGGGAGTGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGAATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGa                                                                                            <  1:1743381/150‑1 (MQ=255)                                                         gCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGGTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGAtg                                                                                          >  2:44066/1‑150 (MQ=255)                                                             ggAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCg                                                                                      <  2:1593533/150‑1 (MQ=255)                                                              gAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGt                                                                                     <  1:1261140/150‑1 (MQ=255)                                                               aGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTa                                                                                    <  1:964011/150‑1 (MQ=255)                                                               aGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTa                                                                                    >  2:1801455/1‑150 (MQ=255)                                                                    cACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTg                                                                               <  2:1514469/150‑1 (MQ=255)                                                                     aCTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTgg                                                                              <  1:283822/150‑1 (MQ=255)                                                                         aTGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCGGGGCATTAAATCGGGTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCg                                                                          >  2:165876/1‑150 (MQ=255)                                                                                agtgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGaaaa                                                                   >  1:388210/1‑150 (MQ=255)                                                                                agtgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGaaaa                                                                   <  1:1802600/150‑1 (MQ=255)                                                                                agtgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGaaaa                                                                   <  2:204845/150‑1 (MQ=255)                                                                                 gtgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAAc                                                                  >  2:1395295/1‑150 (MQ=255)                                                                                  tgatagtTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACa                                                                 <  1:1493872/150‑1 (MQ=255)                                                                                           aaCAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCcac                                                        <  2:1757527/150‑1 (MQ=255)                                                                                                  gcgTCAGCGTCCCATCTTTCATTAACTGGCTAACTCCCGCTTTCCGGCCTTTAAATCGTTTTTCTTCTTCAGCCATTATTAGCACGGCATCCAGTGTCGAGTGTAGATGGGCGTATGCTGGAGCGTCGACAACATGTGGCCACAGCGGAt                                                 <  2:365434/150‑1 (MQ=255)                                                                                                   cgTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATc                                                <  1:1395295/150‑1 (MQ=255)                                                                                                       aGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATcgcag                                            <  1:1204506/150‑1 (MQ=255)                                                                                                          gTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATAAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATcgcagcgc                                         >  1:1311727/1‑150 (MQ=255)                                                                                                             gCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCgcagcgcagc                                      <  1:1105633/150‑1 (MQ=255)                                                                                                                acTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCAt                                   >  1:911277/2‑150 (MQ=255)                                                                                                                   gTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGc                                >  1:37231/1‑150 (MQ=255)                                                                                                                        tAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCag                           >  2:338016/1‑150 (MQ=255)                                                                                                                             ggCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGtat                      >  1:226216/1‑150 (MQ=255)                                                                                                                             ggCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGtat                      >  1:22088/1‑150 (MQ=255)                                                                                                                              gaagaaaaCCGATTTAATGCCCTGAAAGCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTcc                             >  1:713826/5‑142 (MQ=38)                                                                                                                              gaagaaaaCCGATTTAATGCCCTGAAAGCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTcc                             <  2:713826/138‑1 (MQ=38)                                                                                                                               cTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGGATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTatat                    >  1:360081/1‑150 (MQ=255)                                                                                                                               cTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTatat                    >  1:1696328/1‑150 (MQ=255)                                                                                                                               aagaaaaCCGATTTAATGCCCTGAAAGCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTatat                    >  2:294179/4‑150 (MQ=255)                                                                                                                                 gaaaaCCGATTTAATGCCCTGAAAGCCGTTTTCTTCATCACGCATTATTTGCACGGCATCCAGTGTCGATTCTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCAACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTATATTc                  >  2:740372/2‑150 (MQ=255)                                                                                                                                    aaCCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTATATTCCgg               <  1:1272945/150‑1 (MQ=255)                                                                                                                                      ccGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTATATTCCGGTg             <  2:339285/150‑1 (MQ=255)                                                                                                                                         cTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTATATTCCGGTGCtt          >  1:918981/1‑150 (MQ=255)                                                                                                                                           ttCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTATATTCCGGTGCTTCg        >  1:385548/1‑150 (MQ=255)                                                                                                                                               gggCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTATATTCCGGTGCTTCGCCaa    >  1:50681/1‑150 (MQ=255)                                                                                                                                                 gCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTATATTCCGGTGCTTCGCCAAtt  >  1:246615/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                      AGCCAGGCTTCGGTAAATTCAAACCAGCCGCGCGGTGTCAGTCGCACACCGTTCACGCCGGGAGGGGCACTAATGCGATAGTGATAGTTGAACAGCGGCGTCAGCGTCGCATCTGTCATTAACTGGCTAAAAACCGCTTTCAGGGCATTAAATCGGTTTTCTTCATCAGGCATTATTTGCACGGCATCCAGTGTCGATTGTAGATGTGCGTATGCTGGAGCGTCGAAAACATGTGGCCACAGCGGATCGCAGCGCAGCCATTGCTCCAGAGTATATTCCGGTGCTTCGCCAATT  >  NC_000913/465631‑465924

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