Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,365,454	0	A	C	21.9%	109.7 / 24.9	64	intergenic (‑97/+19)	cutA/dcuA	divalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base A (21/29);  new base C (8/6);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCC  >  NC_000913/4365316‑4365587                                                                                                                                           |                                                                                                                                      caccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGc                                                                                                                                   >  1:560737/1‑143 (MQ=255)   ccgcGACAGACGCGGTATTCGAACTTTGTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGccc                                                                                                                                 <  1:997930/140‑1 (MQ=255)     acGACAGGCTCGGTCCTCGAACTTTTTTCTCCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTAGACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcg                                                                                                                               <  2:322060/143‑1 (MQ=255)     acGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGcaatg                                                                                                                               <  2:745751/143‑5 (MQ=255)      cGACAGACGCGGTATTCGAACTTTTTTCATAAAGCATAATCACATCGCAATGATTCTTTCGTCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgc                                                                                                                              <  2:1004821/143‑1 (MQ=255)        aCAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctt                                                                                                                            <  1:12007/143‑3 (MQ=255)         cAGACGCGGTATTCGAACTTTTTGCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAt                                                                                                                           <  1:803450/143‑1 (MQ=255)             cgcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGAt                                                                                                                       <  2:769181/143‑1 (MQ=255)                 gTATTCGAACTTTTTTCATCAGGCATAAACACCTCGCAAGGACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAtt                                                                                                                   <  1:320060/143‑1 (MQ=255)                  tATTCGAACTTTTTTCATCAAGCACCAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGCACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTa                                                                                                                  <  2:32955/143‑1 (MQ=255)                    ttcgatcatattgaatcaAGCAAAAACACCTAGCAAGAACAGATGGAACGGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACa                                                                                                                <  2:511787/129‑1 (MQ=255)                     tCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGCAAGAAAGCGGGGCGTGCACGCCCCGCGATTGATCATTAccg                                                                                                               >  2:24087/1‑141 (MQ=255)                             gtttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                       <  2:220868/142‑1 (MQ=255)                                   tcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGCTTGATCATTACAGCATGAAGCTACCCa                                                                                                 >  2:649612/1‑143 (MQ=255)                                   tcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTGAGTCGGCGGAAAAGAAAGCGGAGCGTGAACGCCCCGCGATTGACCACTACAGCAACAAGCTACCCa                                                                                                 >  1:327166/1‑143 (MQ=255)                                   tcaAGCATAAACAACTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGCGACCGCCCCGCGACTGATCATTACAGCATGAAGCTACCCa                                                                                                 >  1:364020/1‑143 (MQ=255)                                     aaGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTCAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                                                               >  2:319029/1‑143 (MQ=255)                                      aGCATAAACACCCCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCa                                                                                              <  2:983160/143‑1 (MQ=255)                                       gCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCa                                                                                              <  1:1011187/142‑1 (MQ=255)                                       gCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCa                                                                                              >  2:1011187/1‑142 (MQ=255)                                          tAAACGCCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaa                                                                                          <  2:688560/143‑1 (MQ=255)                                              cacCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGCTTGATCATTACAGCATGAAGCTACCCAGCAcgaccccg                                                                                      >  1:949198/1‑138 (MQ=255)                                                ccgcgcgAGACCCGACGGAACCGAAAGCGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaa                                                                                    <  1:739400/140‑1 (MQ=255)                                                ccTCGCTAGAACAGATGGAACCGCAAGAGAGCGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaa                                                                                    <  1:266535/143‑1 (MQ=255)                                                ccTCGCAAGGACAGATGGAACCTCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaa                                                                                    <  2:101656/143‑1 (MQ=255)                                                 cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                   <  2:755909/143‑1 (MQ=255)                                                   cGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGCTCCTTACAGCATGAAGCCACCCAGCAcgacaccgagcca                                                                                 >  1:545996/1‑139 (MQ=255)                                                    gCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                           <  2:312449/100‑4 (MQ=255)                                                    gCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                           >  1:312449/1‑97 (MQ=255)                                                           cAGATGGAACCGCAAGGGAGAGGTCACTATACAACGGCCGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCa                                                                         <  2:423726/143‑1 (MQ=255)                                                             gagggAACCGCAAGAGAGAGGGCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAgg                                                                       <  2:344317/140‑1 (MQ=255)                                                             gATGGCACCGCAACAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAgg                                                                       <  2:9254/143‑1 (MQ=255)                                                              gtgaaaCCGCCAGGTAGAGGTCACCATAGAAGGGAGGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGCACGCCCCGCGATTGATCCTTACCGCCTGATCCTTCCTAGCACGAAGCCGAAGCAAACGGCCAggg                                                                      <  1:312625/139‑1 (MQ=255)                                                               tGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                     <  2:560737/143‑1 (MQ=255)                                                                      gCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                                                              >  2:385565/1‑143 (MQ=255)                                                                                 gTCACTATACAACGGACGGGGGACGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCAGCGTCCGGTCAGTACAGCATGGAGCTACCCAGCACGCAGCCGAAGCAAACGGCCCGGGCAACACCCAGCGGAGCCgg                                                   >  1:305121/1‑143 (MQ=255)                                                                                 gTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCgg                                                   >  2:253863/1‑143 (MQ=255)                                                                                    aCTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAt                                                >  1:142024/1‑143 (MQ=255)                                                                                     cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATg                                               >  2:572531/1‑143 (MQ=255)                                                                                               gacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                          <  2:989355/58‑1 (MQ=255)                                                                                               gacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                          >  1:989355/1‑58 (MQ=255)                                                                                               gacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGAt                                     <  1:24087/143‑1 (MQ=255)                                                                                               gacggGGGAAGGATGTTAGGCGCCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGAt                                     <  1:851498/143‑1 (MQ=255)                                                                                                    gggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                >  2:935059/1‑143 (MQ=255)                                                                                                     ggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGa                               <  2:194686/143‑1 (MQ=255)                                                                                                          gATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACg                          >  2:303831/1‑143 (MQ=255)                                                                                                           aTGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGGCGAAGCACGGATGGTTGAAGACGa                         >  1:630518/1‑143 (MQ=255)                                                                                                           aTGTCAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGa                         <  1:319029/143‑1 (MQ=255)                                                                                                            tgctagacatcggCAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGaa                        <  1:344928/133‑1 (MQ=255)                                                                                                            tGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGaa                        >  2:834303/1‑143 (MQ=255)                                                                                                            tGTTAGGCTGCGCGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCTGGATGAAGAACGGATGGTTGAAGACGaa                        <  1:410318/143‑1 (MQ=255)                                                                                                               tAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAAttt                     >  2:655129/1‑143 (MQ=255)                                                                                                                aGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTa                    >  1:907730/1‑143 (MQ=255)                                                                                                                 ggcggcggGAAAGAAAGCGGGGCGTGAGCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                   <  2:614776/141‑1 (MQ=255)                                                                                                                 gTCGGCGGGAACGACAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                   <  2:796847/143‑1 (MQ=255)                                                                                                                 gTCGGCGGGAAAGAAGGCGGGGCGTGAACGCCCCGCGATTGATCATGACAGCATGAAGCTCCCCAGCCCGACGCCGTAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                   <  1:763219/143‑1 (MQ=255)                                                                                                                         gaaagaaaGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACgg           <  2:142024/143‑1 (MQ=255)                                                                                                                          aaagaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGACCGGATGGTTGAAGACGAATTTACCGATACggg          >  2:335886/1‑143 (MQ=255)                                                                                                                              atcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCACACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATCTACCGATACAGgtagt      >  2:312701/4‑143 (MQ=255)                                                                                                                               tcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAGGACGAGTTTACCGATACGGgtagta     >  1:685726/3‑143 (MQ=255)                                                                                                                               tcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgtagta     >  2:247889/3‑143 (MQ=255)                                                                                                                               tcgCGGGGCGTTCACGCCCCGCGATGGATGATTACCGCATGAAGCTACCCAGCAAGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGCGAACCCGGGAGGAAGAGCGGACGGTGGAAGGCGCAGATAGGGATACGGggagta     >  2:259784/3‑143 (MQ=255)                                                                                                                                  cGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGCACCCGGGATGAAGAACGGACGGTCGAAGACGAATCTACCGATACGGGCAGTAccc  >  1:587699/1‑143 (MQ=255)                                                                                                                                  cGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTAACCAGAACGAAGAAGAAGCAAACGGGTGGGTCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTAccc  >  2:471252/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                      CACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCC  >  NC_000913/4365316‑4365587

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