Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,564	0	A	C	28.3%	96.6 / 27.2	53	*559E (TAA→GAA)	ybaL	putative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base A (34/4);  new base C (0/15);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 6.20e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGG  >  NC_000913/501426‑501700                                                                                                                                           |                                                                                                                                          ttAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑c                                                                                                                                         <  1:309947/140‑1 (MQ=255)  tAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑cc                                                                                                                                        <  2:697027/140‑1 (MQ=255)           gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTc                                                                                                                                   >  1:656633/1‑136 (MQ=255)           gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTc                                                                                                                                   <  2:656633/136‑1 (MQ=255)           gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGC‑GGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTcaccac                                                                                                                              >  2:994159/1‑140 (MQ=255)              ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCt                                                                                                                            >  2:1596992/1‑140 (MQ=255)              ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCt                                                                                                                            >  1:1797990/1‑140 (MQ=255)                ttCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCa                                                                                                                          >  2:738763/1‑140 (MQ=255)                         cAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGt                                                                                                                 >  2:661570/1‑140 (MQ=255)                          aaaaaaGGCGTCCCCCCCGCATCTGGCTTTTGCCTCTCGATCCCTGATGCGCCGCGGGCGCTTCTTATCATCCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTg                                                                                                                <  2:325165/136‑1 (MQ=255)                            aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGc                                                                                                              >  1:1934680/1‑140 (MQ=255)                            aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCGTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGc                                                                                                              >  2:1047495/1‑140 (MQ=255)                                 gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                         >  1:1446129/1‑140 (MQ=255)                                 gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                         >  2:1117866/1‑140 (MQ=255)                                 gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                         >  1:325165/1‑140 (MQ=255)                                   tttCCCGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCGGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTC‑CCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCca                                                                                                       <  2:1373073/139‑1 (MQ=37)                                    ttccccccgcCTTTGGCTTTTGGCCCCCGATTCCTGAGGCGGCGCGGGCCCCTTTTTTCATCCCTCAATTTTTTTCATTTTTTACCCCCCGCAACCCCCGTTTC‑CCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCcag                                                                                                      <  2:687168/134‑1 (MQ=255)                                    ttCACGCCGCATCTGGCATTTGGCTCTTGATGCCTGGTGCGCCGCTGGGGCGTCTTTTGGTGCCTCCTTTTTTTTGATATTTTCCATCGGCCACCCCCCGTTTC‑CCCCATCCCCCCCTCTCCCGCCGGTGGCGTTTCcag                                                                                                      <  2:646579/140‑1 (MQ=255)                                    ttCACCCCGCTTCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCGGGCGCGTTTTACCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTC‑CCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCcag                                                                                                      <  2:888997/140‑1 (MQ=255)                                    tccACGCCGCTTCTGGTTTTTGGCTCTCCATCCCTGATGCGCCGCGGGCCCGTTTTATCAGCCCTACATTTTTTTCATTTTTTACCTCCGCCACCCCCTGTTTA‑CCCCGTCACCCCCTCAACCGCCGGTGGCGTTTCcag                                                                                                      <  2:8955/138‑1 (MQ=255)                                     tCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                     >  2:1230180/1‑140 (MQ=255)                                      cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                    >  1:1268218/1‑140 (MQ=255)                                        cgccgcATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACATCCGGCACCCCCCTTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGt                                                                                                  <  2:1116823/140‑1 (MQ=255)                                         cccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCGTTTC‑CCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGtt                                                                                                 <  1:1596992/139‑1 (MQ=255)                                            gcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCa                                                                                              >  1:1110386/1‑140 (MQ=255)                                             cttCTGGCTTTTGGCTCTCGATGCCTGATGCGCCGCGGGCGCGTCTTATCATGCCTACTTTTTTTTCATATTTTACCCCCGGCACCCACCTTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                             <  2:1596777/138‑1 (MQ=255)                                                   gCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                       >  1:711280/1‑140 (MQ=255)                                                       ttGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgg                                                                                   >  2:1668582/1‑140 (MQ=255)                                                              tcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                            >  1:306124/1‑140 (MQ=255)                                                                   gCCTGCGGCGACGCTGGCGCTTCTTATCATGCCTACATTTTTTTCATATTTTACCCCCGCCAACCCCCGTTTC‑CCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                       <  2:177676/140‑1 (MQ=255)                                                                       gATGCGACGCTGGCGCTTTTTATCATCCCTACATTTTTTTCATATTTTACACCCGCCACCCCCCGTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                   <  2:164546/140‑1 (MQ=255)                                                                       gATGCGACGCTGGCGCGTCGTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAACCACCGTTTC‑CCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                   <  1:1047495/140‑1 (MQ=255)                                                                        aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCACCGTTTC‑CCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                                  <  1:1644809/140‑1 (MQ=255)                                                                          gCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTC‑CCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCa                                                                <  1:1115407/140‑1 (MQ=255)                                                                           cgccGCTGGCGCTTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                               <  1:1668582/137‑1 (MQ=255)                                                                            gACGCTGGCGCTTTTTATCATCCCTACATTTTTTTAATTTTTTACACCCGGCACCCCCCTTTTC‑CCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                              <  2:306124/140‑1 (MQ=255)                                                                             aCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                             >  1:1876030/1‑140 (MQ=255)                                                                              cGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                            >  2:191217/1‑140 (MQ=255)                                                                                              ctgccttACATTTTTTTCTTTTTTTCCTTCCGGCAACCTCCGTTTC‑CCCCGCCCCCCCCCATCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacca                                            <  2:1934680/134‑1 (MQ=255)                                                                                                tGCCTACATTTTTTTCATTTTTTCCACCCGGCACCCCCCGTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg                                          <  2:1268218/140‑1 (MQ=255)                                                                                                 gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                         >  1:1152914/1‑140 (MQ=255)                                                                                                    tACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTC‑CCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                      <  1:661570/140‑1 (MQ=255)                                                                                                    tACATTTTTTTCATTTTTTACACCCGGCAACCACCGTTTA‑CCCCGCCACCCCCACACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                      <  2:1180499/140‑1 (MQ=255)                                                                                                    tACATTTTTTTCATATTTTACATCCGGCAACCACCGTTTC‑CCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                      <  1:55137/140‑1 (MQ=255)                                                                                                     aCATTTTTTTCATTTTTTACATCCGGCAACCCCCGTTTC‑CCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                     <  1:877066/140‑1 (MQ=255)                                                                                                     aCATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                     >  1:1205486/1‑140 (MQ=255)                                                                                                            tttCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                              >  2:1889450/1‑140 (MQ=255)                                                                                                             ttCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                             >  2:1282535/1‑140 (MQ=255)                                                                                                             ttCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                             >  1:1155330/1‑140 (MQ=255)                                                                                                             ttCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                             >  1:515656/1‑140 (MQ=255)                                                                                                             ttCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                             >  2:1475529/1‑140 (MQ=255)                                                                                                               cATTTTTTACACCCGGCAACCCCCGTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                           <  2:547009/140‑1 (MQ=255)                                                                                                               cATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                           >  2:616055/1‑140 (MQ=255)                                                                                                                    tttACTCCCGGCAACCCCCTTTTC‑CCCCTCCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                      <  2:1876030/140‑1 (MQ=255)                                                                                                                      tACACCCGGCACCCCCCGTTTC‑CCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                    <  2:1110386/140‑1 (MQ=255)                                                                                                                       taccccGGCCACCCCCCTTTCCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                    <  2:415195/136‑1 (MQ=255)                                                                                                                        cctCCTGCACCCACCATTCC‑CCCCGCCCCCCCCCCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATc                  <  2:1742264/138‑1 (MQ=255)                                                                                                                           ccGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc               >  2:1344323/1‑140 (MQ=255)                                                                                                                             ggCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt             >  2:1530037/1‑140 (MQ=255)                                                                                                                              gCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGTGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa            >  2:944663/1‑140 (MQ=255)                                                                                                                              gCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa            >  2:1103928/1‑140 (MQ=255)                                                                                                                                     ccGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg     >  1:1861003/1‑140 (MQ=255)                                                                                                                                      cGTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc    <  1:1746751/140‑1 (MQ=255)                                                                                                                                        tttA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCATGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                    >  2:531222/1‑122 (MQ=255)                                                                                                                                        tttA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCATGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                    <  1:531222/122‑1 (MQ=255)                                                                                                                                        tttA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg  >  1:1177556/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                          TTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGG  >  NC_000913/501426‑501700

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