Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	51.9%	54.3 / 81.1	52	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (22/3);  new base C (0/27);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 8.50e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-08

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993768‑994064                                                                                                                                                         |                                                                                                                                                    gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCcca                                                                                                                                                     >  2:314510/1‑150 (MQ=255) gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCCGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCcca                                                                                                                                                     >  2:1575813/1‑150 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAGCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCcaccac                                                                                                                                                 <  1:1724026/150‑1 (MQ=255)      agccgcgCATCCTGAAACATCATCCGCGTGTCTTCCTGAATTTCACCCAGCGGTGGGGGGCCCGCTAACCC‑‑CGCCGGCGGTTGGCGTTTCCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGc                                                                                                                                              <  1:1606216/147‑1 (MQ=255)               tCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACaaa                                                                                                                                        >  1:157505/1‑149 (MQ=255)                   gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGAc                                                                                                                                   >  1:1428500/1‑150 (MQ=255)                        atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                             >  1:112581/1‑151 (MQ=255)                         tcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                            >  2:180408/1‑151 (MQ=255)                         tcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                            >  2:812827/1‑151 (MQ=255)                           atTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGaa                                                                                                                          >  2:762470/1‑151 (MQ=255)                            ttCGCGTGTCTTCCTGAATTTCACCCGGGGGTGGGGGGCCCGCAAACACATCGCGGGGGGTTGGCTTTCCCAGCCCGCCCAGCAGGCGCAGGGGGGTTTTTTT‑CCCC‑CCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAat                                                                                                                         <  2:1477689/151‑1 (MQ=255)                             tCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAata                                                                                                                        >  1:973007/1‑151 (MQ=255)                             ccgctTGTCTCCCTAAATTCGCGCCGGGGGTGTGGTCCCCCCAACCCCTCCGCCTGTGGTGGGCGTTCCCAGCCCGCCCGGCGGGCCCA‑CCGGGTTCTTTTCCCCC‑CCC‑CCCCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                        <  1:320568/150‑1 (MQ=255)                                cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                     >  1:489984/1‑151 (MQ=255)                                      tcctaaaaTTTAGGCCAGGGGGGGGGGCCCCCCAAACACATCCCCGGGGTTGGGCTTTCCCAGCCCCCCCGGCAGGCGCAGGAGGGTATTTTT‑CCCC‑CCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                               <  2:1350651/145‑1 (MQ=255)                                                    ccgggggTGGGGTCCCCGCAAACCCACCCCCGGGGTTGGGCTTTCCCAGCCCTCCCAGCGGGCCCAGGGGGGTTCTTTT‑CCCC‑CCC‑CCCCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                                 <  2:265211/146‑1 (MQ=255)                                                     aaggggggtggtgCCCGCTACCACCCGGCCTGGGGTGGTATTTTCCAGCCCGGCCGGGAGCCCCAGAAGGGTATTTTT‑CCCA‑CCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                                 <  1:742668/143‑1 (MQ=255)                                                          gtgtGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                           >  2:1326087/1‑151 (MQ=255)                                                          gtgtGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                           >  1:153499/1‑151 (MQ=255)                                                          gtgtGGCCCCCGCAACCCCACCGCCTGCGGTGGGCTTTCCCAGCCCCCCCGGCGGGCGCGGCGGGTTTTTTTT‑GCCC‑CCC‑CCGTTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                            <  1:1132522/150‑1 (MQ=255)                                                                    cGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                 >  1:625688/1‑151 (MQ=255)                                                                      cTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                               >  2:1656363/1‑151 (MQ=255)                                                                      cTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                               >  2:197781/1‑151 (MQ=255)                                                                      cTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                               >  2:1034565/1‑151 (MQ=255)                                                                               cGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                                       >  1:1106276/1‑150 (MQ=255)                                                                                  ctgggtttGGCGTTCCCAGCCCGCCCAGGAGCCGCAGCGGGTTCTTTTT‑CCCC‑CCC‑CCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGAACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                   <  2:649688/145‑1 (MQ=255)                                                                                     gggttggtttccccAGCCCGGCCGGCGGGCGCAGCAGGTTTC‑TTT‑GCCACCCC‑CCGCGGGGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                 <  2:598103/138‑1 (MQ=255)                                                                                        ttggttttcccAGCCCGGCCAGCGGGCGCGGC‑GGGTTCTTTTCCCCC‑CCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                             <  2:1297087/142‑1 (MQ=255)                                                                                         tggctttcacAGCCCTGCCGGACGGCGCAGCGGGGT‑TTTTT‑GACACCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                            <  2:1323243/142‑1 (MQ=255)                                                                                              ttcccaccccTCCCAGCGGGCGCAGCAGGTTTCTTTC‑CCCA‑CCC‑CCTCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                       <  1:812409/143‑1 (MQ=255)                                                                                               ttCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                      >  2:1114368/1‑151 (MQ=255)                                                                                                 ccaccccgcccggcgggCGCAGAAGGGTCCTTTG‑CCCC‑CCC‑CCTGTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                    <  2:1661103/136‑1 (MQ=255)                                                                                                   acaCCGGCAAGCGGGCGCAGCGGGGTTTTTTT‑CCCC‑CCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                  <  2:506060/149‑1 (MQ=255)                                                                                                     cccTGCCAGCGGGCGCAGCAGGTTTTTTTT‑CCCC‑CCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGTCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                  <  2:1395544/148‑1 (MQ=255)                                                                                                                 gcgcAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                    >  2:196784/1‑151 (MQ=255)                                                                                                                             aCTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCCGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                         >  1:229438/1‑150 (MQ=255)                                                                                                                             aCTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                         >  1:142683/1‑150 (MQ=255)                                                                                                                              ttttt‑GCCC‑CCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                        <  2:1258515/150‑1 (MQ=255)                                                                                                                              tctttcccca‑ccc‑CCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                        <  1:491337/144‑1 (MQ=255)                                                                                                                               tctttgcccccca‑ccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                        <  2:1065653/141‑1 (MQ=255)                                                                                                                                    ttcccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  2:567889/142‑1 (MQ=255)                                                                                                                                    tgcccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  2:490157/142‑1 (MQ=255)                                                                                                                                    tcgccccccgccagccGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  1:44618/136‑1 (MQ=255)                                                                                                                                    tccccccccccggggcgcCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  2:640713/137‑1 (MQ=255)                                                                                                                                    tccccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  2:1475595/142‑1 (MQ=255)                                                                                                                                    tccccccccccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  2:1497407/142‑1 (MQ=255)                                                                                                                                    tccccccccccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  2:794306/142‑1 (MQ=255)                                                                                                                                    tccccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  1:947298/142‑1 (MQ=255)                                                                                                                                    ccca‑ccc‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                  <  1:1516846/150‑1 (MQ=255)                                                                                                                                      cccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                   <  2:1307298/143‑1 (MQ=255)                                                                                                                                      cccccccccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                   <  2:353538/143‑1 (MQ=255)                                                                                                                                         cca‑ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg               >  2:1468027/1‑150 (MQ=255)                                                                                                                                               gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc           >  1:658112/1‑149 (MQ=255)                                                                                                                                                  gCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa      >  2:247045/1‑151 (MQ=255)                                                                                                                                                     ccccACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt   <  2:143773/150‑1 (MQ=255)                                                                                                                                                      ccccccccCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATg  <  2:620737/151‑1 (MQ=255)                                                                                                                                                      ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATg  <  1:914102/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                    GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993768‑994064

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