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%	282.0 / NA	72	G115S (GGC→AGC)	fabR	transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/34);  total (38/34)

GAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCTTCAGTG  >  NC_000913/4161318‑4161613                                                                                                                                                     |                                                                                                                                                    gagCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGc                                                                                                                                                   <  2:511440/151‑1 (MQ=255)     ggTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATTGGTAATAATCTTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCAcc                                                                                                                                                >  1:236464/1‑150 (MQ=255)      gTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTc                                                                                                                                              >  2:500002/1‑151 (MQ=255)          aaTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGc                                                                                                                                           <  2:543424/150‑1 (MQ=255)           aTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTg                                                                                                                                         <  1:121958/151‑1 (MQ=255)             gCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGc                                                                                                                                        >  2:1150505/1‑150 (MQ=255)              cTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGt                                                                                                                                      >  2:415592/1‑151 (MQ=255)               tACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGtt                                                                                                                                     >  1:611057/1‑151 (MQ=255)                        tCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTgccgc                                                                                                                            >  1:140736/1‑151 (MQ=255)                         cATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTgccgcc                                                                                                                           <  2:424157/151‑1 (MQ=255)                            gcgcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGtt                                                                                                                        <  2:1096600/151‑1 (MQ=255)                             cgcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTg                                                                                                                       <  1:1074999/151‑1 (MQ=255)                              gcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTgc                                                                                                                      <  2:537316/151‑1 (MQ=255)                                 aGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGcgcg                                                                                                                   >  2:117604/1‑151 (MQ=255)                                  ggCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGcgcg                                                                                                                   >  2:135976/1‑150 (MQ=255)                                     gcgTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGaa                                                                                                               <  2:918924/151‑1 (MQ=255)                                         cAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTc                                                                                                           <  2:1035896/151‑1 (MQ=255)                                         cAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTc                                                                                                           <  1:883558/151‑1 (MQ=255)                                         cAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTc                                                                                                           >  2:358543/1‑151 (MQ=255)                                               aTCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACt                                                                                                     >  2:348898/1‑151 (MQ=255)                                                 cGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTc                                                                                                   <  2:162841/151‑1 (MQ=255)                                                  gCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTc                                                                                                   >  1:1058877/1‑150 (MQ=255)                                                   ccAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCAt                                                                                                 <  1:304054/151‑1 (MQ=255)                                                      aaGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGc                                                                                              >  2:4738/1‑151 (MQ=255)                                                          cGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCgg                                                                                            >  1:873856/1‑149 (MQ=255)                                                                      ccGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCAc                                                                                                      <  2:390808/127‑1 (MQ=255)                                                                      ccGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCAc                                                                                                      >  1:390805/1‑127 (MQ=255)                                                                      ccGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTa                                                                               <  2:799102/150‑1 (MQ=255)                                                                      ccGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTa                                                                               <  2:489234/150‑1 (MQ=255)                                                                      ccGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTAt                                                                              <  1:344721/151‑1 (MQ=255)                                                                            cTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGa                                                                         >  2:240284/1‑150 (MQ=255)                                                                               ggTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTc                                                                     >  1:75804/1‑151 (MQ=255)                                                                               ggTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTc                                                                     <  2:176299/151‑1 (MQ=255)                                                                                 tctcCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCg                                                                    >  1:433874/1‑150 (MQ=255)                                                                                  ctcCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGa                                                                   >  2:172619/1‑150 (MQ=255)                                                                                   tcCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGaa                                                                  >  2:454824/1‑150 (MQ=255)                                                                                    ccACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGaaa                                                                 >  1:126711/1‑150 (MQ=255)                                                                                       caTTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAAcc                                                              >  2:604327/1‑150 (MQ=255)                                                                                        aTTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCat                                                            >  2:251895/1‑151 (MQ=255)                                                                                          ttATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCata                                                           <  1:660617/150‑1 (MQ=255)                                                                                           tATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATg                                                         >  1:257799/1‑151 (MQ=255)                                                                                                  ttCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCgcgtg                                                  >  1:413576/1‑151 (MQ=255)                                                                                                    cATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCgcgtgcg                                                <  1:567342/151‑1 (MQ=255)                                                                                                     aTCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCgcgtgcgt                                               <  2:1058902/151‑1 (MQ=255)                                                                                                      tCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGtt                                              >  2:1042714/1‑151 (MQ=255)                                                                                                       cGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGttt                                             <  1:1022998/151‑1 (MQ=255)                                                                                                       cGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGttt                                             <  2:663768/151‑1 (MQ=255)                                                                                                       aggTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGttt                                             >  1:1068371/2‑151 (MQ=255)                                                                                                          taataaTCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACt                                          >  2:991942/1‑151 (MQ=255)                                                                                                              aatCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTAc                                           <  2:571123/146‑1 (MQ=255)                                                                                                              aatCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTAc                                           >  1:571119/1‑146 (MQ=255)                                                                                                                  cTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAgcgc                                   >  2:1157977/1‑150 (MQ=255)                                                                                                                   tAACGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCa                                  >  2:319604/1‑150 (MQ=255)                                                                                                                    aaCGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCa                                  <  1:1069656/149‑1 (MQ=255)                                                                                                                    aaCGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCa                                  <  2:1003438/149‑1 (MQ=255)                                                                                                                      cGCCTTCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGACAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGcc                              >  1:403572/1‑151 (MQ=255)                                                                                                                           tCCGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGc                         >  1:741800/1‑151 (MQ=255)                                                                                                                            ccGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCa                        <  2:691260/151‑1 (MQ=255)                                                                                                                             cGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGaaa                                                                 >  2:12726/1‑109 (MQ=255)                                                                                                                             cGGTTATTATTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGaaa                                                                 <  1:12726/109‑1 (MQ=255)                                                                                                                                ttattaTTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATgg                    >  2:508415/1‑151 (MQ=255)                                                                                                                                ttattaTTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATg                     <  2:303665/150‑1 (MQ=255)                                                                                                                                 tattatTGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATgg                    <  2:1199878/150‑1 (MQ=255)                                                                                                                                  attattGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTg                  <  2:421142/151‑1 (MQ=255)                                                                                                                                     attGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGAc                <  2:1150016/150‑1 (MQ=255)                                                                                                                                       tGCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACaa              <  2:595144/150‑1 (MQ=255)                                                                                                                                        gCGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAAtt            >  1:339547/1‑151 (MQ=255)                                                                                                                                         cGGGAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTg           <  1:669391/151‑1 (MQ=255)                                                                                                                                          gggAACGCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGt          >  1:506760/1‑151 (MQ=255)                                                                                                                                                gCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCTTc      >  1:97179/1‑149 (MQ=255)                                                                                                                                                gCTCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCTTCAg    <  1:488616/151‑1 (MQ=255)                                                                                                                                                  tCCAGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCTTCAGTg  >  2:1072545/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTATCTGGAACTCGAAAACCATATGCCGCGTGCGTTTACTGAAGCGCAAGCCGAAGCAATGGTGACAATTGTCTTCAGTG  >  NC_000913/4161318‑4161613

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