Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,531,440	0	A	C	21.0%	131.6 / 10.9	62	Y334D (TAT→GAT)	yhgE	putative transport protein YhgE
Reads supporting (aligned to +/- strand):  ref base A (31/18);  new base C (0/13);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 4.96e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATGCGCATACTGTCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGAACACATACAGCGGAGCAACCAGCAAAGCGGT  >  NC_000913/3531306‑3531578                                                                                                                                       |                                                                                                                                           aaTGCGCATACTGTCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCa                                                                                                                                       <  2:725579/140‑1 (MQ=255) aaTGCGCATACTGTCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCa                                                                                                                                       <  1:273054/140‑1 (MQ=255)   tGCGCATACTGTCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGc                                                                                                                                     <  1:651194/140‑1 (MQ=255)   tGCGCATACTGTCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGc                                                                                                                                     >  1:362929/1‑140 (MQ=255)        aTACTGTCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACca                                                                                                                                >  2:68809/1‑140 (MQ=255)         tACTGTCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACcac                                                                                                                               <  1:262996/140‑1 (MQ=255)          aCTGTCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACcaca                                                                                                                              <  1:317372/140‑1 (MQ=255)          aCTGTCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACcaca                                                                                                                              <  2:119293/140‑1 (MQ=255)              tCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAg                                                                                                                          >  1:418769/1‑140 (MQ=255)                      cGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCa                                                                                                                  >  2:312003/1‑140 (MQ=255)                        gCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGAc                                                                                                                                  >  1:23548/1‑122 (MQ=255)                        gCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGAc                                                                                                                                  <  2:23548/122‑1 (MQ=255)                        gCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGc                                                                                                                >  2:681261/1‑140 (MQ=255)                                  aagCACCAGTATCACCAACCCAAATAACGAAACGGCAGGGTTCATTTGCCCCTAAAAAACAACGGGTTTTGCCCCTTTCCGCCAGACGATGCTGACAAAATACCCCAGCGCCCACACCGTTAACACCAGCACCGCCAGCa                                                                                                      <  2:665695/140‑1 (MQ=255)                                      aCCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCa                                                                                                      <  2:263319/136‑1 (MQ=255)                                      aCCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCa                                                                                                      >  1:263319/1‑136 (MQ=255)                                               accaacaccaaTAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGc                                                                                              >  2:192296/1‑135 (MQ=255)                                               accaacaccaaTAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGc                                                                                              <  1:192296/135‑1 (MQ=255)                                               accaacaccaaTAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCg                                                                                         >  1:480624/1‑140 (MQ=255)                                                    caccaaTAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACTAGCACCGCCAGCACGCCCTGCAAGCGGTCGa                                                                                    <  1:176653/140‑1 (MQ=255)                                                    caccaaTAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGa                                                                                    >  2:509488/1‑140 (MQ=255)                                                     accaaTAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGAc                                                                                   >  1:214232/1‑140 (MQ=255)                                                     accaaTAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGAc                                                                                   >  2:720984/1‑140 (MQ=255)                                                          tAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTcc                                                                              >  1:489781/1‑140 (MQ=255)                                                            aCGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAg                                                                            <  1:720984/140‑1 (MQ=255)                                                               aaaCCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCg                                                                         >  1:616403/1‑140 (MQ=255)                                                               aaaCCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCg                                                                         >  1:217151/1‑140 (MQ=255)                                                               aaaCCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCg                                                                         >  2:145635/1‑140 (MQ=255)                                                                 aCCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAGATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTa                                                                       >  2:212321/1‑140 (MQ=255)                                                                            cACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGCCGATGCTGACAAAATCCCCCAGCGCCCACACCAGTAACCCCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGAccc                                                            <  2:232205/140‑1 (MQ=255)                                                                            cACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGAccc                                                            >  2:81030/1‑140 (MQ=255)                                                                             aCTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCg                                                           >  2:218365/1‑140 (MQ=255)                                                                              cTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGc                                                          <  2:480624/140‑1 (MQ=255)                                                                                    ccTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACCCCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCAt                                                    <  2:418769/140‑1 (MQ=255)                                                                                    ccTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCAt                                                    <  1:751435/140‑1 (MQ=255)                                                                                    ccTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCAt                                                    >  2:260077/1‑140 (MQ=255)                                                                                        aagaaCAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACg                                                >  2:409279/1‑140 (MQ=255)                                                                                          aaacaacGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGcc                                              <  1:654246/139‑1 (MQ=255)                                                                                             caacGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAg                                           >  1:141476/1‑140 (MQ=255)                                                                                             caacGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAg                                           >  1:445451/1‑140 (MQ=255)                                                                                              aacGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAgg                                          <  1:473388/140‑1 (MQ=255)                                                                                                cGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACCCCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAggcg                                        <  1:769491/140‑1 (MQ=255)                                                                                                cGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAggcg                                        >  2:448696/1‑140 (MQ=255)                                                                                                  ggATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAggcggc                                      >  1:123165/1‑140 (MQ=255)                                                                                                  ggATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAggcggc                                      >  1:598190/1‑140 (MQ=255)                                                                                                    tttttGCCCTTTCCGCCAGACGATGCTGACAAAATCCCCCAGCGACCACCCCAGTAACCCCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcga                                    <  2:500694/139‑1 (MQ=255)                                                                                                    aTTTTGCCCTTTACGCCAGACGATGCTGACAAGATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcga                                    <  1:212321/140‑1 (MQ=255)                                                                                                    aTTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcga                                    <  1:133354/140‑1 (MQ=255)                                                                                                    aTTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcga                                    <  2:616403/140‑1 (MQ=255)                                                                                                    aTTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcga                                    <  1:81030/140‑1 (MQ=255)                                                                                                    aTTTTGCCCTTTACGCCAGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcga                                    <  1:681261/140‑1 (MQ=255)                                                                                                    aTTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcga                                    <  1:260077/140‑1 (MQ=255)                                                                                                    aTTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcga                                    <  1:68809/140‑1 (MQ=255)                                                                                                      tttGCCCTTTACGCCAGAGGATGCTGACAAAATCCCCCAGCGACCCCACCAGTAACCCCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcgacg                                  <  2:214232/140‑1 (MQ=255)                                                                                                        tGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGaa                                >  2:755279/1‑140 (MQ=255)                                                                                                          cccTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcgacg                                  >  1:360668/1‑136 (MQ=255)                                                                                                          cccTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGcgacg                                  <  2:360668/136‑1 (MQ=255)                                                                                                           ccTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGAacac                             >  2:191824/1‑140 (MQ=255)                                                                                                                aCGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGAACACATACa                        <  1:312003/140‑1 (MQ=255)                                                                                                                     aGACGATGCTGACAAAATCCCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGAACACATACAGCGGa                   <  1:165565/140‑1 (MQ=255)                                                                                                                      gACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCTCAGGCGGCGACGAACACATACAGCGGAg                  >  2:665292/1‑140 (MQ=255)                                                                                                                       aCGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGAACACATACAGCGGAGc                 <  1:180340/140‑1 (MQ=255)                                                                                                                           tGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGAACACATACAGCGGAGCAAcc             >  2:747578/1‑140 (MQ=255)                                                                                                                                  aaaaTACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGAACACATACAGCGGAGCAACCAGCAAAg      <  1:145635/140‑1 (MQ=255)                                                                                                                                  aaaaTACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGAACACATACAGCGGAGCAACCAGCAAAg      <  1:448696/140‑1 (MQ=255)                                                                                                                                      tACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGAACACATACAGCGGAGCAACCAGCAAAGCGGt  <  1:755279/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                           AATGCGCATACTGTCCAGCACCGGCGAATTAAGAAGCACCAGTATCACCAACACCAATAACGAAACCGCAAGGTTCACTTTGCCCTGAAGAACAACGGGATTTTGCCCTTTACGCCAGACGATGCTGACAAAATACCCCAGCGACCACACCAGTAACACCAGCACCGCCAGCACGCCCTGCAAGCGGTCGACAGTCCAGCCGTACTGAGCGACCCGCAGCCATAACGCCCAGGCGGCGACGAACACATACAGCGGAGCAACCAGCAAAGCGGT  >  NC_000913/3531306‑3531578

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: ‑