Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	46.8%	84.6 / 22.2	62	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (30/3);  new base C (0/29);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.22e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCG  >  NC_000913/907452‑907744                                                                                                                                                       |                                                                                                                                               aCGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTa                                                                                                                                                <  2:155534/151‑1 (MQ=255)  cGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTAc                                                                                                                                               >  2:612141/1‑151 (MQ=255)    tAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCg                                                                                                                                             >  2:1027621/1‑151 (MQ=255)    tAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCg                                                                                                                                             >  2:1035234/1‑151 (MQ=255)    tAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCg                                                                                                                                             >  1:402709/1‑151 (MQ=255)       cgcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTc                                                                                                                                          >  1:620173/1‑151 (MQ=255)         cgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCAc                                                                                                                                        <  1:422410/151‑1 (MQ=255)                atcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCAt                                                                                                                                  >  2:58666/1‑150 (MQ=255)                 tcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATa                                                                                                                                 >  2:1098037/1‑150 (MQ=255)                        gaaattCGTGTTTCAGCCCCTAAATCACCGCCCGGGGGGGGGGGGGCGGTGCGAAAGTAATCACATTAAAAACCCCCCCCAAAACCCAGGGGGGGGGGGGGGGAATGGGGATCACCCGGCGGCTCTTCCCGCTCACCGCCATAAAATGtt                                                                                                                          <  1:935522/145‑1 (MQ=255)                          agATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCa                                                                                                                        >  1:161956/1‑150 (MQ=255)                           gATCGTTTTTCAGCCCCTAAATCACCGCCGTGGGGGGGGGGGGCGGTACGGAAGAAATCACATTGAAAACCCACACCAAAACCCGGGGGGGGGGGGGGGGGAGGGGAATCAACCAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCaaa                                                                                                                      <  2:346190/151‑1 (MQ=255)                              cGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTg                                                                                                                   >  2:651027/1‑151 (MQ=255)                                tgtTCCACCCCCTAAATCAACGCCCGGGGGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAACCCACCCCGAAACCCGGCGGGGGGGGGGGGGGAGGGGGATCAACCAGGGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACt                                                                                                                    <  2:1032956/148‑1 (MQ=255)                                   ttCAGCCCTGAAATCACCGCCCGGGGGGGGGGGGGCGGTAGGGAATAATTCAGATTGAAAACCCACACCGAAATCCAGGGGGTGGGGGGGGGGAGGGGAATCAACCGGCGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGc                                                                                                               <  2:170464/150‑1 (MQ=255)                                   ttCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCt                                                                                                              >  2:913868/1‑151 (MQ=255)                                         cccTGAATCAACGCCCTGGGGGTGGTGGGCGGTACGGAAGTAATCACATGGAGAACCCCCACCGAAACCCGGGGGGTGGGGAGGGGGAGGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATGAAATGTTCAAACTGTTGCTGATAAc                                                                                                        <  2:12098/151‑1 (MQ=255)                                              acccaccccccGGGCGGGGGGGGGGGGTACGAAATTAATCCCATTAAAAACCCCCACCAAAACCCGGGGGGGGGGGGGGGGGATGGGTATCACCCAGGCGCCTTCCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                                      <  1:125647/140‑1 (MQ=255)                                                 cAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAcc                                                                                                >  1:903246/1‑151 (MQ=255)                                                  accGCCCGGGGGGGGGGGGGGGGTAGGGAAGTAATCACATTGAAAACCCCCACCAAAACCCAGGGGGGGGGGGGGGGGAGGGGGATCAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAc                                                                                                 <  2:1142948/147‑1 (MQ=255)                                                  aaCGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAcc                                                                                                >  2:647041/1‑150 (MQ=255)                                                  aaCGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACct                                                                                               >  2:287770/1‑151 (MQ=255)                                                    cGCCGTGGCGGTGGGGGGCGGTACGGAAGTAATCACATGAAGAAACCACCCCGAAAGCCAGGGGGGGGGGGGGGGGAGGGGGATCACCCGGCGGCGTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                             <  2:60944/151‑1 (MQ=255)                                                    cGCCCGGGGGGTGGGGGGGGGTAGGAAAGAAATCACATGAAAAACCCACCCCGAATGCCGGCGGGGGGGGAGGGGGAGGGGGATCAACCAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctc                                                                                              <  2:1132849/150‑1 (MQ=255)                                                    cGCCCGGGGGGGGGTGGGGGGTACGGAAGTAATCACTTTGAAAACCCCCCCCAAAACCCAGGGGGGGGGGGGGGGGAGGGGGATCAACCGGCGGCTCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                             <  2:697349/151‑1 (MQ=255)                                                      cccTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTgg                                                                                           >  1:1114753/1‑151 (MQ=255)                                                      cccGGGGGGTGGTGGGCGGAACGGAAGAAACCACATGAAAAACCCACCCCAAAATCCAGGGGGGGGGGGGGGGGAGGGGAAACAACCGGGGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTgg                                                                                           <  2:530105/151‑1 (MQ=255)                                                          ggCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                                       >  2:636626/1‑151 (MQ=255)                                                                   gggggtaaggaaattAATCACATGGAAAACCCACACCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAACCAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTc                                                                               <  1:899240/136‑1 (MQ=255)                                                                   gggCGGTACGGAAGTAATCACAGTGAAAACCCACCCCAAAATCCAGGGGGTGGGGGGGGGGAGGGGGATCAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCg                                                                              <  1:946176/151‑1 (MQ=255)                                                                        gTACGGAAGTAATCACATTGAGAACCCACCCCGAAACCCGGGGGGGGGGGGGGGGGATGGGGATCAACCAGGGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGAt                                                                         <  2:421372/151‑1 (MQ=255)                                                                         tACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATg                                                                        >  1:174655/1‑151 (MQ=255)                                                                         tACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATg                                                                        >  1:883349/1‑151 (MQ=255)                                                                          agggAAGAAATCAAATTGAAAACCCACCCCAAAACCGGGGGGGGGGGGGGGGGGAGGGGGATCAACCGGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGc                                                                       <  2:899939/149‑1 (MQ=255)                                                                            ggAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCgg                                                                     >  1:1159032/1‑151 (MQ=255)                                                                             gAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgt                                                                    >  1:1083570/1‑151 (MQ=255)                                                                              aaGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCgg                                                                     >  1:689462/1‑149 (MQ=255)                                                                                 tAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                >  2:463887/1‑151 (MQ=255)                                                                                  aaTCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCg                                                               >  2:374479/1‑151 (MQ=255)                                                                                  aaTCAAATGGAAAACCCCCACCAAACCCCGGGGGGGGGGGGGGGGGGGGGGTATCACCCGGGGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                <  1:1118209/150‑1 (MQ=255)                                                                                     cacatgaaaaaaCCACCCCGACCCCCGGGGGGGGGGGGGGGGGGGGGGAATCACCCGGCGGCGCTCCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCt                                                             <  1:325563/141‑1 (MQ=255)                                                                                         ttGAAAACCCACCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGAATCAACCGGGGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTgg                                                           <  2:117809/148‑1 (MQ=255)                                                                                         tggagaACCCCCACCAAAACCCGGGGGGGGGGGGGGGGGATGGGAATCACCCAGCGGCGCTTACCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGcc                                                         <  1:764518/148‑1 (MQ=255)                                                                                         tagagaAACCCCACCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAACCAGGGGCGTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGc                                                          <  1:418564/147‑1 (MQ=255)                                                                                            aaaacccACCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGAACAACCAGGGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTg                                                            <  2:375399/139‑1 (MQ=255)                                                                                              aaaCCACCCCGAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatg                                                    <  2:505954/150‑1 (MQ=255)                                                                                              aaaCCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGCCGCTTACCGCTCACCGCCTTAAAATGTTCAAACTGTTGCTGTTAACTGAGAACCTCTGGCCCGGCGGCTTCGAAGATGCGGGGTTGGCGGGCCGgatg                                                    >  2:470453/1‑150 (MQ=255)                                                                                               aaCCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatgat                                                  >  2:1013445/1‑151 (MQ=255)                                                                                                  cacaCCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCt                                                >  1:555571/1‑150 (MQ=255)                                                                                                      cggAAACCCGGGGGGGGGGGGGGGGGAGGGGGTTCACCCAGGGGCTCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaa                                           <  2:487014/149‑1 (MQ=255)                                                                                                         aaacccgggggggggggggggggAGGGGGATCAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGc                                        <  2:1197055/142‑1 (MQ=255)                                                                                                          accccgggggggggggggggggAGGGGAACCAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                       <  1:19460/143‑1 (MQ=255)                                                                                                                       ggggggggggggggAATCACCCGGCGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgca                          <  1:1032499/151‑1 (MQ=255)                                                                                                                         gAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagc                        >  1:903160/1‑151 (MQ=255)                                                                                                                              ggatgggaatCACCCCGGGGCCCTTCCCTCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt                   <  2:434928/151‑1 (MQ=255)                                                                                                                              ggagggggacAACCCGGCGGCCCTTCCCCCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt                   <  1:328867/151‑1 (MQ=255)                                                                                                                                  gggtacccaCCAGCGGCGCTCCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtttt                 <  2:1098699/141‑1 (MQ=255)                                                                                                                                   gggacCAACCGGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCa              <  2:903291/151‑1 (MQ=255)                                                                                                                                    ggatCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc               >  1:1048646/1‑149 (MQ=255)                                                                                                                                           ccggcggcgCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAAt       <  2:994985/147‑1 (MQ=255)                                                                                                                                             aGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATggg    >  2:524334/1‑151 (MQ=255)                                                                                                                                             aGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATggg    >  1:709546/1‑151 (MQ=255)                                                                                                                                               cggcgCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCg  >  1:964306/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                               ACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCG  >  NC_000913/907452‑907744

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