Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,381,728	0	A	C	46.7%	87.7 / 14.8	60	G197G (GGT→GGG)	frdA	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (29/3);  new base C (0/28);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.36e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGG  >  NC_000913/4381589‑4381878                                                                                                                                            |                                                                                                                                                       tcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGcc                                                                                                                                             >  2:87976/1‑151 (MQ=255)  caGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGcc                                                                                                                                             <  2:583370/150‑1 (MQ=255)  caGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGcc                                                                                                                                             >  2:734940/1‑150 (MQ=255)               ccgggCAGACCGGTTGGGTGATCCTGAACGAATCCCATGTCCCGCGGCGAAACGCCGTGGCTTAGCGCAATCCCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTTTAAGAATAACCGCGACCCGCACCGCCAGTAGCCATAACg                                                                                                                                <  1:478375/147‑1 (MQ=255)               ccTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTCGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACg                                                                                                                                >  2:888668/1‑150 (MQ=255)                      gcccGGTTGGGTGATACTGAACGAATCCCATGCCACGCGGCGGAACGCCGGGGCTTAGCGCCATCCCCATCCCGTCCCCGGTACCGATGCCGCCGTGGGTGTTGTAACAAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGt                                                                                                                         <  1:123708/148‑1 (MQ=255)                        ccGTTTGGGTAATCCTGACCGAATCCCAGGCCCCGCGGCGAACCGCCGTGGCTTAGCGCCATACCCATCCCGCCCCCGTAACCGAGGCCGCCGTTGGTGTGGAAACGATAAACGCGCCCCGCTCCGCCAGTAGCCATAACGACcgc                                                                                                                           <  1:876537/146‑1 (MQ=255)                        ccGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGtt                                                                                                                        >  2:573494/1‑149 (MQ=255)                        ccGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGtt                                                                                                                        >  1:833947/1‑149 (MQ=255)                          ggTTGGGTGATACTGAACGAATCCCATGTCACGCAGCGGACCGCCGTGGCTTAGCGCCATCCCCATCCCGTCCCCGGTAACGATCCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                                                                                                    <  1:967367/151‑1 (MQ=255)                           gtggggTGATACTGAACAATTTCCATTCCACGCGGCGGACCGCGGGGCTTTAGCGCCATCCCCACCCCGCCCCCGGTAACGATGCCGCCGTTGGTGTTGTAAGGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACcgc                                                                                                                           <  1:403309/140‑1 (MQ=255)                              gggTGATACTGAACAAATTCCATTTCCCGCAGGGGAACGCGGGGGTTTGGCCCCATACCCATCCCGTCCCCGGAAACGATCCCGCCGTGGGTTTTGTAACAAAAAACGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACg                                                                                                                  <  2:609279/149‑1 (MQ=255)                                      cTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAcc                                                                                                        >  1:567920/1‑151 (MQ=255)                                       tGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAcc                                                                                                        >  1:358048/1‑150 (MQ=255)                                        gAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAcc                                                                                                        >  2:1353/1‑149 (MQ=255)                                         aaaaatgtCCGTGCCACGCAGGGAAACGCCGTGCCTTAGCCCCATCCCCACCCCGCCCCCGATACCGATGCCGCCGTTGGTGTTTTAACAAGAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAcc                                                                                                        <  1:401282/141‑1 (MQ=255)                                           cGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGt                                                                                                   >  1:477951/1‑151 (MQ=255)                                           cGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCAGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCg                                                                                                    >  1:764377/1‑150 (MQ=255)                                                 ccATGTCACGCAGCGGACCGCCGTGGTTTAGCGCCATCCCCATCCCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCtt                                                                                              <  2:186290/150‑1 (MQ=255)                                                         cGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATCCCGTCCCCGGAAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATg                                                                                      <  2:949696/150‑1 (MQ=255)                                                          gCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt                                                                                    >  1:469747/1‑151 (MQ=255)                                                             gCGGAACGCGGTGGCTTAGCGCCATCCCCATCCCGTCCCCGGTAACGATGCCGCCGTTGGTGTTGTAACGAAAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt                                                                                    <  2:375883/148‑1 (MQ=255)                                                             gCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCAt                                                                                 >  1:616332/1‑151 (MQ=255)                                                                 acggcgggggtttACCGCCATACCCATCCCGTCCCCGGAACGGAGGCCGCGGTTGGTGTTGTAACGAAAAACGGGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGc                                                                              <  1:427476/139‑1 (MQ=255)                                                                       gTGGTTTAGCCCCATACCCATCCCGTCACCGATAACGATGCCGCCGTTGGTTTTGTAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTAcc                                                                          <  2:525990/148‑1 (MQ=255)                                                                        tGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGc                                                                      >  1:305600/1‑151 (MQ=255)                                                                           cTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcgc                                                                   >  1:985449/1‑151 (MQ=255)                                                                           cTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcgc                                                                   >  2:220909/1‑151 (MQ=255)                                                                            ttAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCAGAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCACCACGTTCATTGCTACCAGGCcgcg                                                                  >  1:498937/1‑151 (MQ=255)                                                                             tAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcgcg                                                                  >  1:215438/1‑150 (MQ=255)                                                                              aGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGa                                                                 >  1:546829/1‑150 (MQ=255)                                                                                cccctcccccaccccGTCCCCGGAAACGGTCCCCCTGTGGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACa                                                              <  1:399642/138‑1 (MQ=255)                                                                                       cccatcccGTCACCGGTAACGAGGCCGCGGTTGGTGTTGTAACAATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcat                                                       <  2:928835/151‑1 (MQ=255)                                                                                       ccaatcccgcccccGGTACCGAGCCCGCCGTGGGTGTTGAAAGAAAAACCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACa                                                              <  1:126390/132‑1 (MQ=255)                                                                                        ccaccccGTCCCCGGAACCGAGGCCGCCGTTGGTGTTTTAACGATAACCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcatc                                                      <  1:781857/146‑1 (MQ=255)                                                                                          ataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcatcat                                                    >  2:236422/1‑151 (MQ=255)                                                                                             ccGTCCCGGGTAACGATGCCGCCGTTGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCAtcatca                                                  <  2:492224/150‑1 (MQ=255)                                                                                             ccGTCCCCGGTAACGATGCCGCCTTGGGTGTTGTACCGATAAGCCCACCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCaa                                                 <  1:734634/151‑1 (MQ=255)                                                                                                   cgggTAACGTTGCCCCCTTGGGGGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAg                                             <  1:160197/147‑1 (MQ=255)                                                                                                   ccGTTAACGGTCCCGCCGTGGGTGTTGAACCGAAAACCGCACCCCCCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGa                                            <  1:760146/150‑1 (MQ=255)                                                                                                        aaCGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATcc                                      >  2:473059/1‑151 (MQ=255)                                                                                                         aCGAGGCCGCCGTTGGTGTTGTAACGATAACCGCGCCCCGCGCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCa                                     <  1:725020/151‑1 (MQ=255)                                                                                                           gATGCCGCCGTTGTTGTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAg                                    <  2:987404/150‑1 (MQ=255)                                                                                                            aTGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGc                                   >  2:971973/1‑150 (MQ=255)                                                                                                               ccgcctTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCa                                  <  2:666946/148‑1 (MQ=255)                                                                                                                   cgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATg                           >  1:634499/1‑151 (MQ=255)                                                                                                                             gtAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCaaa                  >  1:966505/1‑150 (MQ=255)                                                                                                                             gtAACGAAAACCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCaaaa                 <  1:888321/151‑1 (MQ=255)                                                                                                                              tAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAAc                >  1:646341/1‑151 (MQ=255)                                                                                                                               aactaaaaaCGCGCCCCCCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCaaaa                 <  2:175647/143‑1 (MQ=255)                                                                                                                               aaCGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCaaaa                 <  1:644908/149‑1 (MQ=255)                                                                                                                               aaCAATAAACGCCCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAAc                <  2:35372/150‑1 (MQ=255)                                                                                                                                aCGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGc              >  2:479212/1‑151 (MQ=255)                                                                                                                                 cGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCt             >  1:881994/1‑151 (MQ=255)                                                                                                                                   aaaaccgcgCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTgg           <  2:710178/146‑1 (MQ=255)                                                                                                                                   aaaaaCCCGCCCGGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTg            <  2:591100/148‑1 (MQ=255)                                                                                                                                     aaaCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTg            <  2:370618/148‑1 (MQ=255)                                                                                                                                     aaaCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACg               >  1:308589/1‑145 (MQ=255)                                                                                                                                     aaaCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGa          >  2:296903/1‑150 (MQ=255)                                                                                                                                       aCGCGACCCGCACCGCCAGAAGCCACAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCt       >  2:855175/1‑151 (MQ=255)                                                                                                                                            aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCggg  >  1:34197/1‑151 (MQ=255)                                                                                                                                            |                                                                                                                                                       TCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGG  >  NC_000913/4381589‑4381878

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