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%	77.8 / 14.2	60	G197G (GGT→GGG)	frdA	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (22/10);  new base C (0/28);  total (22/38)
Fisher's exact test for biased strand distribution p-value = 7.26e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCG  >  NC_000913/4381580‑4381876                                                                                                                                                      |                                                                                                                                                     cTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGAcc                                                                                                                                                    >  1:821004/1‑151 (MQ=255)   tCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGAccc                                                                                                                                                   <  1:1007486/150‑1 (MQ=255)   tCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGAcc                                                                                                                                                    <  2:360721/149‑1 (MQ=255)        catcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCAcct                                                                                                                                             <  2:300172/151‑2 (MQ=255)        catcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCg                                                                                                                                             >  2:707481/1‑151 (MQ=255)         atcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGc                                                                                                                                            >  2:1019081/1‑151 (MQ=255)             ggATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGt                                                                                                                                        >  2:1176196/1‑151 (MQ=255)                        ccgggCAGCCCGGTGGGGTGATACTGAACGAATTCCATGCCCCGCAGCGGACCGCCGTGGTTTAGCGCCATCCCCAACCCGTCCCCGGTAAC‑GAGGCCGCGGTTGGTGTTTTAACGAAAACCGCGACCCGCACCGCCAGTAGCCATAACGa                                                                                                                             <  2:116004/148‑1 (MQ=255)                         cTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGAc                                                                                                                            >  2:701137/1‑151 (MQ=255)                          gggaagccgggTTGGGTTATACTGAACGAATCCCATGTCACGCGGCGAAACGCGGTGGTTTACCCCCATCCCCATCCCTTCCCGGGTAAC‑GATGCCGCTGTTGGTGTTTTAACAATAACCGCGCCCCGCACCGCCAGTAGCCATAACGAcc                                                                                                                           <  1:664210/142‑1 (MQ=255)                            gCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACcgc                                                                                                                         >  2:1103070/1‑151 (MQ=255)                               gACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGtt                                                                                                                      >  1:752708/1‑151 (MQ=255)                                   ggTTGGGTGATACTGAACGAATTCCATGCCCCGCAGCGGAACGCGGGGGCTTAGCCCCATCCCCATCCCGTCCCGGGAAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                                                                                                  <  1:470734/151‑1 (MQ=255)                                     tggggTGATACTGAACGATTCCCATGTCACGCAGCGGAACGCCGTGGTTTAGCGCCATCCCCACACCGTCACGGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGAAAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCAc                                                                                                                 <  1:243230/148‑1 (MQ=255)                                           gATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTg                                                                                                          >  1:591562/1‑151 (MQ=255)                                            aTACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGc                                                                                                         >  2:711047/1‑151 (MQ=255)                                             tACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGc                                                                                                         >  1:281801/1‑150 (MQ=255)                                                tGACCGATTTCCATGTCACGAAGAGGAACGCCGTGGCTTGGCGCCATACCCATACCGTCACCGGAAAC‑GACCCCGCCGTTGGTGTTGTAAGGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCt                                                                                                           <  1:1106109/145‑1 (MQ=255)                                                  accaattcccATGTCCCGCAGCGGAACGCCGTGGTTTAGCGCAATCCCCATCCCGTC‑CCGGGAACGGAGGCCGCCTTTGTTTTTGAAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGc                                                                                                   <  2:1151409/143‑1 (MQ=255)                                                  aaCGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATCCCCATCCCGCCCCCGGAACC‑GATGCCGCCGTTGGTGTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAcc                                                                                                      <  1:619417/148‑1 (MQ=255)                                                   aCGAATCCCATGTCCCGCGGCGGAACGCCGTGGCTTAGCGCCATCCCCATCCCGTCCCCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCg                                                                                                  <  2:59079/151‑1 (MQ=255)                                                     gAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTg                                                                                                >  2:560156/1‑151 (MQ=255)                                                     aaaTCCCATGCCACGCAGCGGAACGCGGTGGTTTGGCGCCATCCCCATACCGCCACCGGAAAC‑GAGGCCGCCGTGGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTg                                                                                                <  2:750986/150‑1 (MQ=255)                                                      aaTTCCATGTCCCGCAGCGGACCGCGGTGGCTTGGCGCCATCCCCATACCGCCACCGATAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCTCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGc                                                                                               <  1:221549/151‑1 (MQ=255)                                                      aaTCCCATGTAGCGAAGGGGAACGGGGGGTCTTAGCGCCACCCCCAGTCCGTCCCCGGAAAC‑GATCCGGCCGTGGGGGTTGAACCGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGt                                                                                                 <  1:437396/149‑1 (MQ=255)                                                          ccATGTCACGCAGCGGACCGCCGTGGCTTAGCGCCATCCCCATCCCGTCCCGGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTc                                                                                           <  1:848255/151‑1 (MQ=255)                                                          caaTGTCACGAGGCGGACCGCCGTGGTTTAGCGCCATACCCATCCCGCCACCGGTAAC‑GATGCCGCGGTTGGTGTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCtt                                                                                            <  1:687796/148‑1 (MQ=255)                                                           cATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTcc                                                                                          >  1:774236/1‑151 (MQ=255)                                                                cgcGCAGCGGACCGCGGTGGCTTAGCCCAATCCCAAACCCGTCCCCGGAACC‑GATGCCGCCGTTGGTTTTGTAACAATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat                                                                                     <  2:703372/149‑1 (MQ=255)                                                                 aCGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATg                                                                                    >  1:999935/1‑151 (MQ=255)                                                                    cAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTc                                                                                 >  1:482524/1‑151 (MQ=255)                                                                      gCGGAACGCCGTGGTTTGGCGCCATCCCCATACCGTCCCCGGTAAC‑GAGCCCGCCGTGGGTGTTGTAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCa                                                                                <  2:958287/150‑1 (MQ=255)                                                                          aaCGCGGTGGCTTAGCCCCATCCCCATCCCTTCCCGGGTAAC‑GATCCCGCCTTTGTTGTTTTAACAATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCa                                                                                <  1:532624/146‑1 (MQ=255)                                                                             gCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTAcc                                                                        >  2:1060515/1‑151 (MQ=255)                                                                                ggggtttAGCGCCATCCCCATCCCGTCACCGGAAAC‑GATGCCGCCGTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAg                                                                      <  1:1070411/145‑1 (MQ=255)                                                                                gTGGCTAAGCGCAACACCCATTCCGTCACCGTAAAC‑GAGGCGGCGGTTGGTTTTGTAACGAAAACGGCCCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAgg                                                                     <  2:996841/151‑1 (MQ=255)                                                                                 tGGTTTAGCGCCATCCCCACCCCGTCCCCGGAAAC‑GATCCCGCCGTTGGTTTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAg                                                                      <  2:999913/149‑1 (MQ=255)                                                                                  ggCTTAGCGCCATCCCCATCCCGTCCCCGGTAAC‑GATGCCGCCGTTGGTGTTGAAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAg                                                                      <  2:498979/148‑1 (MQ=255)                                                                                  ggCTTAGCCCCATCCCCATCCCGTCCCCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGc                                                                    <  1:1033404/150‑1 (MQ=255)                                                                                  ggCTTAGCCCCACCCCAATCCCGTCCCCGTTACC‑GATGCCGCCGTTGTTGTTTTAACGAAAACCGGGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAg                                                                      <  2:95041/148‑1 (MQ=255)                                                                                    cTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcgc                                                                 >  2:28197/1‑151 (MQ=255)                                                                                     ttggccccatccccatccCGCACCCGGTAAC‑GATGCCCCCGTTGGTGTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcgc                                                                 <  2:858194/144‑1 (MQ=255)                                                                                     tgacccccatccccatccCGCCCCCGGAAAC‑GATGCCGCCGTTGGTGTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcg                                                                  <  2:734151/143‑1 (MQ=255)                                                                                         cgcCATCCCCATCCCGTCACCGGTAAC‑GATGCCGCCTTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAAc                                                             <  1:467842/150‑1 (MQ=255)                                                                                         cgcCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCTCTCATGTTCATTGCTACCAGGCCGCGAAc                                                             >  1:432794/1‑150 (MQ=255)                                                                                         cgcAATACCCACACCGCCCCCGGTAAC‑GTTCCCCCCGTGGGTTTTGTAACGAAAAACGCGCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcgc                                                                 <  2:974245/146‑1 (MQ=255)                                                                                                         ccccgggTAAC‑GATGCCGCTTTTGGTGTTGTAACAAAAAACGCGCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCa                                            <  1:892770/146‑1 (MQ=255)                                                                                                            cgggaaaC‑GATGCCCCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAg                                           <  1:780079/144‑1 (MQ=255)                                                                                                              ggTACC‑GAGGCCGCCGTTGTTGTTGTAACGATAAACGCGACCCCCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAata                                       <  2:820984/151‑1 (MQ=255)                                                                                                              ggTAAC‑GATGCCGCCGTTGGTGTTTTAACAAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAata                                       <  2:785869/151‑1 (MQ=255)                                                                                                               ggaatcaaTGCCCCCGTTGGTGTTGAAACAAAAACCGCAACCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGa                                          <  1:1176492/141‑1 (MQ=255)                                                                                                                 aaC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATc                                     <  1:670163/150‑1 (MQ=255)                                                                                                                             cgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATg                         >  2:274350/1‑151 (MQ=255)                                                                                                                                 ggTGTTGTAACGATAAACGCGCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATg                         <  2:121709/147‑1 (MQ=255)                                                                                                                                     ttgaAACGAAAACCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCa                  <  1:253428/150‑1 (MQ=255)                                                                                                                                          aCGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGc            >  1:229127/1‑151 (MQ=255)                                                                                                                                               aaccgcgCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGAt       <  2:414178/148‑1 (MQ=255)                                                                                                                                               aaaCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGa        >  2:551399/1‑150 (MQ=255)                                                                                                                                               aaaCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGAt       >  2:248230/1‑151 (MQ=255)                                                                                                                                                    cgACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCg  >  1:971910/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAC‑GATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCG  >  NC_000913/4381580‑4381876

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