Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0104734,658,8880AG38.5% 56.2 / 40.0 39intergenic (+143/+277)yjjV/yjjWmetal‑dependent hydrolase/YjjW family glycine radical enzyme activase
Reads supporting (aligned to +/- strand):  ref base A (16/8);  new base G (5/10);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 5.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCG  >  NC_010473/4658688‑4659086
                                                                                                                                                                                                                  |                                                                                                                                                                                                                          
ggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGCGCCGTAGGCCGGATAAGGCGTTGACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑G‑‑‑‑CGA‑CCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGGTCG‑GCGCACAAt                                                                                                                                                                                                                          <  1:8188S7/200‑1 (MQ=21)
 aagtagaTTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGcaca                                                                                                                                                                                                                            >  1:7900S43/5‑200 (MQ=21)
  aTGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACaa                                                                                                                                                                                                                           <  1:15666S7/200‑1 (MQ=34)
            gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGc                                                                                                                                                                                                                >  1:323S38/1‑200 (MQ=33)
                    ttcTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTGACGC‑‑ATTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTc                                                                                                                                                                                                     >  1:8761S59/3‑200 (MQ=17)
                                   acaaattctcgcgccacccGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTCGT‑‑GCG‑CA‑GATGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑AACGCT‑GACGCGTC‑‑‑‑‑‑‑TTATCagac                                                                                                                                                                                >  1:2475S21/18‑197 (MQ=2)
                                          gttTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CAGACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACaa                                                                                                                                                                          <  1:764S90/200‑1 (MQ=17)
                                           tttAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑g                                                                                                                                                                        >  1:24229S36/1‑200 (MQ=32)
                                            ttAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGGCATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑g                                                                                                                                                                        >  1:19333S16/1‑200 (MQ=17)
                                                  acccGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTGGCATGCTTCGCAGTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACaaattct                                                                                                                                                                     >  1:31922S16/3‑195 (MQ=2)
                                                                         cACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGT‑C‑GAACCGTAGGCCGGATAA‑GGCg                                                                                                                                    >  1:22316S5/1‑200 (MQ=17)
                                                                          aCGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTACGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAAAGT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCg                                                                                                                                    >  1:12584S25/1‑200 (MQ=17)
                                                                              cg‑cATCCGGCAGTCGT‑‑GCG‑CA‑GATGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑TGC‑GTGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCAcg                                                                                                                              <  1:32027S6/200‑1 (MQ=2)
                                                                                       tccggcagtcgtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTA‑AA‑TT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCa                                                                                                                         <  1:1959S76/182‑1 (MQ=2)
                                                                                           ttG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGAATGCCGACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGCGCTAAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATc                                                                                                                       <  1:32573S10/200‑1 (MQ=17)
                                                                                             g‑‑‑gCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAg                                                                                                                 >  1:26551S85/1‑200 (MQ=32)
                                                                                             agtcgtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCg                                                                                                                     <  1:22890S9/188‑1 (MQ=2)
                                                                                               tcgtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGc                                                                                                                   >  1:280S3/11‑200 (MQ=2)
                                                                                               tcgtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGc                                                                                                                   <  1:24231S32/190‑1 (MQ=2)
                                                                                                 gtgcgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAg                                                                                                                 <  1:3180S58/192‑1 (MQ=2)
                                                                                                 gcg‑cACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGt                                                                                                                >  1:6875S15/1‑200 (MQ=32)
                                                                                                    cgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑g                                                                                                            <  1:30953S3/195‑1 (MQ=2)
                                                                                                    cgcagaTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GACGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑g                                                                                                            <  1:28620S14/195‑1 (MQ=2)
                                                                                                      acaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑gc                                                                                                         >  1:25996S7/1‑200 (MQ=32)
                                                                                                         aTGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑ATCCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑gtgcg                                                                                                        >  1:285S11/1‑196 (MQ=2)
                                                                                                          tGCCTGATGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAA‑AT‑‑TT‑ACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑gtgctggga                                                                                                    <  1:27323S24/200‑9 (MQ=2)
                                                                                                              tGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATg                                                                                                >  1:17914S23/1‑200 (MQ=32)
                                                                                                                aTGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATTCTCGCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑GTGCGC‑agatgcc                                                                                               <  1:496S2/200‑7 (MQ=2)
                                                                                                                 tGCGACGCT‑AACGCGTCTTATCAGGCCTACAAATT‑C‑CCGCGCCATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGACCTACAAATT‑‑CTCACGC‑C‑‑ATCCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑gtgctgggatgggc                                                                                               >  1:10160S18/1‑187 (MQ=0)
                                                                                                                            aCGCGTCTTATCAGGCCTACAAATT‑C‑GCGC‑C‑ACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TCGTGCGCA‑GATGCCTGA‑TGC‑GACGCT‑AACGCGTC‑‑‑‑‑‑‑TTATCAGGCCTACAA‑ATTCCC‑GCGC‑C‑‑ATCCGTAGGCCGGATAAAGGCGTTCACGCCGCATCCGGCAGTC‑‑GTGCGC‑A‑GATGCCTGATGCGACGCt                                                                                  <  1:19201S149/200‑1 (MQ=2)
                                                                                                                             cgcgTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCtt                                                                                 >  1:8183S8/1‑200 (MQ=32)
                                                                                                                                           ccTACAGGTT‑T‑GTGC‑CGGACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTATCa                                                                   <  1:15240S23/200‑1 (MQ=17)
                                                                                                                                            cTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCCTATCAt                                                                  >  1:11331S4/1‑200 (MQ=17)
                                                                                                                                              aCAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGc                                                                <  1:8104S94/200‑1 (MQ=32)
                                                                                                                                               cAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGcc                                                               <  1:27075S58/200‑1 (MQ=32)
                                                                                                                                                              aaccgttgggCCGGATAAGGCATTAGCGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTCTTGTGGTTTATTATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTCCGG‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGTCTAc                                                            >  1:4332S65/9‑200 (MQ=14)
                                                                                                                                                                            aTAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGcc                                     <  1:4727S8/200‑1 (MQ=32)
                                                                                                                                                                                 gCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATa                                >  1:11548S43/1‑200 (MQ=32)
                                                                                                                                                                                   gTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑AAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg                              >  1:21091S46/1‑200 (MQ=17)
                                                                                                                                                                                       aCGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGt                          >  1:9452S39/1‑200 (MQ=32)
                                                                                                                                                                                                            gcgcACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAACGCCCGACGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGt       <  1:32720S16/200‑1 (MQ=17)
                                                                                                                                                                                                                  aaTGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑A‑GGGCGTTAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGgcg  >  1:26832S39/1‑200 (MQ=17)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                          
GGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCG‑CATCCGGCAGTTG‑‑‑GCG‑CACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTT‑T‑GTGC‑CGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAG‑TTG‑GCGCACAATGCCTGA‑TGC‑GACGCTTGTCGCGTC‑‑‑‑‑‑‑TTATCATGCCTACAA‑GT‑‑CT‑GTGC‑C‑GAACCGTAGGCCGGATAA‑GGCGTTCACGCCGCATCCGGCAGTC‑‑G‑GCGC‑ATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCG  >  NC_010473/4658688‑4659086

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 49 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑