Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,166	0	G	A	40.0%	73.6 / 22.4	56	S312F (TCC→TTC)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base G (25/8);  new base A (0/22);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 4.79e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGT  >  NC_000913/2921016‑2921316                                                                                                                                                        |                                                                                                                                                       ccGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCgg                                                                                                                                                        <  1:427828/151‑1 (MQ=255)  cGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCgg                                                                                                                                                        >  2:145879/1‑150 (MQ=255)                 ccATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTCGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAg                                                                                                                                        <  2:857009/151‑1 (MQ=255)                     aaaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGcaca                                                                                                                                    <  1:512217/151‑1 (MQ=255)                     aaaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGcac                                                                                                                                     <  1:117944/150‑1 (MQ=255)                       aGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAcc                                                                                                                                  >  1:530420/1‑151 (MQ=255)                                  atcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaac                                                                                                                       >  1:429712/1‑151 (MQ=255)                                           acgtggaggtgtttacAGAACACCCTCCCCACCGGCGGTACCATCCCCATTACGTTCGGTTTTGAAGGCGTAGTGTGCGGGGTCCCCGT‑GCCGCGCACCGGCCAGTCTGAAATAATCCCCCCCACCACCCCGCCGATAaaaccgcaaacc                                                                                                               <  2:363731/135‑1 (MQ=255)                                                 aCGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAAccgccg                                                                                                           >  1:1126867/1‑148 (MQ=255)                                                 aCGTAGTGGCAAACCCCCACCACCATCAACAGCAACATGCCCATTCCGATCGTTTTTTTCCCCGCAATGTACAGCGACCCCGT‑GCGGCCCACCACCCAGCCGAAAAAAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCcgga                                                                                                         <  2:154237/150‑1 (MQ=255)                                                    tAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgg                                                                                                     >  1:637990/1‑151 (MQ=255)                                                     aGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgg                                                                                                     >  2:897887/1‑150 (MQ=255)                                                      gTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGCTTTCCGCGCATTGTTCGGCGAGCCCGT‑GCGGCGCATCAGCCAGTCGCAATTAATCCCACCCAGCACACCGCCGATAAGCCCGCACACCGCcggaacgg                                                                                                     >  2:801538/1‑149 (MQ=255)                                                        tGCAGACCCCCACCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGACCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgga                                                                                                    <  1:517303/148‑1 (MQ=255)                                                         gCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGcc                                                                                                >  1:480470/1‑151 (MQ=255)                                                           aGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcac                                                                                              >  2:991052/1‑151 (MQ=255)                                                               caccaccaccatCCACAGCAACATGCCCTTTACGCTGGTTTTTTCACGCGAAATTTACAGAAACCCCGT‑GCGGGTCATCACCCAGTCGAAAAAATCCCCCCCCAAAAACCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaaa                                                                                           <  2:1061450/150‑1 (MQ=255)                                                                     caccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCt                                                                                    >  2:564119/1‑151 (MQ=255)                                                                      accatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCtt                                                                                   >  1:384057/1‑151 (MQ=255)                                                                       ccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCttt                                                                                  >  1:497891/1‑151 (MQ=255)                                                                       ccaccaacacCAACATCCCCATTACGATGGTTGTTTCCCGCGAAATGTCCAGCGACCCCGC‑GCGGCACATCACCCAGCCGAAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCttt                                                                                  <  2:958043/145‑1 (MQ=255)                                                                            gACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGa                                                                              >  1:1062559/1‑150 (MQ=255)                                                                              caacAACATCCCCTTTCCTCTCGTTGTTTTACGCGCATTGTCAAGCTA‑CCCGTGGCGGCCAACCACCCACCCGAAAATAACCCCACCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCa                                                                                <  2:1099201/143‑1 (MQ=255)                                                                              cAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                           >  2:1014748/1‑151 (MQ=255)                                                                                 cAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc                                                                        >  2:587649/1‑151 (MQ=255)                                                                                  aaaaTGCCCATTACGACCGGTGTTTTACGCGAAATGTCCAGCGACCCCGG‑GCGGCGCACCAGCCGGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGa                                                                         <  1:35068/146‑1 (MQ=255)                                                                                  aaCATGCCCATTACGTTCGGTTTTTTCCGCGAAATGTCCAGCGACCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                           <  1:1107558/147‑1 (MQ=255)                                                                                  aaCATGCCCATTACGATCGGTTTTTTCCGCGCAATGTTCAGCGACCCCGT‑GCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCg                                                                          <  1:674935/148‑1 (MQ=255)                                                                                   aCATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACAt                                                                      >  2:410437/1‑151 (MQ=255)                                                                                    cagccccATCACGACGGGTTTTTAACGCGAAATTTCCAGCCACCCCGG‑GCGCCCCACAGGCCAGCCGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACAt                                                                      <  2:89431/146‑1 (MQ=255)                                                                                      tGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATccc                                                                   >  1:117331/1‑151 (MQ=255)                                                                                        cccATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                 >  1:1096575/1‑151 (MQ=255)                                                                                        cccATTACGATCGGTGTTTTACGCGAAATGTCCAGCGACCCCGT‑GCGGCGCATCAGCCAGTGGAAAATAATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                 <  1:669783/151‑1 (MQ=255)                                                                                        cccATTACCACCGGTGTTTTCCGCACAATGTCCAGCGACCCCGT‑GGGGCGCATCACCCAATCGAAAATATCCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcg                                                                  <  1:635858/150‑1 (MQ=255)                                                                                              aCGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCt                                                           >  2:146890/1‑151 (MQ=255)                                                                                               cGATCGGTTTTTTACGCGCAATGTCCAGCGACCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCt                                                           <  1:310582/150‑1 (MQ=255)                                                                                               cGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGCGCCTg                                                          >  2:506140/1‑151 (MQ=255)                                                                                                gctCGGTGTTTTTAGCACAATGCGCACCGCTCCGGT‑GCGGCGCAACAGCCGCTCGAAAAAAATCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCt                                                           <  2:220724/147‑1 (MQ=255)                                                                                                 agcGGTGTTTTCCGCGCAATGTTCAGCGACCCCGT‑GCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCa                                                        <  1:679460/149‑1 (MQ=255)                                                                                                         tttCCGCGAAGTTTTCAGCGACCCCGT‑GCGGCGCACCAGCCAGTCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATa                                                 <  2:389143/150‑1 (MQ=255)                                                                                                          ttACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaa                                                >  2:196677/1‑150 (MQ=255)                                                                                                            aCGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAcc                                             >  2:307424/1‑151 (MQ=255)                                                                                                             cgcgCAATTTCCAGCGACCCCGT‑GCGGCGCATCAGCCAGTGGGAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCg                                            <  1:262732/151‑1 (MQ=255)                                                                                                                    tGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAcc                                      >  2:513595/1‑150 (MQ=255)                                                                                                                     ttcaggCGACCCGGG‑GCGGCCCATCGACCAGTCGAAACCAACCCCACCCACCATGCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCa                                     <  2:497243/145‑1 (MQ=255)                                                                                                                      tccAGCGACCCCGT‑GCGGCCCACAAGCCATCCGAAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAg                                    <  2:629426/148‑1 (MQ=255)                                                                                                                           cgcccccGG‑GCGGCGCATCCCCCAGTCGAAAAAAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGta                                 <  2:530406/144‑1 (MQ=255)                                                                                                                            gATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaata                              >  2:839969/1‑150 (MQ=255)                                                                                                                            cgacccgggggggccCAACAGCAAGTCGAAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAGCCGGGAACCAGGtaat                               <  2:544963/139‑1 (MQ=255)                                                                                                                              ccccctgccggccCACCACCCAGTCGAAAATAACCCCCCCCAACACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGta                                 <  2:637967/138‑1 (MQ=255)                                                                                                                               cccGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaataa                             >  2:1054842/1‑148 (MQ=255)                                                                                                                                 ccaagcggcgGATAAGCGTGTTGAAAATATCCCCCCCCAGCACACCGCCGATAATACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTACCTGCACCAGATAAACCGGGAACCAGGTAATAaagaa                         <  1:858490/147‑1 (MQ=255)                                                                                                                                  gT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaag                        >  2:100423/1‑150 (MQ=255)                                                                                                                                    ggagggacaccAGCAAGCCGAAAAAAACCCCACCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagta                      <  2:1097040/141‑1 (MQ=255)                                                                                                                                          gcATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTc                  >  2:547752/1‑149 (MQ=255)                                                                                                                                              cAGCCAGTCGAAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGc             <  1:934806/150‑1 (MQ=255)                                                                                                                                                 ccAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGtt          <  1:182657/150‑1 (MQ=255)                                                                                                                                                  cggTCGAAAATAACCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaagaag                        <  2:556010/133‑1 (MQ=255)                                                                                                                                                       ggAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATAc     >  1:328493/1‑149 (MQ=255)                                                                                                                                                        gAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt  >  2:62734/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       CCGAAGAAGGCCAGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGT‑GCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGT  >  NC_000913/2921016‑2921316

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