Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,161,466	G→A	G115S (GGC→AGC)	fabR →	transcriptional repressor of fabA and fabB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,466	0	G	A	100.0%	257.2 / NA	68	G115S (GGC→AGC)	fabR	transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/36);  total (32/36)

GAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCTTCAG  >  NC_000913/4161318‑4161611                                                                                                                                                     |                                                                                                                                                  gagCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGc                                                                                                                                                 <  1:477797/151‑1 (MQ=255)    cGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCAcc                                                                                                                                              >  1:120562/1‑151 (MQ=255)        ttAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTcc                                                                                                                                           <  1:270816/150‑1 (MQ=255)             gCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCg                                                                                                                                     <  1:133051/151‑1 (MQ=255)                aCGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGttt                                                                                                                                  >  1:1059292/1‑151 (MQ=255)                         cATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCCCATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTgccgcc                                                                                                                         >  2:569197/1‑151 (MQ=255)                         cATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTgccgcc                                                                                                                         >  1:161607/1‑151 (MQ=255)                             cgcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTg                                                                                                                     <  2:253391/151‑1 (MQ=255)                              gcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTgc                                                                                                                    <  2:120560/151‑1 (MQ=255)                               ccAGGCGCGTCCGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTGCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGCTTCGTGCCGCCGTTgcg                                                                                                                   >  2:154573/1‑151 (MQ=255)                                cAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTgcgc                                                                                                                  >  1:172888/1‑151 (MQ=255)                                cAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTgcgc                                                                                                                  <  1:244723/151‑1 (MQ=255)                                   gcgcGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTg                                                                                                               <  1:907072/151‑1 (MQ=255)                                   gcgcGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTg                                                                                                               <  2:936141/151‑1 (MQ=255)                                       gTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAAt                                                                                                           >  1:602399/1‑151 (MQ=255)                                               aTCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACt                                                                                                   <  1:800307/151‑1 (MQ=255)                                                    ccaaGGCGGGAGTGTGATCCGCCCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATCATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCGCCTCCGCTGCGTTTCGTGCCCCCGTTGCGCGTGAAATTCAGCACTTCa                                                                                                <  2:34545/147‑1 (MQ=255)                                                    cAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCAtt                                                                                              >  1:226210/1‑151 (MQ=255)                                                      aaGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTg                                                                                             >  2:952303/1‑150 (MQ=255)                                                        ggCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCgg                                                                                          >  1:864215/1‑151 (MQ=255)                                                          tgggAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGAAATAATCCTAACGCCTTCCGGTTATTATTGAGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGa                                                                                         >  2:705902/2‑150 (MQ=255)                                                          cGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGaa                                                                                        <  2:1022544/151‑1 (MQ=255)                                                           gggAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAAc                                                                                       >  2:38373/1‑151 (MQ=255)                                                            ggAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACt                                                                                      <  2:304018/151‑1 (MQ=255)                                                              aGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTg                                                                                    <  1:248087/151‑1 (MQ=255)                                                               gtgtGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGc                                                                                   <  2:370606/151‑1 (MQ=255)                                                                 gtgATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCgg                                                                                 >  2:821111/1‑151 (MQ=255)                                                                 gtgATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCgg                                                                                 >  1:36210/1‑151 (MQ=255)                                                                      ccGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTAt                                                                            >  1:1050415/1‑151 (MQ=255)                                                                         cACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCt                                                                          <  2:513647/150‑1 (MQ=255)                                                                           ccTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTgg                                                                        >  1:433422/1‑150 (MQ=255)                                                                            cTCGGTCTGCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGCTTATTATCGCGGGACCGCTCAAGCACCTACGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCTGACCTTGCGGACTATCTgg                                                                        <  2:337649/149‑1 (MQ=255)                                                                             tCGGTCTACACATTTATGGAGTTAATAAGTAATAATCCTAACGCCTTCCGGGTATTATTGCGGGAACTCTCGAGCACCTCCGCTACGTATCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCAGAACTTGCGGATTATCTg                                                                         <  2:603931/147‑1 (MQ=255)                                                                                  ctcCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGa                                                                 >  1:521730/1‑150 (MQ=255)                                                                                    ccACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGaaaa                                                              <  2:148036/151‑1 (MQ=255)                                                                                    ccACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGa                                                                 <  1:102609/148‑1 (MQ=255)                                                                                          ttATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCatat                                                        >  1:470465/1‑151 (MQ=255)                                                                                          ttATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCata                                                         <  1:648506/150‑1 (MQ=255)                                                                                            aTGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGc                                                      <  1:323300/151‑1 (MQ=255)                                                                                              ggAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCcg                                                    >  2:988287/1‑151 (MQ=255)                                                                                                  ttCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCgcgtg                                                >  2:320768/1‑151 (MQ=255)                                                                                                     aTCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTTCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCgcgtgcgt                                             <  2:509985/151‑1 (MQ=255)                                                                                                       cGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGttt                                           >  2:399819/1‑151 (MQ=255)                                                                                                            ataatCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTg                                       >  1:759860/1‑150 (MQ=255)                                                                                                              aatCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGaa                                     <  1:933395/150‑1 (MQ=255)                                                                                                              aatCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGa                                      >  2:762780/1‑149 (MQ=255)                                                                                                              aatCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGa                                      >  1:659733/1‑149 (MQ=255)                                                                                                                tCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAgcg                                  >  2:1941/1‑151 (MQ=255)                                                                                                                tCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGATTCGTGCCGCCGTTGCGCGGGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGGTTACTGAAgcg                                  >  2:293645/1‑151 (MQ=255)                                                                                                                   tAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCa                                >  1:543256/1‑150 (MQ=255)                                                                                                                      cGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGcc                            <  1:678037/151‑1 (MQ=255)                                                                                                                         cTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGa                          >  1:574525/1‑150 (MQ=255)                                                                                                                           tCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGaa                         <  2:261472/149‑1 (MQ=255)                                                                                                                            ccGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGc                       <  1:444344/150‑1 (MQ=255)                                                                                                                              ggTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAAt                    >  2:279386/1‑151 (MQ=255)                                                                                                                                 tattatTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATgg                  <  2:859092/150‑1 (MQ=255)                                                                                                                                  attattGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGt                 <  1:1036132/150‑1 (MQ=255)                                                                                                                                  attattGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTg                <  1:116884/151‑1 (MQ=255)                                                                                                                                      ttGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACa             >  2:889534/1‑150 (MQ=255)                                                                                                                                       tGCGGGAACGCTCCAGCACCTCCGCTGAGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGCACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGAAATGGTGACAAt           >  2:967490/1‑151 (MQ=255)                                                                                                                                        gCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAAtt          <  2:775948/151‑1 (MQ=255)                                                                                                                                          gggAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGt        <  1:938762/151‑1 (MQ=255)                                                                                                                                          gggAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGt        <  1:502266/151‑1 (MQ=255)                                                                                                                                             aaCGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCtt     <  2:758248/151‑1 (MQ=255)                                                                                                                                              aCGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGGCGAAGCAATGGTGACAATTGTCTTc    <  2:809002/151‑1 (MQ=255)                                                                                                                                               cGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCtt     <  2:981155/149‑1 (MQ=255)                                                                                                                                                gCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCTTc    >  2:188548/1‑149 (MQ=255)                                                                                                                                                gCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCTTCAg  <  2:161605/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCTTCAG  >  NC_000913/4161318‑4161611

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